#7) TER MITWIRKUNG VON L. COLLATZ : H.GÖRTLER-- J. HEINHOLD - K. KLOTTER 
% MARGUERRE - H. NEUBER - L. SCHMETTERER - K. SCHRÖDER . H. SCHUBERT 
W. TOLLMIEN - H.UNGER UND C. WEBER . HERAUSGEGEBEN VON H. HEINRICH, DRESDEN 


BAND4. SEITE 337— 396 HEFT9 SEPTEMBER 1961 
e s Be AUS DEM INHALT: 
E | | | HrA: US POT ASUS Re 


e F. I. Frankl: Untersuchungen auf dem Gebiete der 
1 schallnahen Strömungen / M.M. Stanisic: On the 
E, Aerodynamic Forces Acting on a Delta-Wing in 
> Supersonic Flow with Supersonic Leading Edges / 
R. Oba: Linearized Theory of Supercavitating 
Flow through an Arbitrary Form Hydrofoil / 
”- H. Ehrmann: Schranken für Schwingungsdauer - 
; und Lösung bei der freien ungedämpften Schwin- 
gung / K.-W. Gaede: Über die Lösung linearer 
Integralgleichungen mit einer elektronischen Ana- 
logierechenanlage 


n 


KL, ESEN TERM IT TE UN IGSEEE 
BUT GAHTBSESSERIE/C HYU ENGE 
EINGEGANGENE BÜCHER 


NWAZTE Fin «Ro: mn C Fre 


BER DEMIE-VERLAG GMBH - BERLIN 


Bd. 41 ie Nr. 9 337—396 Berlin, September 1961 


INHALT: 


"Hauptaufsätze Seite 
F. 1. Frankl: Untersuchungen auf dem Gebiete der schall- 
nahen.Strömungen. . = le var ee a 337 


M.M. Stanisic: On the Aerodynamic Forces Acting on a 
Delta-Wing in’ Supersonic Flow with Supersonic Leadingh 


Eee PD a Re me I 342 
R. Oba: Linearized Theory of Supercavitating Flow throug 
an Arbitrary Form Hydroföil . ... 2»: 22... + 354 


H. Ehrmann: Schranken für Schwingungsdauer und Lö- 
sung bei der freien ungedämpften Schwingung . . .. . 364 


K.-W. Gaede: Über die Lösung linearer Integralgleichun- 
gen mit einer elektronischen Analogierechenanlage . . 370 


Kleine Mitteilungen 


A. Charkavorti: Flexure of a Compound Cylindrical 
Beam Composed of Non-Homogeneous Material . . . 390 


Buchbesprechungen . ... 2: 2.2202 n0 0. 391 
Eingegangene Bücher: 1. NR: 0.0 395 
Winchrichten. 0... aegere Mr Ba 396 


Wir bitten, alleManuskriptsendungen direkt andenHerausgeber, 
Prof. Dr.-Ing. H. Heinrich, Dresden A 27, Friedrich-Hegel-Str. 31, 
zu richten. Zu den Arbeiten, die als Hauptaufsätze bestimmt 
sind, erbitten wir auf gesondertem Blatt eine kurze Zusammen- 
fassung des Inhalts, nach Möglichkeit in deutscher, englischer 
und russischer Sprache; falls die Übersetzungen nicht geliefert 
werden können, ist wenigstens die Angabe spezieller Fachaus- 
drücke in den verschiedenen.Sprachen erwünscht. Die Arbeiten 
sollen in klarer Schrift,- möglichst mit Schreibmaschine, weit- 
zeilig und einseitig geschrieben sein und die nötigen Hinweise 
für den Setzer betreffend Schrifttypen (z.B. griechisch, Fraktur), 
Sperrungen o. a. enthalten. Zur Beschleunigung des Drucks 
und zur Vermeidung von Satzfehlern empfiehlt es sich, umfang- 
reiche und unübersichtliche Formelausdrücke durch Einführung 
von Abkürzungen zu vermeiden. Bilder sollen als Tusch- 
zeichnungen auf Transparentpapier (Beschriftung mit Bleistift) 
oder als saubere Bleistiftskizzen ausgeführt sein. Bildunter- 
schriften sind am Schluß des Textes anzufügen. Sollte die Arbeit 
bereits an anderer Stelle verbreitet sein (als Dissertation, For- 
schungsbericht, Manuskriptdruck o. ä.), so ist dies auf der ersten 
Textseite in einer Fußnote anzugeben. Für die Zusammenstellung 
der zitierten Literatur, amSchluß der Arbeit, bitten wir, sich einer 
einheitlichen Notation entsprechend folgendem Müster zu bedie- 
nen, z.B. W.Schmeidler, Über die Wärmespannungen in ei- 
nem Körper, ZAMM 28 (1948), S. 54—59 oder G. Bürger- 
meister und H.Steup, Stabilitätstheorie I, 1. Aufl., Berlin 1957, 
Akademie-Verlag, S. 142—147. Die Autoren erhalten von den 
Hauptaufsätzen 75, von den Kleinen Mitteilungen 25 Sonder- 
drucke ohne Berechnung, darüber hinaus bis zu 250 Sonder- 
drucke gegen Berechnung. 


Der Verlag behält sich für alle Beiträge das Recht der Verviel- 
fältigung, Verbreitung und Übersetzung vor. 


EN EEE TE EIER EEE ET HEERES: USREEFEREEN 


Ab Band 8 lautet der Reihentitel — e 
SCHRIFTENREIHE 
DERINSTITUTE 
FÜR MATHEMATIK 


BEI DER DEUTSCHEN AKADEMIE 
DER WISSENSCHAFTEN ZU BERLIN 


und wird herausgegeben x 


v4 


u %4 


von u: 
Prof. Dr. Heinrich Grell, Prof. Dr. Josef Naas, Prof. Dr. on 


Achilles Papapetrou, Prof. Dr. Hans Reichardt, Dr. 
Willi Rinow, Prof, Dr. Erhard Schmidt und Dr. 


herausgegeben von Prof. Dr. J. NAAS und Prof. Dr. 
KURT SCHRÖDER 

1957. VIIL, 317 Seiten — 22 Abbildungen — 9 Tafeln — gr. 8 — 
DM 38,— (vergriffen) 


Lösung des allgemeinen Randwertproblems für eindimen- 
sionale Wellen bei harmonischem Zeitgesetz 
Von Dr. KARL BORKMANN und SIEGFRIED 
OBERLÄNDER 
1955. 99 Seiten — gr. 8° — DM 12,— (vergriffen) 
Von A. W. POGORELOW 
(Übersetzung aus dem Russischen) j 


Wissenschaftliche Redaktion: Prof. Dr. J. NAAS 
1956. 79 Seiten — gr. 8° — DM 5,50 
Bemerkungen über die Stabilitätsuntersuchungen der 
Wirbelstraßen 
Von Prof. Dr. Bl. DOLAPTSCHIEW 
1957. 28 Seiten — gr. 8° — DM 3,80 
Die Verbiegung konvexer Flächen 
Von A. W. POGORELOW 
(Übersetzung aus dem Russischen) 
Wissenschaftliche Redaktion: Prof. Dr. E. REMBS 
1957. 135 Seiten — 26 Abbildungen — gr. 8° — DM 18,50 . 
in die Theorie stationärer Zufallsfunktionen 
Von A.M. JAGLOM 
rsetzung aus dem Russischen) 
Deutsche Übersetzung unter wissenschaftlicher Redak- 


tion von Dr. HERBERT GOERING 
1959. VIII, 177 Seiten — 8 Abbildungen — gr. 8° — DM 24,— 


Tabellen von Exponentialfunktionen und -integralen zur 
Anwendung auf Gebieten der Thermodynamik, Halbleiter- 
theorie und Gaskinetik 


Von SIEGFRIED OBERLÄNDER 
1959. VIII, 8 Seiten Text — 142 Tabellenseiten — 1 Abbildung — 
gr. 8° — Ganzleinen DM 35,— (vergriffen) 


Konstruktion ganzer, rationaler und reeller Ordinalzah- 
len und die diskontinuierliche Struktur der transfi» 
niten reellen Zahlenräume . 


Von DIETER KLAUA 
1961. 141 Seiten — gr. 8° — DM 28,— 
Qualitative Methoden beim n-Körperproblem der Him- 
melsmechanik 
Von F, CHILMI 
(Übersetzung aus dem Russischen; vom Autor durch- 
gesehene und ergänzte Übersetzung. In deutscher 
Sprache herausgegeben von J. O. Fleckenstein) 
1961. 116 Seiten — 4 Abbildungen — 1 Tabelle — gr. 8° — 
M 23, — 
Über eine Bewertungstheorie der Algebren und ihre Be- 
deutung für die Arithmetik 
Von HERBERT BENZ 
1961. 149 Seiten — 5 Abbildungen — gr. 8° — DM 29,50 


Fordern Sie bitte unseren Fachkataiog MATHEMATIK an, der 
Ihnen kostenlos geliefert wird. . 


Bestellungen durch eine Buchhandlung erbeten 


AKADEMIE - VERLAG.BERLIN 


LED vie 


F Strömung enden. 


-3 


ne We 2 ZAMM 4 (1961) Heft 9, Seite 33731 


„ur u 
1. ur en u > Beier > schallnahen Strömungen”) 
£ Von F.1. Bass ; 


Um: Be eines le bei een im elchen 
2. Die Frage der transsonischen Strömungen mit Verdichtungsstößen, die innerhalb der 


a Direkte Randwertaufgaben 


i _ Unter direkter Aufgabe verstehen wir die Aufgabe der Bestimmung der transsonischen 
"Strömung i in einer Düse mit gegebenen Wänden bzw. die Aufgabe der Bestimmung der Umströ- 


mung eines gegebenen Profils. 


Es wird dabei die Störungsmethode weiter ausgearbeitet, die erstmalig von A. A. NIKoL- 
sK1J [1] vorgeschlagen worden ist. 
“ Sei CAEBD der ‚Hodograph der Wände L,, L, der Ausgangsdüse L, O das Zentrum der Strö- 


— mung (d.h. der Punkt, in dem sich Schallinie und Stromlinie unter rechtem Winkel schneiden), OC, 


OD Charakteristiken, A, B Punkte der Schallinie; E entspricht dem unendlich fernen Punkt 
_ stromaufwärts (Bild 1). 


bz 
Bild1 


Die Stromfunktion der Ausgangsströmung y genügt dabei der Gleichung 


DV RAR STERE  ISEAD 
wobei # der Neigungswinkel der Geschwindigkeit und o, K!) die bekannten, von S. A. TscHA- 
PLYGIN eingeführten Funktionen des Geschwindigkeitsbetrages w sind. 

Die variierte Düse L möge sich von L wenig unterscheiden; der Abstand zweier entspre- 
chender Punkte von L und L sei ön; dieser Abstand ist in der Richtung von L Dach: L gemeint, 
normal zu L und, wenn man stromab sieht, von rechts nach links. 


*) Der Beitrag stellt den Inhalt eines Vortrages dar, der auf Einladung der Tagungsleitung aufder wissen- 
schaftlichen Jahrestagung 1960 der GAMMin Freiberg/Sa. gehalten werden sollte. Der Verfasser ist leider am 
7.4.1961 in Naltschik (USSR) verstorben. 


1) o = 0 entspricht der kritischen Geschwindigkeit. 
i 23 


+. 
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Weiter führen wir die Funktion » ein, die aus der Stromfunktion y durch Berührungs- 
transformation entsteht: 


o=y— „(uy— 2) er 2 2 


(x, y — die kartesischen Koordinaten in der Strömungsebene, u, v — die Geschwindigkeits- 
komponenten, g — die entsprechende Dichte, 0, — die Dichte im Ruhezustand). 
Die transformierte Stromfunktion » genügt der Gleichung 
1 dK 
ED en u RE EN: 
K d o Vg 0 ( ) 
Dann ist die Variation öw auf dem Hodographen von L in erster Näherung bekannt, und 
zwar haben wir 


Kwoo + os — 


RER Fi wen +A80) ri = — MO HF On - nee . (4). 
0 “ 


Dabei sind ög, , die Variationen der gewöhnlichen Stromfunktion y an den Wänden und ö9, der 
Winkel, um den man die Düse zusätzlich drehen muß, damit sich der Neigungswinkel der Ge- 
schwindigkeit im Zentrum der Strömung nicht ändert; im Falle einer symmetrischen Düse ist 
dieser Winkel natürlich gleich null. Was die Variationen der Stromfunktion ög,,, betrifft, so hat 
natürlich nur ihre Differenz öq = ög, — 69, d. h. die Variation der Durchflußmenge, einen physi- 
kalischen Sinn. 

Daß die Randwertaufgabe (4) für die Gleichung (3) nicht mehr als eine Lösung haben kann, 
hat vor einigen Jahren C. S. MoRAwETZz [2] im Zusammenhang mit einer anderen gasdynamischen 
Frage bewiesen; den Existenzbeweis für eine schwache Lösung erbrachte unlängst der Verfasser 
[3] mit einer Methode, mit der vorher MoRAwErz das entsprechende Problem für die Gleichung 
(1) gelöst hat [4]?). 

Beim Existenz- und Unitätsbeweis spielt aber der Wert der zusätzlichen Durchflußmenge 
keine Rolle; dasselbe gilt für den Winkel ö9,. Infolgedessen vermutete ich eine Zeit lang, daß sich 
die Größe ögq wirklich nicht eindeutig bestimmen lasse und daß sie sich in gewissen engen Grenzen 
ändern könne; dies widerspräche aber sowohl der elementaren hydraulischen Theorie der Laval- 
düse, als auch den Versuchsdaten (wie A. A. NIKOLSKIJ und G. I. TaGavow unlängst auf der 
Mechanikertagung der Sowjetunion mitteilten, sind sie von der Vermutung der Nichtunität der 
transsonischen Strömung in der Lavaldüse abgekommen, nachdem die von ihnen zur Prüfung 
dieser Hypothese angestellten Versuche mißglückt waren). 

Gegenwärtig halte ich es für wahrscheinlich, daß sich die Größen ög, 68, durch Regulari- 
tätsbedingungen eindeutig bestimmen lassen. Man muß nämlich fordern, daß die Größe dw,/o 
in den Punkten A, B, O endlich bleibe: im Gegenfall würden die entsprechenden Punkte der 
Strömungsebene ins Unendliche rücken. Diese Forderung führt zu der Gleichung 


2 = = w(ön + nö0) (Differentiation = längs L) ee RENTE 
in den Punkten A, B (ausführlicher siehe [7], [23]). Es ist sehr wahrscheinlich, daß diese Forderun- 
gen voneinander abhängig sind: im Gegenfall kämen wir zu dem Schluß, daß eine stationäre 
stoßfreie transsonische Strömung nicht in einer jeden beliebigen Düse möglich wäre — sogar, 
wenn sie sich beliebig wenig von einer solchen Düse unterscheidet, in der eine derartige Strömung 
möglich ist. Aber eine solche Annahme widerspricht, soweit mir bekannt ist, der Erfahrung. 

Eine zweite unabhängige Gleichung für öq, 60, bekommen wir auf Grund der Forderung 
der Endlichkeit der Größe öw,/o im Punkte O, der dem Zentrum der Strömung entspricht. 

Eine mathematisch strenge Lösung dieser Frage ist bis heute noch nicht gelungen. 

Im Falle der direkten Aufgabe der Umströmung eines Profils mit Schallgeschwindigkeit ist 
die Lage eine ähnliche. Wenn man die Störungsmethode von NIKOLSKIJ anwendet, so ist ein un- 
bestimmter Parameter vorhanden, der mit dem asymptotischen Verhalten der Strömungs- 
geschwindigkeit im Unendlichen zusammenhängt (ausführlicher siehe [7]). Wir benutzen ihn dazu, 
um die Forderungen der Endlichkeit der Größe öw,/o in beiden „Schallpunkten‘“ des Hodogra- 
phen des nicht variierten Profils zu erfüllen. Aus Gründen, die dem Fall der Düse analog sind, 
hängen diese Bedingungen wahrscheinlich voneinander ab. 

Im Falle der Umströmung mit Überschallgeschwindigkeit, wenn sich zwischen Kopfwelle 
und Profil ein Unterschallbereich bildet, entstehen vier Regularitätsbedingungen, und zwar die 
Forderungen der Endlichkeit der Größe öw,/o in den Schnittpunkten des Hodographen des nicht 


öde = 


Bu 5 
2) Eine Näherungs hode zur Lös N a ne b Er 
er a Hr ne a zu1 Lösung der Aufgabe habe ich in [5] vorgeschlagen; im einfachsten Fall, 
sich dıe Auigabe auf das Trıcomı-Problem reduziert, gelingt die Lösung durch Quadraturen [6]. 
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variierten Profils sowie der Stoßpolare mit der Schallinie. Um die Übereinstimmung mit den Ver- 
suchsdaten zu erhalten, welche für Existenz und Unität der Lösung sprechen, muß man anneh- 
men, daß von diesen Bedingungen nur zwei unabhängig sind (eine auf dem umströmten Profi, 


eine auf der Stoßpolare). In der Tat, da die Randbedingung für dw auf der Stoßpolare eine 


Differentialbedingung zweiter Ordnung ist (siehe [7]), stehen uns offenbar zwei Parameter zur 
Verfügung, die aus den oben genannten Bedingungen bestimmt werden können. Wir erinnern 
im Zusammenhang damit daran, daß O. M. BELOTsERKOVsK1J [8] für den Fall einer sehr hohen 
Überschallgeschwindigkeit durch Verwendung einer ähnlichen Regulatitätsforderung eine einzige 
Lösung bekommen hat, die gut mit den Versuchsdaten übereinstimmt; BELOTSERKOVSKIJ arbei- 
tete direkt in der Strömungsebene und berücksichtigte die Wirbel, die hinter der Kopfwelle 
entstehen. 

Zum Schluß erinnern wir noch daran, daß in den einfachsten Fällen die Unität der Lösung 
der oben genannten Randwertaufgaben schon vor langem bewiesen worden ist. Es sind dies der 
Fall einer Düse, deren konvergierender Teil geradlinige Wände hat, und der Fall der symmetri- 
schen Umströmung eines rhombischen Profils. Beide Aufgaben führen auf das TRıcomı-Problem 
für die gewöhnliche Stromfunktion, was im ersten Fall der Verfasser [9], im zweiten GUDERLEY 
und YosHIHARA [10] und gleichzeitig L. V. OvsJsanNıKov gezeigt haben. Der zweite Fall ist der 
Grenzfall der Aufgabe der Umströmung eines rhombischen Profils bei Überschallgeschwindigkeit 
mit abgelöster Kopfwelle, für die die Unität ebenfalls der Verfasser gezeigt hat [9]. 


2. Transsonische Strömungen mit Verdiehtungsstößen, die innerhalb der Strömung enden 
Solche Strömungen werden, wie bekannt, beobachtet, wenn sich ein Körper in der Luft mit 
einer Geschwindigkeit bewegt, die nur wenig kleiner ist als die Schallgeschwindigkeit. Die Rela- 
tivgeschwindigkeit ist in einem solchen Falle nur in einem kleinen Gebiet höher als die Schall- 
geschwindigkeit, und der Überschallbereich endet meist mit einem Verdichtungsstoß. 
Ein theoretisches Beispiel einer ebenen Strömung dieser Art ist erstmalig vom Verfasser 
konstruiert worden [11]. Man geht dabei davon aus, daß die Stromfunktion in der Nähe der Schall- 
geschwindigkeit näherungsweise der TrıcomI-Gleichung genügt: 


TE A RR 2% (6), 


wobei n die seinerzeit vom Verfasser [12] eingeführte Funktion der Geschwindigkeit ist. Im ge- 
gebenen Fall ist die Stromfunktion ein Binom: 


DE Rn. 7 ee (7. 


Die Stromlinien und ihre Hodographen sind in Bild 2 gegeben. Der Verdichtungsstoß ist in diesem 
Falle ein gerader. Die vordere Hälfte der Stromlinie, die durch das Ende des Verdichtungsstoßes 
geht, fällt mit der Schallinie zusammen. 
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| 
| 
Stromlinien 
| 


IR const. 


RN Verdichtungsstoß 


af DE 
Grenzlinien 
Bild 2 


Ein Mangel dieser Lösung besteht darin, daß sich im Überschallbereich ein kleiner Streifen 
bildet, in dem die Geschwindigkeit dreideutig wird; dieser Streifen ist von zwei Linien begrenzt, 
die am Ende des Verdichtungsstoßes enden und ihn dort berühren. 

Infolgedessen ist diese Lösung von P. GERMAIN angegriffen worden [13]. GERMAIN weist 
darauf hin, daß der Exponent » im Ausdruck einer homogenen Lösung der Trıcom1-Gleichung 


pe ee ee I) 
zwischen 0 und 5/, liegen müßte, um diese lokale Dreideutigkeit zu vermeiden. Es gelang ihm aber 


nicht, auf solche Weise zu erreichen, daß der Überschallbereich vollständig zwischen der vorderen 
23% 
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Hälfte der durch das Ende des Verdichtungsstoßes gehenden Stromlinie und dem Verdichtungs- 
stoße liege, wie es faktisch beobachtet wird. Wahrscheinlich ist das prinzipiell unmöglich. In 
einem vom Verfasser konstruierten Beispiel [14] (mit » = 2/3) liegt der Unterschallbereich vollstän- 
dig zwischen dem Verdichtungsstoße und der hinteren Hälfte der durch das Ende des Verdich- 
tungsstoßes gehenden Stromlinie; wir haben es also hier mit einem ganz anderen Strömungs- 
typus zu tun als dem, der bei Umströmung von Profilen bei hohen Unterschallgeschwindigkeiten 
beobachtet wird. 

Die genannten Beispiele [11, 14], in denen der Verdichtungsstoß ein gerader ist, hat I. Bıs- 
Bossunow [15], [16] verallgemeinert und Beispiele konstruiert, in denen der Verdiehtungsstoß 
gekrümmt ist. 

Im Beispiel [11] kann man die lokale Dreideutigkeit der Geschwindigkeit loswerden, wenn 
man eine zusätzliche Sprunglinie einführt. Von den vier aus den Erhaltungssätzen hervorgehen- 


den Stoßbedingungen (siehe [13]) 


N, do’ es do!’ Bach —n' —n" 7 ade 1 7 n2 7 r 
sind die ersten drei genau erfüllt. Statt der vierten (der „‚Stoßpolarenbedingung‘‘) haben wir aber 
nach I. BrsBossunow und F. ABuTALIEW [17] 


tt _. HrN2 
m ai, da — —0,9882620 . . 2.2.2 2.3 (10), 
a 
so daß der Fehler im ganzen 1,2%, ausmacht. Man bemerkt aber, daß diese zusätzliche Sprung- 
linie kein Verdichtungs-, sondern ein Verdünnungsstoß wäre, was, streng genommen, dem zwei- 
ten Hauptsatz der Thermodynamik widerspricht. 

Der Ausweg aus diesen Widersprüchen besteht wahrscheinlich darin, daß Strömungen vom 
betrachteten Typus nicht streng stationär sein können. Die theoretischen stationären Strömungen 
[11], [15] mit lokaler Dreideutigkeit des Geschwindigkeitsfelds sind aber zweifellos gute Annähe- 
rungen an die Wirklichkeit. 

7 


Bild 3 


Von diesem Standpunkt ausgehend, hat der Verfasser eine Reihe von Randwertaufgaben 
in der Hodographenebene formuliert, die auf Umströmungen von Profilen mit lokalen Überschall- 
bereichen führen, welche mit Verdichtungsstößen enden. 

ur LM dieser Aufgaben entspricht der symmetrischen Umströmung eines symmetri- 
schen Profils mit geradem Verdichtungsstoß [18], [24]. Die I i e i 

g58 ‚[24]. Die Randbedingungen sind folge 
(siehe Bild 3): - er 


El) auf ABC a er ee 
y=1(9) "auf! Arm 10 Se 
Be / 
y=oT-U2sin Sara e=1/9 +(s— sp), Se En nr, 
> 3 = (13), 
s—-sg=ocost, d=osint 
v(0,,n) = yv(dy, — N), voldn, n) = 0 auf DE Were 


®) EF ist eine Charakteristik, ED= DC. 


Be = 
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E fi 


Die Bedingung (13) bedeutet die Eindeutigkeit des Geschwindigkeitsfeldes im Unendlichen, 
dem in der Hodographenebene der Punkt B entspricht. Die Bedingungen (14) charakterisieren 
den geraden Verdichtungsstoß. 


Unter gewissen Einschränkungen haben A. W. Bırsapsz [19] und J. W. Drwingrar [20] 
die Existenz und Unität der Lösung dieser Randwertaufgabe bewiesen. 


Die Aufgabe läßt sich auf die Fälle von Zirkulationsströmungen und von gekrümmten Ver- 
dichtungsstößen verallgemeinern [18], [21]. 


Wenn man auf diese Aufgabe noch die Störungsmethode von A. A. NIKOLSKIJ anwendet, so 

kann man auch eine entsprechende direkte Randwertaufgabe formulieren, d. h. eine solche, bei 

der das umströmte Profil in der Strömungsebene gegeben ist [22]; in den vorher formulierten 

Aufgaben [18], [21] sind die Hodographen eines Teils des umströmten Profils, sowie des Ver- 
dichtungsstoßes gegeben, und das umströmte Profil selbst wird gesucht. 

Existenz und Unität der Lösungen der erwähnten Randwertaufgaben sind einstweilen nur- 


in Spezialfällen bewiesen; auch die Ausarbeitung von effektiven Lösungsmethoden erfordert noch 
weitere Untersuchungen. Br: 
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On the Aerodynamic Forces Acting on a Delta-Wing 
in Supersonic Flow with Supersonic Leading Edges 


By Mıromir M. STanısıd 


Es wird eine Methode zur Bestimmung der aerodynamischen Kräfte dargelegt, welche auf einen Delta- 
flügel mit Überschall-Vorderkante in Überschallströmung wirken. . Bee 

Das Wesentliche bei der Methode ist, die Lösung für eine stationäre Überschallströmung bei beliebiger 
Verwindung und Geometrie des Deltaflügels zu finden. Die Lösung für die stationäre Strömung im Falle 
eines ebenen Tragflügels erscheint dann als Spezialfall. Bine funktionale Beziehung zwischen stationärer 
und nichtstationärer Strömung wird mit Hilfe der Transformationsmethoden der Operatorenrechnung ge- 
wonnen. Die Lösung für die nichtstationäre Überschallströmung wird in Form eines geschlossenen I ntegrals 
angegeben, wobei der Integrand die Lösung für den stationären Fall enthält und keinerlei pragmatische 
Schwierigkeiten auftreten. H 

Die Methode kann ohne Schwierigkeiten auf andere Flügelformen ausgedehnt werden; dabei kann die 
Strömung in bezug auf die beiden Flügel symmetrisch oder — wie im Schiebeflug — unsymmetrisch sein. 


A method is presented in order to find the aerodynamic forces acting on a delta-wing with supersonic 
leading edges in supersonic flow. . 

The basic postulations of the method consist of the determination of the solution for a steady supersonie 
flow passing a delta-wing of arbitrary warping and chamber geometry. The solution for the steady flow for 
the flat lifting wing appears as a special case of this study. 

A functional relation between steady and unsteady solution is obtained using the transform techniques 
of operational caleulus. The solution for the unsteady supersonic flow is presented in the closed integral 
form, in which the integrand contains the solution of the steady flow, and no diffieulties of a pragmatical 
nature arise. 

The method can be extended without any difficulty to the wings of other shapes, in either symmetrical 
or yawed flow. 


B pa6oTe HalaraeTcaA MeTon MIA ONpeMeleHHA Aa3apONMHAMHYeCKHX CHI, KOTOPbIe Meii- 
CTBYEPT B CBEPX3BYKOBOM HOTOKE HA CTPEJOBHAHLIE KPbLWIbBA C CBEPX3BYKOBOH HepenHeli KPoM- 
Koi. 

CyINecTBo MeToAa 3aKJIoYAaeTCca B TOM, YTOÖBI HAÄTH PellleHHe AJIA CTAUHOHAPHOTO CBEPX3BYKO- 
BOTO NHOTOKA NPH AI60H 3AKPYTKe MU TeOMETPHH CTPENOBHAHBIX KpblibeB. Permenme AA CTa- 
UINMOHAPHOTO MOTOKA B CIIyYae INIOCKOTO KPbLAlA HOJIYYAeTCA KaK CHEINMAAIbHBIa cayyai. C 
IIOMOINbIO METONOB IIPe0o0pa330BaHNA ONEePATOPHOTO HCYHUCHEHNHA VCTAHABAHBACTCH DYHKUH- 
HAJIBHOe COOTHOIIEHHE MEKAY CTAIMOHAPHBIM H HEeCTAUHOHAPHBIM IIOTOKAmM. Peimiennme ua 
HEeCTalMoOHAPHOTO CBEPX3BYKOBOTO MOTOKA MaeTCcA B BUNE HHTETPAaNa, Te HONBIHTETPANBHOE 
BbIpasKeHNMe CoMep>KUT PelleHme MIA CTAIIHOHAPHOTO NOTORA. 


ITOT METON MOHKET ÖbITb TIPHMEHEH Öe3 TpyAa K APYTUM POoPMaM KpblllbeB; IIPH TOM IIOTOR 
IIO OTHOINIeHWIO K OÖ00HM KPbIJIBAM MOHteT ÖbITb CHMMETPHYHEIM HAM KAK IIPM CKOAB3AIIEM 
IIOJIETe HECHMMETPHUYHBIM. 


1. Introduetion 


Recently'), a method has been developed for the determination of aerodynamic forces 
acting on a delta-wing with subsonic leading edges, in supersonic flow [1], [2]. The method of 
acceleration potential [3] was applied and the problem was reduced to the solution of a singular 
integral equation, which was solved employing the method of perturbation of the domain of 
integration [4]. 

The problem of the determination of the aerodynamic forces acting on a delta-wing with 
supersonic edges has received attention in the development of modern transport jet airplanes [5]. 
With this in mind, a study of the aerodynamic response of a delta-wing with supersonic leading 
edgesin supersonic flow has been made by the author of this paper, upon the suggestion of R. 
SK00G, Senior Group Engineer of the Boeing Airplane Company. 

The solution to the problem of the determination of the aerodynamic forces acting on a 
delta-wing, with supersonic leading edges in the case of steady loading was presented by Busr- 
MANN [6]. BuSEMANN’s solution was extended by Gau [7] and Kovaneva [8] in order to obtain a 
solution for the aerodynamic forces acting on a rectangular wing. KovALEva’s solution was re- 
strieted, however, by the assumption that downwash distribution over the lifting surface is a 
pure time function; this simplified the problem very much mathematically, since the actual con- 
dition required that the downwash distribution be a function of both time and space coordinates. 


However, KovAuevA developed a method for the extension of the steady to unsteady load- 
ing, based on the BUSEMANN steady solution. This extension requires the existence of derivatives 
in the plane direction of the steady-loading function. In the case of a rectangular wing, the ex- 


!) The number in squared brackets refer to the bibliography at the end of this paper. 
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tension gives rise to no difficulty since BusEMANn®’s solution [6, Egq. 15] is a function of plane coordi- 
nates. Moreover, in the case of a triangular wing the BuszmaAnn solution [6, Eq. 17] is a constant 
and therefore, the plane derivatives vanish everywhere. Hence, mathematical diffieulties appear 
but they can be avoided using the ö-DirAc function. 


In that which follows, the solution for steady loading will be obtained in the form of Besser, ' 
functions associated with LEGENDRE Polynomials, as is expected in Gas dynamics, and also as 
suggested by MaGnARADZ [9] and later extended by KovALzvA’s method to unsteady loading. 


2. Formulation of the Boundary Problem 


Let (z, y, z) be rectangular Cartesian coordinates in the field of flow, as shown in Fig. 1, 
with z-axis in the direction of the undisturbed velocity of magnitude U. Furthermore, denote the 
Mach angle by y, so that y is the acute angle for which 


ME ee d), 


XYPlane 


Fig. 1. Geometry of the Delta-Wing with Supersonic Leading Edges 


The partial differential equation describing the disturbance of the flow in the presence of a 
solid body, is given by 


1 
EM’ —1) DE, U, 2, 2) D,,y® U, 2, t) m D,&; U, 2, f) Fu a? [2 U D,&; U, 2, D) nr D,,& U, 2, n] = 0 


2), 
where ®(, y, z, f) is the velocity potential function. 
The force acting on the wing is given by 
h mx 
Be a, REN ZEdHUDEEH Ey 0.2. .:0) 
0 —mz 


where o is the undisturbed density. DER 
The boundary conditions will first be specified. The lifting surface at any time {is given by 


DE en snaleolson Hin Sal); 


where g(z, y) is a regular function everywhere on the surface under consideration and ® is the 
frequency of vibration of the wing. Then, the condition that the normal component of velocity 
on the lifting surface be zero leads to [1]. 
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However, Egs. (5). (6) and (7) lead to vi 


Ur; Z)|.-0 = eF2[U g,(2, y) — io g(&, y)] ; 


a Ki oc 1 a, 
a 


‘ 
= rer. 93e ( ıa<h ) a 
Vz |yl>|m«! s 
Eq. (3) becomes 
0=— 20 Ve [ Mt en) rend + Wenn] aan... | 
0 —mx j iu 
Eq. (10) can be written as aut | 
vi2,9 a = Fr Hr - we (13), 
where 
a )=U0Ge N) laden) 5 Sn wre ee (14). 
Note that g(z, y) is a regular function, and therefore has a TAYLor series expansion in the form >a 
N PT (15). 
je, y) - Due +u 2] a he ) 
Or 
rs ir "0, 0) , 0"f(0, 0) 16 
I, y) = Pa: [2 () T Da y dyr FE a BET ( ) 
Therefore, cs a 
BESSER TE "0, 0) 00, 0) ARE r=012..N 1. 
Ka = nl)” ori r ey ee (17) 
Let 
mt a8 eis ee a (18) 
Thenij=s-+r. Hence, e 
s + r\ 0°f(0, 0) "f(0, ek 
fa, y) = BPN- Era | i 5 ae (19). 


Evidently, 


I, y) -2 2 an Kl; 1er. che Te N (20), 


where a,, are coefficjents depending on the constants U and » and given by 
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Hence, 


2 1.7.9°7(0,.0)°07f(0;:0) 
Tg a (23). 


Eq. (20) can be written in the form 


Denote: kml, 
mir: 7 A uU; EDER EN ee (25), 
where 
eg). \ 
a EEE EEE Se (26) 


Then Egs. (8), (10), (11), (25) and (26) lead to 
— an a) H+ Fun 92 2) + Pu Ye» 1) + Ya, a)=0 ... NM 


and 
v, : r 2-0 = —— Bz, s 2 | 1 jW ” 28 
‚(1 Yı> Zu) ]z,=0 Au N yı) — mus, sy Smuz, > 
Pr,y,2)=0; ( <= AR Er N (29) 
; Y lo, 2)=0; Il>|mua| 
with 
k 
N ee (30) 


Eqs. (27), (28) and (29) represent the boundary value problem under consideration. 


3. Method of Solution 


In order to solve the boundary value problem, a steady loading function *(x,, y,, 2) 
must be found characterized by: 


= — Fa + rd er)... .- (31) 
Pa Yo Zu)la-0 = > era fa, Yı); ( lee (32), 
u —muz sy zmu%, 
Fe 3 In En ) Be re dr (33). 
FE Ye, 2) = 0; \yl>Imuzı| 


Let this problem be referred to as the first auxiliary problem. Clearly, Y*(x,, y,, 2) will be deter- 
mined as soon as the solution is obtained for the second auxiliary problem, given by 


— Ya + Fr dt Fr m) .» > (34), 
0o<m,<Ah Ei 
Yy z » y»» lz Soll; RER AR S RE: Sr 
0, (& Yu Zu)la=0 ee Mi (35) 
es ern | te RR art (36). 
Porn Yu 2) = 0; yl>|maz| 
Note that Egq. (32) can be written as 
1 1 
Pl Yo Zu) =0 = Au | ler fa, Dh. dr :-: 20. (37). 
0 


Hence, if Y$(x,, y, 2) is found to satisfy Egs. (34), (35) and (36), then the solution to the 
first auxiliary problem, Y*(x,, y,, z) at any &,, is simply given by 


ia f a 
Pr, Y Z,)|2,=o = Au | PEr Yo Z)|2,=0 je Pa fi, — &, Yıllz, U. = 138). 
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Integrating by parts, Eq. (38) can be written as 


Pla lan = Zu File Dr Alan LTE I 


+ Pr Ye lan er I rc) BI. 
0 


Or, since the free term vanishes by virtue of the boundaries, then 
7 
1 E at —_ 
P* (2, Yo Zu)l=0o = Au | Perle Yo la — Er y)di - > (AO). 
0 


It must be mentioned that as soon as the steady loading function Y*(x,, Y,, 2) is obtained 
than the unsteady loading function (x, Yı, 2) will be obtained by a procedure to be shown later. 


0-7) . 


( 
Ah 


Fig. 2, Transformation of Delta-Wing in Rectangular Form 


Yı 
mn=e! — = 
murx, 


Firstly, 7, Yı, 2) must be determined. For a wing of rectangular form, GaLın and 
KovauevA adopted Busemann’s solution [6, Eq. 15] for which the derivative Pf, (2,, Yı, 2) neces- 
sary for the determination of steady loading function Y*(x,, y,, 2) exists. However, in the case 
of a triangular wing the pressure difference is a constant [6, Eq. 17]. Recognition of this fact 
is necessary for the construction of the solution in this paper; however, it imposes mathematical 
difficulties in the solution of Eq. (40), since 7,1, Yı, 2) vanishes. This difficulty can be sur- 
mounted by using the ö-DirAc function. But this is not necessary, since Y(z,, y,, 2) can be 
obtained by direct solution of the second auxiliary problem. Note that standardized product 
solution applied to the second auxiliary problem leads to difficulty since the leading edges 
are presented in functional form, namely y, = + m ux,. For such purpose the delta-wing will 
be mapped into a rectangular form and the flow must be subjected to the corresponding trans- 
formation. 


Denote 
yı 


u=8; mr zb ner ne SOEBEN 


This transformation transforms the delta-wing into the rectangular form shown in Fig. 2. 
Hence, Egs. (34) and (41) lead to: 
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Let | | 
I a (43). 
Eqs. (42) and (43) lead to 
ka\2 1 
On ’ >= U, ’ he PSP a a er 
EI Ind = (7) Bd (44) 
and 
1 — (man) ana) mn... (45), 
where k# is a constant to be determined later. Denote 
MEI. TERM ZU a u ee (46). 
Then Egs. (44) and (46) lead to 
Last ÄEN ans ZulC) a0 Me (47). 
And 
Inw1 - er 
Xu) + I) + zi] U(E) — Oe MR Aa (48). 
Solution of Egq. (47) is given by 
S 
ZI CE Sn u dA, ren rk (49), 
s=0 
where C{"®) and C{"®) are integration constants. Egq. (48) through substitution 
Kl) 1ED, (Et bar ten, 2 200 ee (50) 
reduces to 
ESEL BEI AH (51) 
with 
le 
a=- |, — m) | IE WERT SE (52). 


Hence, Egq. (51) is a BEsseEL equation the solution to which is 


Ss 
LIEZEN DEP ATTNE (as ehın Bee (53). 
3s=0 
Since at £ —0, 9,(£) is finite, then AR — 0. Hence, Eq. (50) can be written as 
5 = 
X,(E) = & Ans VE 1 PERLE AR 3 Bela Me u RT Be Se a a a (54). 
s=0 
Moreover, Eq. (45) must be solved. The LEGENDRE differential equation is readily recognized. 
Let 
ET I Er (55) 
Then Egs. (45) and (55) reduce to 
*? * Iym* * kr > * 
A—n*?) Kn,nen{n*) —2n* Kuna?) + m) EU) Da BESSERE 56). 
Putting 
* ,\2 
(2) le 1,2 REN (57), 
mu 


Eg. (56) is a LEGENDReE differential equation. This equation for integral values of n has a polyno- 
mial solution in 7*, namely 


ET Re car anderer nn ei (58), 
where 
CI en- NiyreZan 
Be N (59) 


and N = n/2 for n even, and N = (n— 1)/2 for n odd. Edq. (58) is valued for 
Aa Er > (60). 


This solution converges for — (1/m u) - = ne 
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where b,, are roots of the BesseL function under consideration. However, the boundary. condi- 
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Multiplying both sides of Eq. (66) by VE Jo, (ms 9 P,(m un) and integrating over the wing sur- 
face, one obtains 


EN 
Ah mu ” 
FT IE Oo.) Palm de an 
- i 
Ar 
NIS Ah mu 
= 2 Bi CR J S EI Ensd) To, Ans) Pu(lm un) Pu(m un) dEdn (67). 
Note that ; 
Ah E 
| EI (ns) Jokes E) de = |... (68) 
a EEE 0,5 0a 
and 
I 
my Yor an 
P,(m un) Pu(m un) dy = R) Re (69). 
ie VD ee nn | 
mju | 
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Moreover, on the wing surface for which £ —0, W# (& gi is a constant, say K, then 
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Without loss of generality one can write 
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Therefore, 
DI a a 
Thus, the solution P%(x,, Yı, 2) at z, —0 is given by 
INDEIES: 
Pia = I I 0 al Pula) 
n=0 s=0 1 
4. Relation Between Steady and Unsteady Loading Funetion 
Let 
Pe, 1: Dee ee ER 
0 
1 CA RN N a ee ee} 
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From Ecg. (27) it follows that 
e > ee AFFE AT FERNEN - > en 
Eq. (75) can be written as 


F(&, U: 2 == —f Ye, U 2,) de=*2: a a ee a 
0 
Integrating Eq. (78) by parts it follows 


F(6, Yı 2) = — I Y: A a 


a A ee 


Moreover, 


RN 
76,954) = = Fey, 2)” IR | 7.,& Yı, 2) de 
0 


Or, by virtue of boundary conditions, Egq. (80) becomes 


1 oo 
= 91 a) Fa & | Pu Yı» 2) a a a 
ö 


Integrating by parts Egq. (81) leads to 
1 
Fioa,y,a)= = EI EUER SL ELEN, MN: 
DD 


Combining Eq. (77) and (82) yield: 
1 
F(&, Yı 21) = A N [Fu Yo 2) + Pa (Co Ie 2) + Pe Yo )] da 
ö 
Therefore, 
1 
F(&, y, 2) = *} [Fyu.& Y» 21) + Fa Yo 21) + FO y Zur Fear 
Hence, 2 
Fy.u(& Yı 2) + Fz2.(% Yı 21) — (0 — Deßuy.z)= 0. se: 
Using the same technique outlined above, Egq. (76) can be written as 
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F*ra,0,24)= 7 j a weh Pe KL. 7 ee er 
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From Eg. (31) it follows that 
Pin Yo 2) = Pay Yo 2) + else) een ne 


(75), 


(76). 


(77). 


(78). 


(79). 


(80). 


(81). 


(82). 
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Egs. Se and (87) lead to 
Fi, Yn 2) + Fan Yu 2ı zZ) RP, JO ....2 0.0 (88). 
Egs. (85) and (88) imply 
Far, 2) = Frl — 1 gm) Pe (89). 


Flo, yo leo = Fr VO — 1, vll 2 (90). 
Note that Eq. (90) represents the functional relation between the steady and Pr 


loading functions. 
However, from Eq. (32) the following relation holds 


Trial NEE 107 EEE on. 


Therefore, Egq. (40) and (91) lead to 


Or 


P*(z,Y. Z)l-0 = fr: (X — &1 Yu 21) 2,=0 P9.2(&ı » Yo z))2,=0 dei. - A 
N) 
Denote 
S{P6.:&r Yo z)}2,=0 = T(6, Yı z)|2=0 = al ea Pi: (&r Yo Z)la-0 dir - - - (93). 
Evidently, from Eg. (76) it follows that 
SPA Yo Z)}.=0 = Fa Yo Zu) =0 = en PC Yo Z)la=0 de - - - (9). 


Then, Egq. (92) by virtue of Eq. (93) and (94) represent the Faltung integral; namely, Eq. (92) 
can be written as 


P*z, Yo Z)l=o = PElo Yo Z)la=0 O Pal Ur Zla=0 - +» (95), 
where the sign @ denotes the functional relation by means of Faltung Integral Theorem. 
Therefore, 


E * D (£ — &» Y Zu). =0 P%,5€&1 Yı Zu). =0 der = LP: Yo 2)} 2-0 S{P6.5&r Yo )}2=0 


: (96). 
Hence, by virtue of Eq. (92), Eq. (96) leads to 


S{P*T Yo Z)}e=0 = S{P2.(% Yı Z)}.,—o SP: Er Yo 21) Yz=o a ;;: 


Moreover, Egq. (97) can be written as 


Ras i 
F*(o, Yn Z)\z,=0 _ m FE Yı Z1)|2,=0 T(a, U Zu)l,=0 IR An 35T a (98). 


1 
F* (Va — 1, y, Zu)l.-0 = 15 
&X — 


Then from Egqs. (90) and (98) it follows that 


- F2 (Va — 1, Yı, 21)).=-0 TV —1, 9, Z)e-0 (99). 


F(&, Yo Zi)l.,=0 = Trug F£ (Va? — 1, y, Bil T (Vo? —1,, 2) 


z,=0 . . . (100). 
But 


F,(& Yı A), -0 = =] u (101). 


Substituting Eq. (28) into Eq. (101) one obtains 


& 
Pub dh = 7 | Tale er rl (102). 
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Note that from Egq. (90) 


Fe nA ehem een (103). 
Therefore, Eq. (100) and (103) lead to 
1 

F(a, Yı Z)2=0 ai V@e_ı1 H,,(&; U Z)\.,=o T( x — 1; U 2) 20. 5 (104). 


Clearly, Eq. (104) can be written as 


1 
F(&, Yı» Zı)|2.=0 = 5 F,(& Yı, a) 


© = TV 1, yu a)la=o] 22.0 


Moreover, the term in the first bracket on the right side can be modified. Firstly, note 
that F,(&, Yı, 2) is of exponential order; now if there exists a function q(x,, Y,, 2) such that 


1 
{gt Up Allen = er F,@ Us Z)220 ne  C (106), 
then 


MEY: Ar J a ee (107), ° 


provided that the inverse of F',(&, Yı, 2) exists. 
The proof of Egq. (107) is trivial. 


Assume G(X,, Y,, 2,) is a function having a function F,(&, Y,, 2,) for its transform, namely 


G(&,, U A)lz.=o = SE: [F;.(&; U» Z)l.=o ee N En (108). 
Then 
2 (6) C 
Ss! Ca, 2) da, = af ern | [66 yo 2) ds de Mu (109) 
0 z=0 fi) 0 A= 
Denote 
U= [6 Dez) sods: dr ='e ade, ; 
; Bene (110) 
du = G(&,, Yı, Zu)l,=0 dt; ; DE: Gun 


Therefore, by integration by parts 


£ ( G(& Yı 21) 1, 
N) 


e 92 2 2 20 1 n 
=& - IKS Yı» Zu)la=0 ds, + = f er Gl, Yı, 2) 2 
z—0 0 0 


n A=0 
(111). 
Evidently, 
x 1 
e| [ey zdda) =, HE we (112). 
ö 
Eq. (108) and (112) lead to 
x 1 
:| Fox. yo 2) du = ne ar (113). 
0 2%, =0 
Hence, from Egq. (108), (106), and (113) it follows that 
i i ; . 114 
£& An ee Kr, Yı) an & 2.(@, Up 2) N re re ER ee ( ). 
This completes the assertion given by Eq. (106). 
However, the term in the second bracket of Eq. (105) must be developed. 
Denote 
ee MU, ee re (115), 
where 


ee ae (16). 
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But from a table of LarLAce Transforms 


srlle He — +] 0} = 7 Pin I, (a VE—R); t>k . (117), 

R— 
where I, is a modified BzsseL function of the first order. Putting a= 1 and k = u in Eq. (117), 
it follows after taking the LarLacz Transform that 


0; Fe 


<— 1 u en 
eis _ ee) gl 21V) 2.2... 118). 
(e e ) & a 2_ı2 I e “| : ( ) 
Evidently, 
eV—Ie _ ee) —- |e st _—___—_ IL lE Pd » 2... (119) 
\ (HM 
0 - 
with the condition that the integrand is zero for t< u. 
Hence, 
(e-Ver=iu _ een) (u) = en Ve Zi ..... (120). 
u 
0 


Integrating Egq. (120) from O0 to oo with respect u it follows that 


fe N}, oa [Jena at nn. ILVe— ı2)dtdu .. . (121). 


v 


Denote the left side of Eq. (121) by L and right side by R, then in accordance with Egq. (115) 
it follows that: 


T(&, Yı» 
| en u du en (122) 
& jı=0 
0 
and 
ea fu) du = TVer—1, DEN) Er er (123). 
Vo? —1 z=0 
Therefore, 
n- ee —Ina)|, ICna Kar TTE (124). 
Ve — 1 en‘ & 1=0 
However, 
R= | td fu I (125). 
Ve 
But as the integrand is zero when t < u, we have 
[6.} u: 
R=|ru| | 2) Ve ala .... 
| nn ai ) RE K u h  \ 
Evidently, ; 
1 u (u) 
R=— 2 — 2 
z$ | par, Va—w)dul. . or. 0.00 (127) 
Let ı 
ee DW. ee ee (128). 
= Mn 
Be u=VR— 7 


en te it ee 


ven 


BENENNEN EREEERRN 
rt 


m a nr at 
k 


m iu 


- M.M. Sranısıe, On the Aerodynamie Forces Acting on a Delta-Wing in Supersonic Flow ... 353 


Therefore, 


Be & :| | Le) (VER) dr) dee (129). 


Evidently, Eqs. (124) and (129) lead to 


a . 2) aan fon - | / I) RZ) dr) ra 114 0130), 
or 
ee Nana e| [ La) /Ve=7) d) adste 
Therefore, 
— a = 2110 A: Jh) Ve—®) dr) 1321 


Eq. (132) by virtue of Eq. (93) can be written as 


& VI gi EZB 

a1 T( a®—1, Y 20% =£ [76.0 Y» Z)|2,=0 + f 1,() Pl u 2 Yı» ii) 2, —=0 de, (133). 
— () 

Clearly, Eqgs. (105), (114) and (133) lead to 


F(&, Yı, 21) I.=0 


1 T & 
—uY (ar [ er’ f(r, yı) 2 x [eiac. Yo A)la=0 + | IL) Po: (Ve — &3 Y»» 2) 2 (134). 
{) (m) 


Therefore, 


ifo, Y,, Zu)}l.,-—o 


17 a en 
- (2 Se re.) 1) © (viu6 Yo z)la=o + | Ile) 76, VR—E, Yu Z)la-o “) (135). 
0 0 


Finally, 
7%, Yı 21) 


1 2 (u —5 & 
re | | e®" Kr, yo dr Irr.6. Yo z)la-o + | Ile) 7. VRR, vo Z))a=0 2 de, (136). 
0 0 0 


Therefore, Egs. (25) and (136) lead to 
Y(z, y, 2)\.=0 


1 Arte 5 er 
gr 2 | [ ePrf(r,y)dr ris6Aun and + [now Auy, Aus de, 
0 0 0 
(137). 


Eg. (137) represents the aim of this paper; in other words, the velocity potential function, 
Eq. (7) is completely defined. Hence, Eq. (12), which represents the aerodynamic forces on a 
delta-wing with supersonic leading edges can be evaluated for practical needs with no difficulty. 
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Linearized Theory of Supercavitating Flow 
through an Arbitrary Form Hydtrofoil 


By R. Osa 


Für einen Tragflügel mit einer sich bis ins Umendliche erstreckenden Kavitationsblase und be- 
liebiger Profilform wird eine linearisierte Theorie entwickelt. Diese Theorie umfaßt: I. Eine Methode, 
um für eine vorgegebene Druckverteilung die zugehörige Form des Flügelprofils zu erhalten; 2. eine Me- 
ihode, um die Auftriebscharakteristiken eines vorgegebenen Flügelprofils zu berechnen. Einige besonders 
wesentliche Profilformen werden numerisch durchgerechnet. 


For an arbitrary form supercavitating hydrofoil, the author developed a linearized theory which 
included the following methods: (i) a method to obtain the hydrofoil form for a given pressure distribution, 
(ii) a method to estimate the characteristics of a hydrofoil form. Some excellent kinds of hydrofoil-form 
were calculated numerically. 


PasBuBaetcn JIMHeÄaHaA Teopun ANA Kpblla C ÖecKOHeyYHOH KABHTANUHOHHOH IIOJNOCTLP. 
IJTa Teopun CcoNTep>KuT 1. METON HOJAIYYCHHA COOTBETCTByWIMeli POPMEI IpobmuAa Kpblaa mu 
MAHHOTO pPacııpejesieHuHA MaBJIeEHMA, 2. MEeToN pacyeTa XapakTepHCcTHK TONYUSMHOÜ CHJIEI. F 
Ilpneneubl UNCJIeHHbIe PacyöTbI IA HEKOTOPBIX OCOÖEHHO CYINECTBEHHBIX BOpM npobmum. | 


1. Introduetion 


With speeding up and the increased power of hydraulic machineries, there arise three well 
known kinds of trouble due to cavitation: (i) Heavy drop in the hydraulic efficiency of the ma- 
chine; (ii) damage and destruction to the machine from erosion; (iii) vibration and noise evoked 
by cavitation. It can be suggested to use supercavitating hydrofoils in these hydraulic machineries 
to keep them out of these troubles. For, if we adopt these hydrofoils, we may not only get rid 
of the above mentioned troubles (ii), (iii) no doubt [1]}), but also the trouble (i) fairly well, be- 
cause of the superior high-speed-characteristics [3] for a supercavitating model ship-propeller to 
use Tuin-hydrofoil [2]. 

Well, the supercavitating flow through a hydrofoil can be treated as two dimensional dis- 
continuous flow, theoretically. According to literature, it is true that Levi-Civita [4] proposed 
an analytical method for two dimensional free-streamline flow about an arbitrary form body 
but, it is essentially difficult to obtain the solution for a given arbitrary hydrofoil in the physical 
plane or a given arbitrary pressure-distribution. And also TuLın [2] developed a linearized 
theory for supercavitating flow through an arbitrary hydrofoil, but Tuzıw’s solutions are doubt- 
ful, for they include some mistakes in the given boundary conditions and no consideration is 
given to the necessary conditions for the solution. 


!) Numbers in parentheses refer to the bibliography at the end of the paper. 


fe >> + 5 W Wer . ‚ 

i u u A N 2 u N, - 

n erer R el . er "oa 
ar u ’ > ai 


Don . 
Ge 2 
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Accordingly, in this paper, for the purpose of excellent hydrofoil research, the author 


developed a theory which depended upon a linearized method like TuLın’s theory [2] and included 


next matters: 


- () A method to obtain the supercavitating hydrofoil for a given arbi 2 
tribution (Inverse method). s =. EINENSArSILTAEYSRTESSurEzUiG 


(ü) A method to estimate the characteristics for an arbitrary form itati i 
hen) | y supercavitating hydrofoil 
_@üi) The necessary conditions for the supercavitating flow and the upper limit of the lift- 
coefficient and so on, are clarified by numerical calculations. 
(iv) And besides, a few kinds of excellent hydrofoil-form are calculated, of which the theore- 
tical drag-lift ratios are smaller and the thickness-distributions of the trailing cavities, 
that is, the thickness-distributions of the hydrofoils, are better than that of TuLıw’s hydrofoil. 


2. Basie Equations 


To simplify the problems, it is assumed that the chord-length of an arbitrary supercavi- 
tating hydrofoil y(x) in the Z-plane is a unity, the leading edge is at the origin of coordinate, 
the trailing cavity length is infinite and the magnitudes of the disturbance velocities V=u-—iv 
due to the presence of the hydrofoil in the main stream are small in comparison to the magnitude 
of the main stream velocity U,. Then, the next relations hold for the bottom surface of the 
hydrofoil, see Fig. 1: 


dy v 

dx == D, ee Me, u ar Haren Me Tee ee te ar aa Er ee ae (1), 
B>p, 2u w 

(„= 1 KlwaTr ar (2), 


where, the pressure coefficient expressed by Eq. (2) must be always positive, from the necessary 
conditions for no cavitating on the bottom surface. Next, the whole Z-plane where the super- 
cavitating hydrofoil is placed, can be comformed into the under half ö-plane by the well known 
Turın’s transformation [2], > 

VRR NE ee (3), 


where, Z=r +iy,i=$ 5 in. In the Z-plane, the flow is changed into the one through a 
'non-cavitating aerofoil and also is made completely correspond to the one in the Z-plane, see 
Figs. 1 and 2. 


Fig.1. Linearized boundary conditions in the physical Z-plane Fig. 2. Linearized boundary conditions in the mapping Z-plane 


Then, the next relations hold for the bottom surface: 


A TIRPeR N (4,), 
en ER REN; 7a). 
Den RE ER (4,), 
Ve ee TR BEER 


where, Eq. (4,) and Eq. (4,) are equivalent each other and indicate the condition for the com- 
formal mapping. Egq. (4,) expresses the conservative law of the circulation for referring pair 
points in each mapping plane. And dZ/d£ can be given by 2£ for the bottom’ surface in this 
linearized theory. 

What is more important to make an additional remark: In Tuuin’s theory, it is assumed, 
that the point —oo in the £-plane refers to the point —0o in the Z-plane and the velocity at 


the point —oo is equal to U.. Obviously these assumptions are mistaken. 
24* 


T 
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From the above equations, the next relations among the lift-coeffiecient C,, the drag- 
coefficient C, and the moment-coefficient about the leading edge Cy, for the hydrofoil an 
the aerofoil, that is, so called hydrofoil-aerofoil equivalence relationships hold: 


EVER RT a rin 6), 
Cyu(&) = Cys(&) a 5 SR Fr ee eh Er (6), 
Cola) = CHOlEn - . M. 


3, The Direct and the Inverse Method for a Supereavitating Hydrofoil 


In consideration of the conservative law of the eirculation, Eq. (4,), and the sudden change 
in the pressure distribution near the leading edge, the bottom-surface-pressure-distributions in 
the Z-plane and the Z-plane can be expressed as follows: 


E= = 5-10) en ee A Fe (8,), 
u)  ulE 0 eg 
= a ir hu (8;), 
where 029 >—n, 2 
Then, 
A, 
Ge 2a (4 2% or. Ei... 9), 
7 As 
ud) =Z(A+Aı-7) 2 ee (9,), 
7 A, 
Cu) = 35 (540 +74 —74, +3 a (9,). 


Next, the bottom hydrofoil-profile-slopes, Eq. (4,), are 


dy dn = 
Fe Fe Zr N (10). 
In the inverse method, the profile form can be calculated from the pressure-distribution 


which is the given pressure-distribution minus A, ah ‚ because A, isregardless of the profile 


2 
form and is easily estimated from the pressure-distribution near the extreme leading edge. Then, 
ua) u) _ X £ 
Fe N 2 A. sinn 0. 0. SE Pe (8, ), 
ay _d_NS 
har 2 An COS. ee (10). 


On the other hand, the bottom profile form is expressed as follows: 


yo + un MR TER 0 RE (11) 
n=1 : 
where the next relation must be satisfied, in order that the leadi igi 
| » ead a 
ee; ing edge need to be at the origin 
1 A 
0 2 Up n in mtl allen Sa (IR): 


Here, new auxiliary FOURIER-constants, C,„ defined as follows, are introduced into. 


oo oo 
23na,sinnd=—sindYC,cosn®.... (12) 
2 Bi rn. 5 


For C,, An, the next simple relations hold, 


1 
4=— 26-0), | 
1 
an Fi irn (Con-ı — Can+)> N u ee RE (13), 
Ar: | 
an+l (Con — Con+2) ) 


Du (2n+ 1) 


Ä 
N 
ö 
\ 
t 
h 
j 


R. Opa, Linearized Theory of Supercavitating Flow through an Arbitrery Form Hydrofoil 357 


or 


2 2 ®n +1 43241 


Co Bonn hr ul Fe e N (133% 


n 


TE ) 


N 


Then, from Eqs. (10), (11), 


Bf! 1 
Do 2 | Gag Ara] 


2) 
: ee, 
mA (An+ı t An-ı); 1,2 


After all, the profile form which refers to the pressure-distribution given in Egqs. (2), (8), 
(11’),-is determined in Eg. (11) from Egs. (13), (14). On the other hand, the characteristics, 
that is, the pressure distribution, the lift-coefficient and so on, are easily obtained?), for a pro- 
file form given in Egs. (8,), (10°), (11), (11’). 

There are, however, the limits of application in the incidence-angles for these hydrofoils, 
from the necessary conditions for the supercavitating flow. These limits will be discussed in 
the next clause. 


4. The Free-Streamline Form in Upper Profile Side and the Necessary Condition for the 
Supercavitating Flow 


We will calculate the free-streamline form which refers directly to the hydrofoil strength. 
At first, the free streamline form should satisfy the next relation, 


Eu au Be Rs); 


ee ee 0 eh Kan 0 er lmer Lie: 


Ye) = 7; = v(a) dx = | I gas Fr NR (16) 
0 oo 
And v,(£) in Egq. (16) is 
| in = i ie nn = ee SE (17). 
E “; 
r Put, ee a Ni, 


s Besides, putting Egs. (8,), (18) into Eq. (17) 


ve) 1 [ sind 
; om) t— cosd 


[4 il 5 + &Ausinn0lao= As IH ZA, I ale 


where 


Let 1}, 
L,=—-1+21t], 
Bear, her, Rn z2; "Integer. 


For a special case, J, #0, I„ = 0, of a flat-plate hydrofoil, 


ve) _ = er RE N (20). 
Dem Be 


2) When we caleulate An (n > 1) from Eggs. (10’), (14), we determine An as linear equations of A, at 
first and then A, from Eg. (10). 
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pproximate solution [2] obtained by the comformal mapping. 
Then, the free streamline form for a flat-plate hydrofoil also agrees with the exact one [6], in 
the accuracy leaving the errors under 3” order small. The relation between x and yy4,-o 
obtained by Egs. (16), (20), is tabulated in Table 1. The free streamline form has been obtained 
for A, J, in Eg. (19) already. Then, we will calculate the form as A, = 0 and superpose it the 
solution for a flat-plate hydrofoil, in the following calculations. 


This agrees with the Turin’s a 


Table 1. The bottom profile froms y and the free streamline forms y; for various parameter An 


0 0 0 0 10 0 0 0 0 

0.0045 —0.0045 0.0191 0.0041 0.0030 —0.0030 | —0.0020 0.0014 0.0014 
0.0137 —0.0137 0.0415 0.0115 0.0080 —0.0059, | —0.0047 | —0.0013 0.0029 
0.0319 —0.0319 0.0740 |- 0.0243 0.0166 —0.0056 | —0.0088 | —0.0144 0.0038 
0.0625 —0.0625 0.1162 0.0417 0.0291 0.0052 | —0.0139 | —0.0438 0.0069 
0.1082 —0.1082 0.1669 0.0608 0.0454 0.0348 | —0.0197 | —0.0872 0.0090 
0.1707 —0.1707 0.2245 0.0767 0.0651 0.0889 | —0.0259 | —0.1304 | 0.0109 
0.2500 —0.2500 0.2868 0.0833 0.0874 0.1667 | —0.0322 | —0.1500 0.0127 
0.3444 —0.3444 0.3514 0.0749 0.1115 0.2591 | —0.0384 | —0.1253 | 0.0143 
0.4503 —0.4503 0.4160 0.0474 0.1362 0.3502 | —0.0443 | —0.0510 | 0.0157 
0.5625 —0.5625 0.4779 0.0000 0.1605 0.4219 | —0.0494 0.0563 | 0.0168 
0.6747 —0.6747 0.5348 —0.0642 0.1832 0.4602 | —0.0540 0.1634 ' 0.0177 
0.7797 —0.7797 0.5846 —0.1383 0.2033 0.4605 | —0.0578 0.2372 ı 0.0184 
0.8705 —0.8705 0.6254 —0.2124 0.2200 0.4301 | —0.0609 0.2662 | 0.0190 
0.9406 —0.9406 0.6557 —0.2757 0.2324 0.3857 | —0.0632 0.2442 0.0194 
0.9849 —0.9849 0.6744 —0.3183 0.2402 0.3480 | —0.0645 0.2141 0.0197 
1.0000 — 1.0000 0.6807 —0.3333 0.2427 0.3333 | —0.0649 0.2000 | 0.0197 


As the form defined from the leading edge till the trailing one is quite enough for our use, 
a next parameter p can be introduced into. 


1 
&=— 7 (1—cosg), ESP Sm.nin RE 


Then, from Egs. (16), (19), (but A, = 0) 


4 ? 
11, ine r ee oo 

Yan J 5 A,1,sin og (1 — cosp) dp = | Y Al 3 Imnsinmg) dp . . (22). 
n=1 ö n=1 m=1 


0 


On the other hand, the free streamline form can generally be expressed as follows: 


y4,-0= Z dncosmp + Z eu sinmo Wer De ee 
m=0 m= 
where ar 
EN NR 
m=0 
From Egs. (22), (23) 
— mdun=E Imfss U. sum TE 
n=1 


EN De of the profile form, fn„ are previously able to be calculated, as shown in 
able 23). 


?) /mn are calculated from the FOURTER-analvsis 36 poi i 
D ysis of 36 points summation-method. Whi 
A, can exactly be caleulated from Eg. (16) putting directly Eq. (19) into. So that es 


una? _e_1yye lan er ES 
N 3 | a 3 zeig ler VE 3 
4,[# 2 14 Din 1 
Ya =——2|- —t Be re An Ne, 4 l \ 7 —— t 5 1 | 
2 2 [2 +75 B) 3 Ye N+z+rzY@e—-D+ ZA ZT —lgle+ Yaı) 


The accuracy of these fmn is sufficient enou 
2 ss enough for our use, because the exact i 
the FOURIER-analysis agree with each other, up to the fifth significant ern N. 


Baht Zu al a a an ud dd nn ns nn 


” 


u a 7% w h Lu: 
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Table 2. Parametrie coefficients —fmn/m for the free streamline forms 


n 
re 1 2 | 3 | 4 5 6 7 8 | 9 
m: 
1 1—0.12500 1—0.03375 —0.01003 —0.00332 —0.00124 |—-0.00052 |—0.00025 |—-0.00013 __0.00007 
2 0.01688 | 0 —0.00151 —0.00102 —0.00056 |—-0.00030 —0.00016 |—-0.00009 |—_-0.00005 
g 0.00334 | 0.00100 | 0 —0.00022 —0.00020 |—0.00014 |—-0.00009 |—-0.00006 __0.00004 
4 0.00083 | 0.00051 | 0.00017 | 0 —0.00005 —0.00006 —0.00005 |—0.00003 __0.00002 
5 , 0.000253 | 0.00022| 0.00012 | 0.00004 | 0 —0.00002 —0.00002 —0.00002 |—-0.00001 
6 0.00009 | 0.00011 | 0.00007 | 0.00004 | 0.00001 | 0 —0.00001 —0.00001 —0.00001 
7 0.00003 | 0.00005 | 0.00004 | 0.00003 | 0.00001 | 0.00001 | 0 0 0 
8 0.00001 | 0.00002 | 0.00002 | 0.00002 | 0.00001 | 0.00001 | 0 0 0 
9 0.00001 | 0.00001 | 0.00001 | 0.00001 | 0.00001 | 0.000010 0 ) 


Besides, the profile forms and the free 
streamline forms are calculated for each A, 
The results are shown in Fig. 3 and Table 1. 
Where, A, is given in a negative value as 
the drag-lift ratio is very large for a posi- 
tive A,. 

From these figures, the thickness?) of 
the trailing cavity defined as the vertical 
difference between the free streamline and 
the bottom profile surface, is positive over 
the chord-length only when A,=#0, A,„=(0, 
n>0. Then, the supercavitating flow is not 
necessarily realized, responding to an arbi- 
trary pressure-distribution or an arbitrary 
hydrofoil-form. 

If we give the minimum pressure in 
the pressure-distribution by the trailing-edge 
and A,=(, to obtain the low drag hydro- 
foil as the TuLin’s suggestion, A, and so the 
thickness of the trailing cavity becomes 
negatively large. Then the drag-lift ratio 
obtained in such a way is scarcely significant. 

As we decide whether the supercavi- 
tating flow appears or not, at first, we cal- 
culate the thickness-distribution of the trail- 
ing cavity for a given pressure-distribution 
or a given bottom profile form by the above 
mentioned method. If the thickness is posi- 
tive over the chord, the supercavitating flow 
appears. On the contrary, if negative at any 
part of the chord, no supercavitating flow 
appears. To make the flow to appear for 
later case, then, the incidence angle need to 
be increased by & from the one of the pre- 
liminary design state, till the thickness is just 
positive over the chord. The lower limiting 
case in which the thickness is positive over 


the chord, is the one of the supercavitating Fig. 3. The bottom PR y an the free streamline forms Y 
fl] for various parameters A, 
OW. 


(d) Case of A, = 1.0 


5. The Maximum Lift-Coeffieient 


There is a limit of lift to use in such a supercavitating hydrofoil. The maximum lift- 
coefficient will be discussed in this clause. 

The maximum lift-coefficient5) Cimax can be defined as follows, from the condition in 
which the bottom surface-flow-velocity can not be lowered below zero, and from the consideration 
for the flow separation on the bottom surface and a somewhat bold assumption in which Eq. (8,') 


4) As well known, the upper profile-form isarbitrary on condition that the upper profile-surface is 
within the trailing cavity. 

5) This OL max is the maximum lift-coefficient for zero cavitation number. When the cavitation number 
is not zero, the maximum lift coeffieient can approximately be calculated from the Brrz’s correction [7]. 
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holds good for such a large lift-coefficient, 


> Ao__ As ) 
2 ( u A :34, 
To rer aD Fan a 
sind + —cot — + —sinn 
| A, 2 n=2 A, min 
where, Eq. (25) must not be applied for the velocity-distribution of the extreme neighbourhood 
of the leading edge and is only connected with the velocity-distribution u/C;. 


In an extreme case, C;max for the Turın’s minimum-drag-hydrofoil of which lan 
distribution is expressed by an impulse function, u/C, at the trailing edge, is zero. Such a profile 


is not worth considering. 


ne hand, it is desirable to use the considerably lower lift coefficient than Graa 
reine easy to separate on the bottom-profile surface and then the profile-drag is pe to 
increase for such a large C, a8 Cymax- Of course, it’s degree, is different due to the velocity- 
distribution. On the other hand, if a somewhat larger incidence angle than the designed one 
is given, as shown in Ecg. (5), (9), the larger value than the designed C;max can ‚be selected 
for the profile. The two factors are canceled each other. Then, at least, C;max defined in this 
clause has a right for a measure of maximum lift-coefficient to use. 


6. Numerical Examples 


(1) Turın’s Profile and It’s Family . 

Here, we take up a profile-family defined only by A, A,, A,, which are directly connected 
with the drag-lift ratio, see Egqs. (5), (7), (9,), (95). For two typical examples showing the nature 
of this family, A,/A, = — 1/2 (Turin’s profile [2]) and A,/A, = — 2 (having about minimum 
drag-lift ratio in this family), the profile form, the free streamline form and the veloeity-distri- 
bution for each minimum incidence angle is calculated and is shown in Fig. 4-5. The theore- 
tical drag lift ratios for these profiles are also shown in Fig. 6. 


Eulen n..2 
SARBENEEN 


TILLLBEER 
SET 


Fig. 4. The bottom profile form y, the free streamline form Y Fig. 5. The bottom profile form y, the free streamline form Y 
and the velocity-distribution for TULIN-profile and the velocity-distribution for the case of A,/A, = — 2 


The characteristics of Tuzın’s profile is exceptionally more excellent indeed than the 
ordinal screw-propeller-hydrofoil-profile, for example, the ogival profile of which theoretical 


drag-lift ratio epıate = — Cr, but also, there is a fatal weak point in the form, in which the 
IT 


thickness is very thin®) by the leading edge. And then Turıy’s profile is unfitted in practice, 
Being quite a similar to Tuzın’s profile, a eircular-arc profile which is studied by Wu [8], also 
has the same fatal weak point. 

In the case of A,/A, = —2, the thickness is fairly solid near the leading edge and the 
maximum-thickness-point on the chord goes backwards fairly. So, the additional incidence- 
angle for thickening the profile and the drop in the lift-drag ratio are small. Besides, the drag- 
lift ratio in design state is smaller than the one of Turım’s profile by 7%. Therefore, this profile 
(A,/A, = — 2) should be more excellent than Turıw’s one. 


| %) The incidence-angle need to fairly beincreased from the designed one and then the theoretical drag- 
lift ratio is apt to increase, j 


Kate... 


w 
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In Fig. 6, Tulin (IMG) indicates the theoretical drag-lift ratio for so called Tuuıy’s profile 
designed by Turn, in which no supercavitating flow appears, and it differs considerably from 
the one in supercavitating state. Then, it is mostly meaningless to discuss hydrofoil charac- 
teristics without the necessary condition for supercavitating flow. 


Fig. 6. The relationship between the lift-coefficients C '„ and the theoretical drag-lift ratios e for various profile forms 


Still, there is a most suitable A,/A,-value for the theoretical drag-lift ratio. Because the | 


i drag-lift ratio for A,/A, = — 00 is 0.141 C, and it is larger than 0.138 C, for A,/A, = —2. 
Y (2) Hydrofoil SC1 
4 For an example to the inverse-method, the velocity-distribution in A, = 0 is given as 
. in Fig. 7 and Table 3. 
ns seele Fee ne 
! Cı 10 — N 
EEE 
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> ee 
Cı 0 Ez 
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02 
ee TG LLLECH 
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Fig. 7. The velocity-distribution —2 u/C, [0 } Es an a 
given in hydrofoil SC1 Y 
BE te 
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Fig. 8. The bottom profile form y, the free streamline form Y 
and the velocity-distribution for hydrofoil SC1 


The minimum ineidence-angle «min is 0.113 C, rad., from the necessary condition for super- 
cavitating flow. The bottom surface profile-form, the free streamline form and the velocity- 
distribution for &min are shown in Fig. 8 and Table 3. For &nin, the drag-lift ratio is 0.130 C, 
and the maximum lift-coefficient C;max is 0.66. Then, this profile may be an excellent profile 
too, because the theoretical drag-lift ratio is smaller than that of Tuuin-profile by 13% shown 
in Fig. 6 and the thickness near the leading edge is not very thin, see Fig. 8. 


Well, TuLin touches upon that the profile-characteristics are remarkably affected by the 
profile-form near the leading edge. Then, it is interesting whether the solution can be obtained 
or not, as the solution expressed by good convergent series for the velocity distribution ex- 
tremely changed near the leading edge. Then we will check the above mentioned problem for 
the velocity extremely changed near the leading edge in this example. 
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Table 3. The bottom profile form y, the freestreamline form yand the 
veloeity-distribution for hydrofoil SC1 


mm ee 
x | —2u/Cz- given | y/Cı | wien | —2wer . 

RER N 
0 N) 0 0 oo 
0.00006 0.069 0 0.0001 2.629 
0.00091 0.112 0.0000 0.0008 1.369 
0.0045 0.114 0.0001 0.0028 0.934 
0.0137 0.077 0.0002 0.0063 0.682 
0.0319 0.027 0.0004 0.0119 0.504 
0.0625 0.000 0.0017 0.0196 0.390 
0.1082 0.040 0.0059 0.0293 0.354 r 
0.1707 0.178 0.0158 0.0409 0.415 
0.2500 0.424 0.0330 0.0538 0.575 
0.3444 0.758 0.0560 0.0676 0.814 
0.4503 1.131 0.0793 0.0817 1.089 
0.5625 1.470 0.0939 | 0.0955 1.341 
0.6747 1.698 0.0909 | 0.1084 | 1.503 
0.7797 1.748 0.0661 0.1199 1.522 
0.8705 1.583 0.0236 0.1293 1.365 
0.9406 1.203 —0.0249 0.1364 1.031 
0.9849 0.650 —0.0632 0.1407 0.555 
1.0000 0 —0.0777 0.1422 0 


Table 4. Convergeney of the series, @„, expressing the hydrofoil SC1 


Go | ga Gy, Qs | 2 Q; Ag Qr 


| 
0.07419 | —0.07947 | 0.026538 | 0.058300 | —0.02649 | 0.005331 | —0.00001 | 0.00000 


The series, a,, for the profile form tabulated in Table 4. In Table 4, it may be appreciated 
that the higher terms from a, are converged enough for our use. Then, the convergency of the 
FOURIER-series is sufficient enough, because the series converges enough even in such an extreme 
velocity-distribution in this example. 


3) Hydrofoil SC2 


For a numerical example for direct method, the bottom surface form is given, in the first 
half in a straight line, in the second half in the similar form to hydrofoil SC1 whose excellent 
characteristics have been acknowledged in the former example, as shown in Fig. 9 and Table 5. 


_2u 15 


Fig. 9. The profile form viC, given in hydrofoil SC2 


Fig. 10. The bottom profile form y, the free streamline form Y 
and the veloeity-distribution for hydrofoil SC2 


il From the necessary condition for supercavitating flow, the minimum incidence- N 
is 0.0048 97 rad. For Amin; the drag-lift ratio is 0.123 C, and Cz;max is 0.66. And PR 
surface profile form, the free stream-line form and the veloeity-distribution are shown in Fig. 10 
and Table 5 for Xmin. This profile may also be excellent, because & is the smallest for the present 
the thickness is scarcely thin and the profile form is simple, i 
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Table 5. The bottom profile form 4, the free streamline form Yı 
and the velocity-distribution for hydrofoil SC2 


rer ee ee 
= | wOr-given | cr | wien | —2ujC: 


0 0 0 0 oo 
0.0001 0.0000 0.0000 0.0001 1.666 
0.0009 0.0002 0.0002 0.0006 1.020 
0.0045 0.0008 0.0008 0.0021 0.796 
0.0137 0.0024 0.0023 0.0051 0.621 
0.0319 0.0056 0.0053 0.0100 0.478 
0.0625 0.0110 0.0104 0.0170 0.409 
0.1082 0.0191 0.0180 0.0260 0.431 
0.1707 0.0301 0.0284 0.0367 0.522 
0.2500 0.0440 0.0415 0.0488 0.668 
0.3444 0.0607 0.0572 0.0617 0.880 
0.4503 0.0793 0.0748 0.0751 1.147 
0.5625 0.0939 0.0883 0.0883 1.392 
0.6747 0.0909 0.0850 0.1005 1.506 
0.7797 0.0661 0.0605 0.1114 1.432 
0.8705 0.0236 0.0188 0.1204 1.193 
0.9406 —0.0249 —0.0286 0.1271 0.846 
0.9849 —0.0632 —0.0660 0.1313 | 0.439 
1.0000 —0.0777 —0.0802 0.1327 | 0 


%. Conelusion 


The results obtained in this paper are summarized as follows: 


(1) The author proposed the direct and the inverse methods for an arbitrary supercavitating 
hydrofoil, which depend upon a linearized method, for the purpose of excellent hydrofoil research. 


(2) The Tuzis’s suggestion [2], in which the low drag hydrofoil is considered to be obtained 
if only A, = 0 is given and the maximum pressure is given near leading edge in the pressure- 
distribution, is mistaken. In the supercavitating hydrofoil design, we must think of the necessary 
condition for the supercavitating flow, the maximum lift-coefficient Cimax and the fitness in 
the thickness-distribution of the trailing cavity. 

(3) By the numerical caleulations, it is shown that the ordinal profile, as example, TuLıx- 
profile, has a fatal weak point in the strength, because of the extremely thin thickness by the 
leading edge. And also a number of excellent profiles are caleulated out, in the manner in which 
the profiles must have the lower drag-lift ratio than Tvrı-profile, the suitable trailing-cavity- 
thickness-distribution and so have no weak point in strength. 


In conclusion, the author expresses his cordial gratitude to Prof. F. Numacht, the Director 
of the Institute of High Speed Mechanics, Tohoku University, for his generous guidance through- 
out this study. 
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Schranken für Schwingungsdauer 
und Lösung bei der freien ungedämpften Schwingung 


Von H. EHRMANN 


i * ” * ” Fr Fr nd 

Mit Hilfe eines einfachen Vergleichsverfahrens werden Formeln zur Abschätzung der Lösung u 
im Falle a Lösungen der Abhängigkeit der Schwingungsdauer vom Ausschlag für die freie unge- 
dämpfte Schwingungsdifferentialgleichung hergeleitet. Sodann wird zur genaueren Berechnung dieser Größen 
ein Näherungsverfahren angegeben. n 

By means of an easy method of comparison estimates for the solution and, in the case of periodic solutions, 
for Ihe period Eee on the amplitude are derived for the differential equations of free undamped oscilla- 
tions. Further an approximation method is stated for more exact caleulation of these quantities. 


C NOMOINBIO IPOCTOrO MeTona CPABHeHUA BhIBONATCH POPMyMIbI II ONCHKH penieHun 
nudhepenmmarnbHoro YPaBHeHHA CBOOONHOTO He3aTyxamınero konebanna. B cıyyae mepmo- 
NHYeCKUX PeIIeHMÄ NAIOTCH ONMEHKN 3ABHCHMOCTH TIepmoNa KoleÖaunuA OT AMILIHTYIBL. aa 
YTOUHeHUA ITUX BeImyun Naerch IPMÖNMKEHHBIÄ METoN. 


1. Die Umkehrfunktion der Lösung x(f), x(tj) = X (ty) = X, der Schwingungsdifferential- 
Kerhune Etfd=0 „0.2 ee 


läßt sich bekanntlich geschlossen darstellen durch 


f de 
t—1,= | —— — 
j Yi3— 2 [F(@x) — F(@o)] 


und die Abhängigkeit der Schwingungsdauer T vom Ausschlag a im Falle periodischer Lösungen 
ist unter der Voraussetzung, daß f(x) eine ungerade Funktion ist, gegeben durch 


mit Fa) = f (dd ....:.() 


a 


na - 2/2 | = el... | 


YF(a) — F(«) 


Ist die „Kennlinie‘ f(x) eine ganze rationale Funktion zweiten oder dritten Grades, so stehen 
in (1.2) und (1.3) elliptische Integrale. Im allgemeinen Fall ist jedoch die Auswertung der un- 
eigentlichen Integrale bei der Untersuchung periodischer Lösungen mit einem gewissen Arbeits- 
aufwand verbunden!). Um schnell’einen Überblick über den Verlauf der Lösung und die Größe 
der Schwingungsdauer zu erhalten, ist daher der im folgenden durchgeführte einfache Vergleich 
mit linearen Gleichungen, der Schranken für diese Größen liefert, häufig nützlich. 

2. Grundlegend dafür ist der einfache Satz: 

Gegeben seien die beiden Differentialgleichungen 


RE RE 
Et. ee 


Die Anfangswertaufgaben x;(0) = 0, x;(0) > 0 seien eindeutig lösbar, insbesondere existieren 
die Integrale 


0 


Fi) = [ fla)de (i=1,2) für O<r Sa. 
0] 


Ferner gelten die Gleichungen 


2 Fa) = Fila) A D>U re 
un 

F&) so) <A er Daran 
Dann gilt für die Lösungen x;(l) mit den Anfangsbedingungen 

x(0)=.0,.2(0)=Y24. ((=1,29)..00000 0 re 


Es existieren zwei Konstanten r, und t, derart, daß 
.c)=a, )=0, SH>O Für 0 <t<y (Ü=12) 
gilt, und die Umkehrfunktionen der Lösungen genügen den Ungleichungen 


ba) Sta) für O<saısa, n=hla) <tl)—=n,. 


') Mit Hilfe einer von K. KLoTTER angegebenen Transformation lassen sich diese Inte 


[ L grale näherungs- 
weise nach der Sımpsoxschen Regel berechnen. K. KLOTTER, T g 


echnische Schwingungslehre, Berlin 1951, S. 152ff. 
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Beweis: Es gilt wegen (2.5) &() = 2 (A — F;(x)) und 


1 f dx 
er B 
y2 | YA— F;(@) u 


woraus sich wegen der Eindeutigkeit der Anfangswertaufgabe, (2.3) und (2.4), unmittelbar die 
Behauptung ergibt. 


Ist f;(a) > 0, so sind die Lösungen x;(l) symmetrisch zur Geraden t = r;. Gilt noch I) 
= — f(x), so ist die Schwingungsdauer der betreffenden Lösung T; = 4 r,. 


3. Der Satz bietet eine sehr einfache Möglichkeit, Schranken für die Lösungen und die 
Schwingungsdauer bei der Gleichung (1.1) anzugeben. Betrachten wir zunächst als Vergleichs- 
differentialgleichung die lineare Gleichung 


an a EN GM) 


wobei die Funktion ®?(a) gegeben ist durch 


oa) = Fa) En < (3.2). 


Diese Gleichung besagt wegen 


a 


a 
2 a 
F(a) = | j@) dx = j oXa)adı = a) REM? (3.3), 
0 0 

daß die Flächen unter den beiden Kennlinien f(x) und »?(a) x gleich sind. 

Wir setzen jetzt F(a) > 0 und F(x) < F(a) für 0 <x <a voraus. Gilt dann im ganzen 
IntervallO <x< a stets 

x? w*(a) 
2 

so sind die Voraussetzungen des obigen Satzes erfüllt, und wir haben für die Zeit 7 der Lösung 
xz(f) von (1.1) mit den Anfangsbedingungen x(0) = 0 ,&(0) = Y2 F(a) > 0 vom Nullpunkt bis zum 
ersten Maximum die Abschätzung 


22 w%(a) 


<F(&) bzw. SCHE (3.4), 


NEAR 1 “ 
| EA Eee VR(a) bewe en Te ea (3.5). 
Die Bedingung (3.4) ist erfüllt, wenn die Funktion &%x) für 0 <x <a monoton nicht 
steigt bzw. monoton nicht fällt. 
Besitzt f(x) die erforderlichen Stetigkeits- bzw. Differenzierbarkeitseigenschaften, so kann 
man weiter schließen: 
Die Bedingung (3.4) also 


Fra are seh. zU. md, al rn... : (3.6) 
ist erfüllt, wenn in diesem Intervall stets 
<0 bzw. >00 
gilt, woraus durch Differentiation von (3.3) die Bedingung 
Bra DE bzw. ed: Da ee (3.7) 
folgt und umgekehrt. (3.7) ist erfüllt, wenn 
ande else ZU nn Wal !.. ... m... (3.8) 
gilt. Schließlich folgt (3.8) aus den spezielleren Bedingungen 
f0) =0 und f’(«) <O bzw. .f(0) <O und @)=0 in [I,al ... 9) 


oder auch aus den Bedingungen 
lim f(x) =0 und f(x) monoton nicht zunehmend bzw. 
REEL | on 
lim f(x) <0 und f(x) monoton nicht abnehmend in [0, a] 
z> 70 
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Unter diesen Bedingungen gilt ferner für den Anfangsbogen 0 =! =7, D Er EA 


(2) > i(x), also x(t) = &(t) sin ED —_. : 


wobei &(f) die Lösung der linearen Vergleichsdifferentialgleichung (3.1) und i(x) ihre Umkehr- 


funktion ist. } , 
4. Wir betrachten weiter die Anfangswertaufgabe (1.1), x(0) = 0, &(0) = 2, > 0, lassen 


jedoch die Bedingung &(a) = 0 fallen und fordern nur 
2EW)<en für sroa, EEE (4.1). 
Gelten dann für zwei Funktionen g,(a, x) und 9,(a, x) die Gleichungen 


Flo) = [gla,o)de, 2/@)de<i} (=1,2) für Osz<a... (a) 

und : y h 2 
F(x) <a 2)de und Fe) 2 /g(a,2)ds für Ozsrsa.... (43) 

0 0 


so verläuft der Anfangsbogen der Lösung x(f) der Differentialgleichung (1.1) mit den obigen An- 
fangsbedingungen für 0 <x <a zwischen den Lösungen der Gleichungen 
ut), N) =lt, Ar er rrze (4.4). 
Dies folgt genau wie in 2. leicht durch Einsetzen in Gleichung (1.2). 
Als Vergleichsdifferentialgleichungen kommen in der Praxis nur solche in Frage, deren 
Lösungen leicht angebbar sind: 


Ha, 2) — ala) ala) FE Er u I 
mit der Lösung mit den obigen Angangsbedingungen 
Pa 2 sinol+ 2 0sol—.% mit = c(a) 
[0 


bzw. im Falle c,(a) < 0 mit x = Y— c, 


rl re base I RL 
:)=5[, me ee Te 


Die Konstanten c, und c, sind dabei so zu wählen, daß die Bedingungen (4.2) und (4.3) für 
g = 9, bzw. g = 9, erfüllt sind. Der Spezialfall c, = 0 wurde in 3. behandelt.?) 

5. Analoge Betrachtungen wie in 3. führen bei dem Spezialfall c,(a) = 0 in Gleichung (4.5) 
auf folgende Ergebnisse: 


F(a) 
Die Vergleichsdifferentialgleichung lautet also jetzt 
wi 
ao... 2 
a 
3 ande: Pay Due en MER ; a 
Ihre Lösung ist x(l) = — 3 P+%,t. Bezeichnen wir die Umkehrfunktion mit ?(z), 
so ergibt sich im Falle &, = &, > 0 mit (4.]): 
Gilt im IntervallO <x = a stets 
a Fe) =2.Fla)'S0 (ba a0 I 
so ist dort x) < ix) bzw. x) = x). (5.2) ist erfüllt, wenn in [0, a] stets 
za) — Fo) 20 bzw. ET re 
oder spezieller, wenn 
fe) 20: bzw. SO RS En se 


oder wenn /(x) mit wachsendem x nicht fällt bzw. nicht steigt. 
6. Der Fall c, cs + O in (4.5) läßt verschiedene Möglichkeiten zu. Wir betrachten hier nur 

den Sonderfall, daß mit (4.2) auch zugleich g(a, a) = f(a) gilt, der häufig die günstigste Abschät- 
zung liefert. Man erhält dann nach leichter Rechnung 
2f(a) 2Fla)]. 2 F(a) 
a a |” Bi 


da, = (a) + Re O e 


?) Der allgemeinere Fall (4.1) ohne die Forderung (2.5) läßt sich genauso einfach behandeln. 
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Ist &, > 0 und gilt (4.1), so ergibt sich genau wie oben mit G(x) = ii 9(a, x) de: 
Ist im ganzen IntervallO <x < a stets i 
EERRENOIÄTER EA EN ee (6.2), 


so gilt dort (x) < i(x) bzw. tx) > i(x). 
Mit (6.2) ergeben sich analog (3.7) bis (3. 10) und (5.2) bis (5.4) die Bedingungen?) 


A Sa = 2 lo nn = Olahzusn ea Data (6.3), 
ee Ih ENDEN R) a (6.4), 
x x a a 
ta) e2 ul) > 2 Kr) O9 zw ee een (6.5), 
FE) zZ ÜTEHZWENT GE) ORTE (0:0) 
/(&) mit wachsendem x nicht zunehmend bzw. nicht abnehmend . . . . (6.7). 


Bei der Anfangswertaufgabe (1.1), x(0) = 0, x(0) = x, mit negativer Anfangssteigung 
%, < 0 müssen die Ungleichheitszeichen zum Teil vertauscht werden. Da sich dieser Fall jedoch 
durch die Transformation z(ft) = — x(f) auf den obigen zurückführen läßt, soll er hier nicht mehr 
gesondert behandelt werden. 


x 
7. Wir fassen die Ergebnisse in einer Tabelle zusammen. Dabei ist wieder F(x) = [ f(«) dx. 
0 


Die Bedingungen für f(x), von denen in jeder Rubrik jeweils nur eine erfüllt zu sein braucht, 
finden sich oben unter den Gleichungen (3.6)— (3.10), (5.2)— (5.4) und (6.3)—(6.7). Die Zeichen 
s und > in den 3 Spalten beziehen sich jeweils ohne bzw. mit Klammer aufeinander. 


Abschätzungsformeln für die Lösung dee Problems 
<& + f(x) = 0, x(0) = 0, &(0) = ©, > 0, mit2 F(e)< afür0o se<a 


Falls 2 F(a) = S 
Bedingungen für f(x) in [0, «] Schranken für «{t) re na fx) an 
ist7=47] 
Il. 22Fl(a) — a F(x) <0(=0) EN 2 Fa) 
oder x f(x) —2 F(x) <0 (> 0) zu) < "sinoto-1 = 
oder z f(x) — /«) =<0(=0) 63) 
oder f(0) 20 (< 0) und f’(x) < 0 oder ae 
(= 0) a VZ2F@) 2/2 F 
oder lim f(x) >0 (< 0) und x{) =, sinh #t,u = —— (S) 
zorr0 
f(x) nicht zunehmend (=) 
(nicht abnehmend) fürO se <a 
Il. aFf(&e) —xF(a) <0(=0) — F(a) 
oder & f(x) — F(«) 20 (< 0) 2 >, < 4% __aV? 
oder re) =0(=0) (<) TER) VF@) 
ee, fir 25a (=) 
ENT: ie) = = x + f(a) <(t)> — coswt + "int 
% 
< 1 Dt x 
ee <0(>0) (2) „_2fla) 2 Fa) ?< _arctan = Brionh 
a np? (>) 
oder 2 F(a) 2 Ha} 1 [7 
2 I rn — 
Te a BI, 02 Fa a) 
oder 22 /(x) —2xf(x) +2 F(x) <0| x(t)> "si Pe Tee Fi a fla 
0) i 22 —= — w? 
oder f(x) <0(=0) (=) f 
oder f’(x) nicht zunehmend ER pn 
(nicht abnehmend) Tun 0rem=dg. 


3) Dabei gilt (6.k), wenn (6.k + 1) erfüllt ist, für = 3,4,5. 
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8. Beispiele: Es sei jetzt (— 2) = — /(&). 

Geht dann F(a)/a mit wachsendem Ausschlag gegen Unendlich (harte Feder), so folgt en 5 
daß die Schwingungsdauer T = 4r gegen Null konvergiert. Ist etwa (0) = 0 und f ee = 
für x > 0, so liefert I eine untere Schranke für die erste Viertelschwingung und eine obere Sc + 
für die Schwingungsdauer T im Abhängigkeit vom Ausschlag a und III eine obere Schranke für 
die erste Viertelschwingung und eine untere für T(a). SA 

Entsprechend kann man im Falle /(0) = 0, f’(x) = 0, durch I und III die Lösung nach 
oben und unten abschätzen. Ist außerdem für zwei Konstanten z, >0 und x, > x, die Funktion 
F(&) > 0 für x, <x < x, so treten für diese Ausschläge periodische Lösungen auf, und die Funk- 
tion T(a) läßt sich durch die in der dritten Spalte des Schemas angegebenen Funktionen nach 
unten (I) und oben (III) leicht abschätzen. 


= LM e ne 
Beispil 1: &+2X% + 5 0, Eh 7 + g Ausschlag «a 


EEE 3 & (is +16 + 3ai 
III); ———2—> — arctan | —————— — — 
I und III®) AT = Are pr 
Speziell erhält man für den Ausschlag a = 1 aus 
—— 7 
I. x() > sin Y1,251 und Sm » 1,405, 
II ea Last und ESEL 
\ rn 16 v5 
A a a en abe. 2 
608 = BIN = > — arctan (— Y35) »# 1,314. 
: II. {l) & ee ana see losen und er n (— 35) 
I und III ergeben also T=4r = 5,438 + 0,182. 
# 1 2 a4 $ 
Beispiel 2: +2 zZ" =0, Foeser | 
a T | 9 | 
I ONE = = -H0<a << —, 
= er Nee v3 
16—16@+3at = 
I und III: 2 < £ = m. arctan Ein . ch für0 <a<py2. 
Aa aan a 


Speziell für den Ausschlag a = 1: 


ed ee eek 
1. aDesin 2% % = v3 „old, 
I: cd) > = 2 + e: Y3t; TS 

: 716 a “ 


111.402 098 j +3 sin : 


also 1,814 sr s 2,094. 


—1, rs2arctan Y3 » 2,094, 


Man hat auch hier bei sehr geringer Rechenarbeit schnell einen Überblick. Die Abschätzun- 
gen I und III werden um so genauer je weniger die Kennlinie von der Geraden abweicht. 


9. Bisher wurde die Kennlinie durch eine einzige Gerade ersetzt. Genauere Ergebnisse 
erhält man, wenn sie durch einen Geradenzug ersetzt wird. K. KLOTTEr!) verwendet dafür Poly- 
gonzüge, die stets unterhalb oder oberhalb der Kennlinie liegen, um untere und obere Schranken 
zu erhalten. Es genügt jedoch auch hier nach 2. und 4., die Polygonzüge so zu wählen, daß ihr 
Integral stets oberhalb bzw. unterhalb F(x) bleibt, um Schranken zu bekommen. 

Der einfachste Fall ist ein Treppenzug, bei dem in den einzelnen Intervallen die Vergleichs- 
kennlinien Konstante sind. Analog den Verfahren der mittleren Ordinaten und der mittleren 
Abszissen hat man im ersten Fall im Intervall &,_, <x = x, für die Ordinate c, = AF,/Ar, zu 
setzen und im zweiten Fall mittlere Abszissen x, passend einzuführen. Ist etwa /(&) monoton 
nicht fallend für x = 0, so erhält man bei den „mittleren Ordinaten“ durch Integration eine 


4) Hier ist aretan (— y) = m — arctan y zu setzen. 
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‚obere, bei den ‚‚mittleren Abszissen‘‘ eine untere Schranke für (x) und r. Im ersten Fall gilt dann 
im Intervall [x,_,, x,] das Vergleichsproblem 


i= — Cy <(t,—ı) = %ı—1 en Ti eg also 2 — 1 ei (t > l,—ı) und 
—— Cy 


2 


I ey] = 


(F— a) Hr (24 + Cy t,—_ı) (==4,27) . 


Insbesondere ist&,—= Yy&3 — 2 F,nach Wahl der c, und At, = — and 


v 
Bei den „mittleren Abszissen‘‘ ergeben sich mit /(z,) = f, die Zwischenabszissen aus 


Az, fh) nn 
Af, 


und insbesondere 2, = 2, = (0, %n+ı = x,. Im Intervall [%,_1, %,] ( = 1,2,...) haben wir dann 
das Vergleichsproblem 


AF, = Fi) — Fo) = —%—)) h-ı j“ (.—%) fr, also, = 


EN RN N Br VE 


z v 
Ist g(x) die Ersatzkennlinie (Treppenfunktion) und G, = [ . ga) de= Dh u_), so 
erhalten wir mit /(0) = 0 0 u=1 


im ersten Intervall [0, x] & = &,- 1, also t; = 4(@) =" (&, + 0) und 
Top 


im v-ten Intervall %,_,,%, = 23,3,...,n +1, &— u _ı=—h_ı(lt—t,_ı)und 
t— 1 = ei Bei) ee rear 
Insbesondere gilt AG, = h_1ı u — 1) = h_14% x, = a G, und At, = 7, 2 (HER) 
— & v—1 


und = “1 für h=% 
Tg 

Dieses einfache Vergleichsverfahren soll zur Illustration an dem Beispiel der Differential- 

gleichung £ + x = 0 mit 1, = 2, also dem Ausschlag a = 2 durchgeführt werden. Man erhält 


hier für die obere und untere Schranke von f(x) und r = = die nach den obigen Erklärungen 
leicht verständlichen Schemas A und B: 


A. 
Ly Axy Fy AFy GY= AF, % — Az Abo ty 
Axy 
o=0 u 2 0 
1 ua 0,5 En 06705 0,5350 
einge | 1,73205 0,53590 
1 ERBEN 1,5 | 1,5 1,73205 1,15470 
2 2 - N) 1,69060 


| Au | fr | A | mh | Altufo) | Fo IAF)| ©, | Atv | AG, | @% x» | — Ag Ah ty 
| | 

2. er ae We, RER BEN ER ge u 0 
0 0 EIER a ES a 

2 a N ER I U | 21001057 bee. 2 E a 
1 ı 0,5 A nme _ 
we 1 3 LS TTer 1) 1,41421 ——| 0,83579 
2) 2 ı— 4 —| 2 10,5| 1. I! 1,41421 | 0,7071 — 
— — — — 2 19 0 ——| 1,54290 


Bei dieser groben Einteilung von [0, a] in nur zwei Intervalle erhalten wir also für —=n [2 
die Abschätzung 1,5429 <r = 1,6906. Bei Verfeinerung der Einteilung konvergieren die Schran- 
ken gegen den wahren Wert. 
Manuskripteingang: 12. 1. 1961 

Anschrift: Dozent Dr. Hans Eurmann, Mathematisches Institut A der Technischen Hochschule Stuttgart 
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Über die Lösung linearer Integralgleichungen 
mit einer elektronischen Analogierechenanlage”) 


Von KARL-WALTER GAEDE 


Die Klasse der durch eine elektronische Analogierechenanlage lösbaren (schaltbaren) linearen Integral- 
gleichungen zweiter Art wird durch notwendige und hinreichende Bedingungen für die freie Funktion und 
den Kern charakterisiert, und es wird die Lösungsmethode angegeben. ee nes 

Nach Betrachtung von Spezialfällen werden für den Verfahrensfehler, der durch Approximation 
einer gegebenen Integralgleichung durch eine schaltbare entsteht, mit der Analogierechenanlage bestimmbare 


scharfe Schranken gegeben. b 
Die Lösung linearer Integrodifferentialgleichungen wird kurz behandelt. 


The class of linear integral equations of second kind that may be solved by electronic analog computers 
can be characterised by necessary and sufficient conditions for the frde function and the nucleus. The method 
of solving the equation is given. ‚ | , 

After a discussion of special cases, the problem of determining the error involved when a given integral 
equation is approximated by another one soluble in the sense explained is likewise solved by given definite 
limits that can be determined by the analog computer. 

The problem of solving linear integro-differential equations is also shortly discussed. 


Heo6xoAHMbIMH MU HOCTATOYHBIMH YCJIOBUAMH HA CBOOOAHBIÄ YJIEH H ANLPO XapakTepus3yeTch 
Klacc JINHeÜHBIX MHTETPAJIBHbIX YPABHeHHÄÜ BTOPOTO PONAa, KOTOPble MOTYT ÖbITb PellleHbl 
INEKTPOHHBIM BbIYHCHHTEJIBHBIM YCTPOÜCTBOM HEIpepbIBHOrO Neictsun. Jlaötca meron pe- 
IIIeHHA. 

PaccMaTpuBalwTcH CIIelMAalbHbIe cJIyyan. [IA HOTPeIIHOCTH, BbIBBAHHOÄ AuUmpOKCHMalnmen 
AAHHOTO HHTETPAAIBHOTO YPABHeHHA TAKUM, KOTOPoOeE pemmaeMO Ha 3lIEKTPOHHOM BbIYNCJIHTEJIb- 
HOM YCTPoüÜCTBe HelpepbIBHOTO MefCTBuA, NAI0TCA TOYHBbIe TPAHHLEI, KOTOPbIe MOTYT ÖbITb 
onpenesleHbI 3JIEKTPOHHBIM BbIYHCJIMTENBHBLIM YCTPOÜCTBOM HEHUPEPBIBHOTO HeliCTBuA. 

PaccmaTpnuBaeTtca pemeHue JIHHeÜHBIX UHTETPOAHPPepeHNMAalbHbIX YPaBHeHHü. 


Einleitung, Problemstellung 


1. Seit etwa 1950 haben elektronische Analogrechner verhältnismäßig starke Verbreitung 
gefunden. Es handelt sich hierbei um eine mit einer gemeinsamen Steuerungseinrichtung ver- 
sehene Zusammenstellung von einzelnen Analog-Rechenelementen. 

Den Rechengrößen sind als Analoga zeitlich veränderliche elektrische Spannungen zuge- 
ordnet. Die Rechenelemente erzeugen, gesteuert von den in sie hineingesteckten Spannungen, 
Ausgangsspannungen, so daß die Maßzahlen von Eingangsspannungen und Ausgangsspannung 
bestimmte mathematische Relationen erfüllen. 


In den handelsüblichen Analogrechenanlagen sind folgende Rechenelemente zu finden: 

(Alle Rechengrößen x, x, sind Funktionen der unabhängigen Variabeln s, die wir „Ma- 
schinenzeit‘“ nennen wollen.) 

1. Summatoren: 


N 
Eingänge: 2, ..., %,, Ausgang: = — Y k,x,, wobei k, > 1 bestimmte (von Gerät zu 


v=1 


Gerät verschiedene) ganze Zahlen sind. 
2. Potentiometer: Eingang x,; Ausgang = cx,;c= const; Oo <ce<i1. 


. ur j Me 
3. Integratoren: Eingänge x, ..., 2,; Ausgang x = — (% SIEB dt) 
Ov=1 
k, wie unter 1.; a, ist ein am Integrator einstellbarer Anfangswert. Es steht für alle 
Integratoren nur dieselbe Integrationsveränderliche zur Verfügung. 
4. Multiplikatoren: Eingänge x,, 23; Ausgang x = x, : 2. 
Durch Zusätze oder Umschaltung auch zur Division verwendbar: 
Eingänge 7, %,; Ausgang x = x /x, (x, + 0). 
5. Funktionsgeber: Eingang x; Ausgang x — (<); f fest vorgewählt; 
und eventuell 
6. Zeitverzögerer: Eingang x,(s); Ausgang x = 2% ( —T); tr > 0 fest vorwählbar. 


Den weiteren Betrachtungen legen wir eine Analogrechenanlage zugrunde, die Rechen- 
elemente 1. bis 5. (in genügender Zahl) enthält. 


ee eh BEER 3 : . 
Von der Fakultät für Allgemeine Wissenschaften der Technischen Hochschule München angenom- 


mene Habilitationsschrift. 1. Referat: Professor Dr ö) 5 
a : Professor Dr. Joser HEINHoLD, 2. Referat: Professor Dr. ROBERT 


—_ 


“ur ee 4 ee 


ud a Zr län ah Lzr NETTER WESDER. 
b ’ ai 2 Kl s 
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Sollen auch Zeitverzögerer zugelassen sein, so wird besonders darauf hingewiesen. 
Mit Hilfe einer derartigen Anlage ist man durch geeignetes Zusammenschalten der einzelnen 
Rechenelemente („Rechenschaltung‘“) in der Lage, Anfangswertprobleme für Systeme von ge- 


wöhnlichen (linearen oder nichtlinearen) Differentialgleichungen bequem zu lösen. 


Weiter lassen sich unter gewissen Bedingungen für die Koeffizientenmatrix lineare alge- 
braische Gleichungssysteme lösen. ; 

Wie bei allen Analogrechengeräten müssen die Rechengrößen im betrachteten Intervall der 
unabhängigen Variabeln s beschränkt sein. Bezüglich näherer Angaben über den Aufbau der 
Rechenelemente und ihre Benutzung in Rechenschaltungen siehe man etwa [1], [2], [3]3). 


2. Versucht man, mit einem Analogrechner des oben festgelegten Typs eine Integralgleichung 
gs) = fs) + J K(s, ti) g(t) dt (K und freell) (SV) Ed) 


zu lösen, so stößt man auf zwei Schwierigkeiten. Erstens steht kein Funktionsgeber für Funk- 
tionen von zwei unabhängigen Variabeln zur Verfügung, und hätte man ihn, so könnte man i. a. 
das Integral in (1) nur für fest vorgewählte diskrete Werte von s berechnen und würde damit g(s) 
sicher nicht kontinuierlich als Funktion der Maschinenzeit erzeugen können. 

Da aber die kontinuierliche Erzeugung gesuchter Funktionen das dem elektronischen 
Analogrechner des obigen Typs am besten angepaßte Aufgabengebiet ist, wird man sagen, 
(1) ist mit ihm gelöst, wenn g(s) kontinuierlich erzeugbar ist. Wir führen deshalb folgende 
Bezeichnung ein: 

Wenn mit Hilfe einer elektronischen Analogrechenanlage eine Lösung g(s) von (1) konti- 

nuierlich als Funktion der Maschinenzeit s erzeugt werden kann, nennen wir (1) „schalt- 

bar“. 

Bisher wurden Integralgleichungen (1) nicht auf ihre Schaltbarkeit untersucht, sondern man 
hat für ihre Lösung mit dem Analogrechenprinzip besondere Apparaturen konstruiert. (Z.B. [4]). 

Ihre wesentlichen Bauelemente sind Funktionsspeicher für Funktionen von zwei, bzw. einer 
unabhängigen Veränderlichen. Dabei sind die gespeicherten Funktionen Stufenfunktionen, die 
die gegebenen, bzw. während der Rechnung erzeugten Funktionen approximieren. Mit Hilfe 
von Rechenelementen 1. bis 5. und während des Rechenablaufes automatisch gesteuerten Schal- 
tern wird nun (1) mittels der üblichen oder leicht modifizierten Iterationsverfahren gelöst. (Im 
abgeänderten Iterationsverfahren geht man entsprechend dem sogenannten Einzelschrittver- 
fahren für lineare algebraische Gleichungen vor.) 

Man benutzt also eine Methode der Digital-Rechentechnik und verwendet dazu Analogie- 
Rechenhilfsmittel. 

Da die Berechnung des Integrals für ein festes s, sehr schnell erfolgen kann, ist es möglich, 
das Ergebnis der Integraltransformation (g(s) aus einem Funktionsgeber entnommen) 


Tg: =) — J Ks, 1) g() di 


auf dem Schirm eines Kathodenstrahloszillographen sichtbar zu machen. Durch Vergleich von 
Tg mit der vorgegebenen Funktion /(s) läßt sich durch systematisches Probieren g(s) so verändern, 
daß Tg = f wird und damit eine Lösung g(s) von (1) gefunden ist. (Siehe [6]; Zusammenfassung 
von [4], [5], [6] in [7]). 

3. Wir werden im Folgenden den Typ aller schaltbaren Integralgleichungen (1) angeben 
und zeigen, wie mit Hilfe einer elektronischen Analogrechenanlage die Lösung zu gewinnen ist, 

Approximiert man eine Integralgleichung durch eine andere des schaltbaren Typs und 
löst letztere, so werden für den Betrag der Differenz zwischen der Lösung der schaltbaren und 
der durch sie approximierten Integralgleichung scharfe obere Schranken angegeben, die wiederum 
mit Hilfe einer elektronischen Analogrechenanlage ermittelt werden können. 

Zum Abschluß wird auf die Behandlung von Integrodifferentialgleichungen der Art 


nlo) = K9) + J KG, Yalg) di ee ee 2) 


eingegangen, wobei A, und A, lineare Differentialoperatoren sind und die Ordnung von 4, nicht 
kleiner als die von 2, ist. 

4. Da wir im Folgenden nicht auf die sog. Dimensionierung eingehen werden, sollen 
gegebenenfalls — ohne daß wir dies jeweils besonders erwähnen — in den Rechenschaltungen 
nicht die in der Formel auftretenden Rechengrößen selbst, sondern nur ihnen proportionale 
Größen (mit geeigneten Proportionalitätskonstanten) vorkommen. Ebenso braucht, falls nötig, 
die verwendete Maschinenzeit nicht gleich, sondern nur proportional s zu sein. 


1) Zahlen in eckigen Klammern beziehen sich auf das Literaturverzeichnis am Schluß. 
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$ 1. Notwendige Bedingungungen für die Schaltbarkeit Hr 
i ä j der Kern K(s, !) und die freie 

1. Wir werden zunächst untersuchen, welche Bedingungen 
Funktion f(s) von (1) notwendig erfüllen müssen, damit (1) schaltbar ist. (Alle auftretenden 


Funktionen seien reell.) . 
Zur bequemeren Schreibweise führen wir folgende Begriffe ein. 
Sei /\ die Klasse der Funktionen, die in [0, a] stückweise glatt sind. x 
Da die in einer Rechenschaltung auftretenden Funktionen zu obiger Klasse gehören 


müssen, folgt 6) Ar ee (1.a). 


Mit I, bezeichnen wir die Klasse derjenigen Funktionen von zwei Veränderlichen, für 
die aus Rel‘, und hel, folgt 


[RsOhOdeN Le re 32 Er 
u u 
Durch (1.1) ist eine lineare Transformation von I, in sich gegeben, die wir abkürzend 
it 
E Rh... vun SD WEI Bu 
bezeichnen. 


Als Rechenregeln werden definiert (siehe z. B. [9] S. 137) c Wh =e Rh; (R, + N) h 
=NALRARN,) = N, (N, h) und als Abkürzung (NR... RNh=NRh. 


a-mal 
Wir werden den zu einer Funktion aus /‘, (die wir immer mit großen lateinischen Buch- 
staben benennen) gehörigen Operator mit demselben deutschen Buchstaben bezeichnen. 
Wird der Operator vor eine Funktion aus J', geschrieben, so soll er auf ihr erstes Argument 
wirken. 
Dann schreibt sich (1) 


g=LHt RE: re 


Da die Lösung g(s) von (1) als Funktion der Maschinenzeit s kontinuierlich erzeugt werden 
soll, muß gelten 
Ks, De I SR 2 


Welche Forderungen müssen nun noch an K gestellt werden, damit (1.3) schaltbar wird ? 

Wie in der Einleitung unter 1. ausgeführt, können die Integratoren nur folgende Aufgaben 
lösen: - 

Zu einer Funktion h(s) e /) der Maschinenzeit, die als Eingangsgröße für den Integrator 
zur Verfügung steht, wird an seinem Ausgang die Funktion 


u) = — (00 + Jh) a) (a, = const) . 22 Fe BEE 


erzeugt. 

Wesentlich hieran ist, daß der Integrand als Funktion der Maschinenzeit zur Verfügung 
stehen muß. Hiervon kann maschinengemäß nicht abgegangen werden. 

In dieser Richtung ist also eine Verallgemeinerung des Integrationsprozesses nicht erreich- 
bar. Es sind höchstens noch die Integrationsgrenzen abänderbar. 

Soll also 


vs 


(s)= no dt 


us) 
berechnet werden, so muß nach (1.4) als Eingangsgröße 


2(s) = hu(s)) v’(J)— hluls)) ul) - » . 2. oo 2, (1.5) 
zur Verfügung stehen; d.h. man benötigt die vier rechts stehenden Funktionen. 
Sind die Grenzen Konstante, so ist z/(s) = 0. Dann kann man aber 


b 
Shit) dt 0zsa<b 


bestimmen, indem man zu den Maschinenzeitpunkten a und b den Wert des Integrals fh dt 
0 


abliest; man braucht also nur A(s) als Eingangsgröße. Aus Vorstehendem folgt: 
Damit g kontinuierlich erzeugbar ist, kann ® g nur aus Bestandteilen der Form 


vs) 
k(s) J Kt) g(t) dt; Kl,u, ve]. 2 (l.c) 
us) 


linear kombiniert sein; (weil (1.3) auch linear ist). 


A x x . = 
. 
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Die notwendigen Bedingungen (1.a)—(1.c) wurden aus der Wirk j i 
_ Rechenelemente gewonnen. = Buy SS 
Pe Nun muß noch untersucht werden, ob das Zusammenschalten der einzelnen Rechen- 
elemente zu einer Rechenschaltung nach dem Rückkopplungsprinzip noch weitere Beschrän- 
kungen nötig macht. ET 
x Bis auf unwesentliche Modifikationen und Bezeichnungsänderungen muß die Rechen- 
schaltung wie in Bild 1 aussehen. Man stellt sich auf den Standpunkt, g(s) schon in der Ma- 
Rechenelementen f+ g erzeugt, wie in Bild 1 
durch das Rechteck symbolisiert. 
oFfls)+Rg 


. soll darin bedeuten, daß an die durch den 99) i 
Punkt freigehaltenen Stelle die Funktion am Ein- 
gang der Rechenschaltung zu setzen ist. Eingang 


und Ausgang der Rechenschaltung sollen nach (1.3) II 
- gleich sein. Dies erzwingt man durch Rückkopplung Bild 
111,1. 


Zur Erzeugung von g müßten beim Auftreten von Gliedern (1.c) nach (1.5) auch g(v(s 
und g(u(s)) aus g(s) erzeugt werden können. ne Lo 

Dies ist mit den Rechenelementen 1.)—5.) aber nicht möglich. Aus (1.5) erkennen wir 
weiter, daß die Grenzen von (l1.c) nur lauten dürfen 


u(s) = s,v(s) = const; oder u(s) = const, v(s)=s; oder u(s) = const, v(s) = const (1.d). 
Zusammenfassend erhalten wir: 

Damit (1.3) schaltbar wird, müssen notwendig die Bedingungen (1.a)—(1.c) erfüllt sein 
‘und die Integrationsgrenzen in (1.c) durch (1.d) gegeben sein. 

2. Aus (1.c) und (1.d) folgt für die Bauart des Kerns K(s, t): 

Man zerlege das Quadrat B= {(s,4)|0<s,t<a} in (endlich viele) sich nicht über- 
lappende Bereiche ®,, so daß % Flächen (8,) = Fläche (B) und die Ränder der B, nur aus 


u 

Strecken auf Geraden s= const, f= const, {=s bestehen. K(s, f) ist dann gegeben durch 
DER = ‚ N 

(B) K(s, en) nt) für (Sl, Mae ee NO) 


und k, leT', (damit gehört K(s, t) natürlich zu I%). ((1.6) ist Folge von Bedingung (1.c), die 
Art der Zerlegung von ® ergibt sich aus (1.d.)) 
Eine andere Bauart darf der Kern K(s, ft) nicht haben, wenn (1.3) schaltbar sein soll. 
Damit haben wir folgenden 
Satz 1: Damit die Lösung der Integralgleichung 


96) = IK) + f K(s, da) di 


kontinuierlich als Funktion der Maschinenzeit s mit einem Analogrechner des in der 
Einleitung festgelegten Typs erzeugt werden kann, „damit die Integralgleichung 
schaltbar ist“, ist notwendig, daß f(s) stückweise glatt ist und K(s, f) den obigen Be- 
dingungen (B) genügt. | 
Es fragt sich nun, ob die in Satz 1 angegebenen Bedingungen auch hinreichend für die 
Schaltbarkeit der Integralgleichung sind. 
Diese Frage beantworten wir im nächsten Paragraphen, indem wir angeben, wie man die 
Lösung der in Satz 1 charakterisierten Integralgleichung mit dem Analogrechner ermittelt. 


$ 3. Schaltung der Integralgleiehung 


Wir formulieren das Ziel dieses Paragraphen zunächst als 

Satz 2: Die in Satz 1 charakterisierten Integralgleichungen sind schaltbar, d. h. ihre Lösung 
läßt sich mit einem Analogrechner des betrachteten Typs kontinuierlich als Funktion 
der Maschinenzeit s erzeugen. 

1. Der Beweis erfolgt konstruktiv, indem wir ein Lösungsverfahren angeben. 

Um die Darstellung zu vereinheitlichen, verfeinern wir die Zerlegung von ® gemäß Satz 1 
durch Einziehen neuer Trennlinien auf Geraden f = const. und s = const. solange, bis ® von 
einem Netz aus achsenparallelen Rechtecken und rechtwinkligen Dreiecken mit der Hypothenuse 
auf {= s und den Katheten parallel zu den Koordinatenachsen überzogen wird. 

Die hier angegebene Netzverfeinerung ist stets möglich, da es nach Voraussetzung nur 


endlich viele ®, gibt. 


schine zur Verfügung zu haben. Daraus wird mit Hilfe einer Kombination von geeigneten 
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Die Teilbereiche ®,, die wir auf diese Weise erhalten, numerieren wir so, daß die Dragch? 
mit Scheitel in s > (ihre Anzahl sei m) ®, (u = 1,..., m,) heißen. Wir setzen m, = 2 m, 


0e-> 

und nennen die Dreiecke mit dem Scheitel in s< f) (Anzahl m) 9, a =1-+m,»...; m,) und 
die Rechtecke (Anzahl m;) ®, (u=1-+Mmy...,M;). 

Nach Satz 1 ist in jedem Teilbereich ®, (u =1,...,m + m; + m;) K(5, if) durch eine 
Summe der Form ; k,(s) I,(l) gegeben. Die Anzahl der Summanden für ®, sei n,. en 

Wir führen die Indexfunktion g(v) ein, die festlegt, welches der Index des Teilbereiches 
9, ist, in dem k, bzw. I, an der Definition von K(s, t) beteiligt ist. 

Wir setzen dazu 


et: u z 
eb) = u für IH 2nSsy<Yn  Vere (2.1) 
e-1 De 


(v = HE) 


und verwenden zur Festlegung von Indexmengen falls Verwechslungen ausgeschlossen sind, statt 


oo) = u} einfach op 2, Be (2.1c). 
Mit diesen Bezeichnungen legen wir fest: 
Kid =, 2 Kiss) für (,DEB, Due (2.1d). 
Yb)=u 


Wir können ohne Beschränkung der Allgemeinheit annehmen, daß die k,(s) für jeden Teilbereich 
linear unabhängig sind, ebenso die /,(t). 

Um in der gegebenen Integralgleichung die Integrationsgrenzen angeben zu können, 
nennen wir 


Ss, = up le Ss, mi (), 
(DER, HC, 
EEE TEN - a2 
= sup(d), = uf © 
(EB, UEB, 


und setzen fest, daß k,(s) = 0 für s € [s,g#)» Sogw]- Dann heißt die Integralgleichung (1.3) 


mı s 7 top) Ms !og() 
=) + I ke) SEOWd+ Z ke) SLOgWd+ 8 Ks) SL) gddt (23). 
en=1l Lupe) P)=1+m, & eW=1+m,  fupl) 

Die erste Summe rührt von den Dreiecksbereichen mit Scheitel in s >? her, die zweite 
von den Dreiecksbereichen mit Scheitel in s< ft, die dritte Summe kommt von den Rechteck- 
bereichen. 

Wir führen jetzt folgende Funktionen und Konstanten ein 


h,(s) = SL gl) dt für seh) —l..., nur Renee (2.4a), 
up) 
logo) 
hs) = — S LO gl) dt für vsebeoß)=1-+m,...;M} » . . (2Ab), 
Kr0) 
= / Lg) dt für vepbeobß)=1+tm...,M} ... (2.4c). 


Fupl) 


Aus (2.4a) und (2.4b) folgt: (* = Ableitung nach s) 


R(s) = L(s) 96) b ne n) rt a (25). 
o=1 
Aus (2.5) zusammen mit (2.3) ergibt sich nach Einsetzen von (2.4) ein lineares Differential- 
gleichungssystem für die h,, das sich durch Einführung der Spaltenvektoren mit 


n= 3 n, Zeilen 
o=1 
h(s) = {h,(s)} ; I(s) = {1,(s)} ee (2.6a) 
und der quadratischen (n, rn) Matrix 
A(s) = {+ 1(s) k;(s)} (+ für Tel de 5 N. Ar et 3 Ne 3: n,) (2.6b) 
D=l = —_ 
ergibt zu $ en ra PR 
K)=AJKIHSÜÄ)L+. 3 KU) ..... AN (a7 


Pe)=1-+m, 


a de u nd ud Du U m ern u u 


Aus @. 4) folgen für ve) noch eheneiigdagen, wobei zu beachten ist, daß in den Drei- 
En erschen 5 2 Sy. und & „= So „Ist, was in (2. 4a) und (2.4b) eingesetzt, liefert 


AlSuoo)) = für 2 20) =1,...,m}, | 
Also) = 0 für vebloß@)=1-+m,...,m,} 
Damit haben wir die Lösung der in Satz 1 beschriebenen Integralgleichung als folgendem 
Problem äquivalent erkannt: 


Man löse das Differentialgleichungssystem @. 7) unter den Bedingungen (2.8), wobei die 
in (2.7) noch nicht bestimmten Konstanten c, #)p@) =1 + m,...,m,;) den Donau 
(2.4c) genügen müssen. 

Unter einer Lösung von (2. 2): wollen wir dabei folgendes verstehen: 76) bezeichne die 


rechte Seite von (2.7), dann heißt : h(s) Lösung von (2.7), wenn es in [0, a] bis auf die (endlich 


vielen) Unstetigkeitsstellen 5, $, von T(s) (die nur von erster Art sein können) differenzierbar ist. 
In den 5, soll hs) stetig sein und rechts- bzw. linksseitige Ableitungen besitzen, die gleich 
7T6&+ 0) bzw. 7 (s,— 0) sind. 

Gibt man an einer Stelle , nO<s<a his.) vor, so existiert aan eine Lösung von. 


- (2.7) in obigem Sinne, die in s, den Wert his.) annimmt. 


(Dies ist wegen der gegebenen ES BCH Ups Vorsalui Folge des bekannten Existenz- 
und Eindeutigkeitssatzes für stetiges 7(s).) 

Diese Definition einer Lösung wurde gewählt, weil bei der Schaltung eines Anfangwert- 
problems auf der elektrischen Analogierechenanlage (sofern die c, bekannt sind) un diese 
Lösung von (2.7) erzeugt wird. 

Sehen wir die Integralgleichung aus Satz als gelöst an, wenn eine stückweise glatte 
Funktion g(s) gefunden ist, die die Integralgleichung erfüllt, so ist jede Lösung von (2.7), (2.8), 
(2.4c) vermittels 


Ma Ms 
a) = 1 + ” KOhO— 2 ki) + % ak6) 
Pb)=14+m, PW)=1-+m;, 
auch Lösung der ER es Umgekehrt erfüllt jede derartige Lösung der Integral- 
gleichung auch (2.4c), (2.7), (2.8), wie aus der Herleitung dieser drei Formeln folgt. 

2. Nach dieser Umformung des gestellten Problems bleibt noch zu zeigen, wie man das 
Differentialgleichungssystem unter den gegebenen Nebenbedingungen mit der elektronischen 
Analogierechenanlage lösen kann. 

Direkt schalten und lösen läßt sich ein solches System nur, wenn Anfangsbedingungen 
vorgegeben sind. Wir führen nun unser Problem auf die Lösung von Anfangswertaufgaben und 
einem System linearer algebraischer Gleichungen zurück. 

Wir erzeugen uns zunächst n linear unabhängige Lösungen des homogenen Differential- 
gleichungssystems P Beh 

ET EN ee (2.9), 


hy); Ra); 2.» 5 N3(5) 
bezeichnen. Dies erfolgt etwa so, daß man als h,(s) die Lösung von (2.9) in [O, a] betrachtet, 
für die A,(0) gleich dem »-ten Einheitsvektor ist. 

Damit sind die hs) mit Hilfe des Analogrechners als Lösung eines Anfangswertproblems 


ohne weiteres berechenbar und als wohlbekannt anzusehen. 
Sodann bestimmen wir, etwa indem wir als Anfangsvektor den Nullvektor wählen, eine 


partikuläre Lösung des Systems 


die wir der Reihe nach mit 


Ri) = Ac) hs) + I) Ko) 
die d(s) heißen möge und weiter je eine partikuläre Lösung der Systeme 
Ds) = A(s) h(s) + k,(s) Us) für veblob)=1+Mm,...,m} . . (210), 
die wir h,() (v wie in (2.10)) nennen. 


Nach dem Überlagerungsprinzip läßt sich jede Lösung h(s) von (2.7) aus den hier berech- 
neten Funktionen darstellen als 


— N — — 
Keys D his) Fdlsye = „2. ee RN 
ul 
wobei a 
N=2n 
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Wir setzen zunächst (2.11) und (2.4c) in (2.3) ein und erhalten für g(s): 


96) = /6) + 3 c, Ks) h,(s) + 5 kJ +E)AU).. - - - - (2.12) 
‚1 2 


o(#)=1-+m; 


wobei k(s) der n-zeilige Vektor 


{+ k(s)} 


ist, und T ‚‚transponiert‘ bedeutet. 
Dieser Ausdruck ergibt in (2.4e): 


+ für vepeb)=1,...,M} 
für veßleß) =1+ My..., Mg} 


N ed 5 m; logo) togto) u ) 
s=3co/ KORG 3 ET Kar [od + Kr dd] 1, dt 


1 Eugplo) eW)=1+M, fupto) fupte) 
für oefop@) =1+M,...,m} > 


Aus (2.8) erhalten wir nach Einsetzen von (2.11) noch die folgenden Bedingungsgleichungen 
für die c 


(2.132). 


N 
ex Cr N EN u d,(S4 90) = 0 für vE vo) = r ech m,} r 
a (2.13b), 


N ” - 
a Ug h,(50£6)) + d,(Sog)) N) für ve vle) 1 Mur. m;} 
o=1 


wobei h,, die »-te Komponente von h, ist. 

(2.13a, b) ist ein lineares algebraisches Gleichungssystem mit N Gleichungen für die N 
Unbekannten c,. Die Koeffizienten des Systems sind alle mit Hilfe des Analogrechners bestimm- 
bar. Die Lösungsmannigfaltigkeit von (2.13) bestimmt das Lösungsverhalten von (1.3): Man 
ersieht, ob eine Lösung von (1.3) existiert und von welcher Dimension die vollständige Lösungs- 
mannigfaltigkeit ist. 

Zur Übersicht fassen wir die Schritte, die zur Lösung von (1.3) führen zusammen: 

a) Zerlegung von B in die Dreiecks- und Rechtecksbereiche 

b) Aufstellung des Differentialgleichungssystems (2.7) mit den Nebenbedingungen (2.8) | 

und (2.4c) 

c) Ermittlung von N + 1 Lösungen von (2.7), so daß die vollständige Lösung von (2.7) 

sich in der Form (2.11) schreiben läßt. 

d) Aufstellung des linearen Gleichungssystems (2.13) für die in (2.11) noch unbestimmten c,. 

e) Bestimmung der c, aus (2.13) und Einsetzen in (2.12) liefert die gesuchte Lösung. 


Anmerkungen: 
a) Das Einsetzen der c, gemäß e) geschieht praktisch in folgender Weise: 
Man setzt in (2.7) zunächst die nun bekannten c,(|p@) =1 + m,...,m,) ein. Dann 
sind in (2.7) keine freien Parameter mehr vorhanden. Mit Hilfe der übrigen c, sind für (2.7) 
— N _ .. 
Anfangswerte vorgegeben. Es war nämlich nach (2.11) A(s) = Sc, A,(s) + d(s) und außer- 
ol 


dem 1,(0) gleich dem o-ten Einheitsvektor. Wählen wir, um etwas bestimmtes vor Augen zu 
haben, h,(0) = 0 für W)g@) =1 + my...,ms) und d(0) = 0, so ergibt sich 


Ma 
KO)= 3 6,8,=tl=1ch (Vektor mit n Komponenten) 

po)=1 
(mit €, = v-ten Einheitsvektor). D.h. wir müssen das Anfangswertproblem (2.7) mit h(0) = € 
lösen und erhalten in der Schaltung (in der, wie man in (2.5) sieht, g(s) explizit auftritt) gleich- 
zeitig die gesuchte Lösungsfunktion g(s). 
b) Lassen ‚sich die , aus (2.13) nicht eindeutig bestimmen, so sei die Maximalzahl der 
ER a Lösung des zu (2.13) gehörigen homogenen Problems r und Cıo Cags +++, CNo 
ER I: seien dav on r linear unabhängige Lösungen. Co Ca +++, Cno Sei irgend eine 
Lösung von (2.13). Dann ist mit b, = 1 und b, beliebig fest wählbar, in 


r N Mi r N Ma r Ms 
KEIL I IK — Ed DC % KA) + u 2 0.Kkla) 
o=(Ü o=1 o()=1 0=0 o=1 eb)=1+m, o=0 ?G)=1+m, 


jede Lösung von (1.3) enthalten. 
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4. Ohne den gesamten Formelapparat herzuleiten, sei noch angegeben, welche Klasse von 
Integralgleichungen schaltbar ist, wenn im Analogrechner Zeitverzögerer (Rechenelement 6) der 
Einleitung) benutzbar sind. 

Dann dürfen in dem Netz, das ® in Teilgebiete zerlegt, neben den Strecken auf t=s 
und ?= const., s= const. auch noch Strecken auf {=s—r(r >0) vorkommen. Der Kern 
muß in jedem Teilbereich von 3 wieder in der Form 3) k,(s) I,(t) gegeben sein und die Lösung 
erfolgt im Prinzip nach dem oben geschilderten Verfahren. 

5. Da der Aufbau eines schaltbaren Kerns sich nicht ändert, wenn man in K(s, t) s.und t 
vertauscht, folgt, daß mit der Integralgleichung (1.3) auch stets die Integralgleichung 


96) = fl) + J Kies) dr ee ware 1 08 (2.14) 


schaltbar ist. Wegen.der Schaltbarkeit von (2.14) läßt sich der zu K(s, ) inverse Kern, den 
wir H(s, t) nennen, in einer für spätere Zwecke günstigen Form ermitteln. 


Für H(s, t) gilt: 
HGD KH + [Hs Kedde= Kt) + |K,DH6Dd .. @l5) 


Halten wir hierin s fest, so liegt für H(s, t) eine Integralgleichung der Form (2.14) vor, 
die schaltbar ist. Wir erhalten also den zu K(s, t) inversen Kern für festes s als Funktion seiner 
zweiten Variabeln f, wobei hier ausnahmsweise die Maschinenzeit { genannt wurde. 


$ 3. Wichtige Spezialfälle 


1. Der Kern K(s, t) ist entartet: K(s, t) = 3 k,(s) Lt) in ®B. Dann ist in unserer Be- 
v=1 


zeichnungsweiss m =0; m =0; m, =1 und es ist = %,. Damit wird das Differential- 
gleichungssystem (2.7) mit den Bedingungen (2.8) gegenstandslos, da alle h, verschwinden. Von 
(2.13) bleibt nur (2.13a) übrig und es ist /,ı = 0; 4ı =a; N =n,. Damit heißt (2.13a) für die 
betrachtete Spezialisierung: 


SE [EOLDELI Od, 01... ma. ea. 
v1 0 0 


Es ergibt sich also die für eine entartete Integralgleichung übliche Lösungsmethode. 

Wie man sieht, wird die Analogierechenanlage hier nur zur Berechnung der bestimmten 
Integrale in (3.1) verwendet, (evtl. noch zur Lösung von (3.1)). 

Damit werden die Möglichkeiten der Analogierechenanlage nicht voll ausgeschöpft, be- 
sonders, da das Rückkopplungsprinzip nicht verwendet wird. Andererseits können die in (3.1) 
auftretenden Koeffizienten alle nacheinander berechnet werden, so daß nur eine geringe Anzahl 
von Rechenelementen vorhanden zu sein braucht. 

2. Es sei B, das Dreieck s >1;0 <s,t< aund B, das Dreieck s<t; 0ss,tsa. (D.h. 
m=m,=1 m, =). Der Kern ist dann gegeben durch 


Ks) Li) fü  6DdeR, 
v1 


Ks, = 


N a N 


v=1+n, 


Das Differentialgleichungssystem wird zu 


N, NN, 
KEIM EL KIEL) kA Ra 2: 
vol v=1+n 
Da keine Rechteckbereiche existieren, entfallen die Bedingungen (2.13a) und es bleibt 
nur (2.13b) übrig. Hierin ist zu setzen .ı=0; 05 =aundN=n-+n, 
Damit wird aus (2.135) 


n,+tu 


u SE a a re I Wr 
7 74 REte 3:3), 
Nı-+Na 
I ca h,sla) + d,(a) = 0 für vy=1-+n,...M+Nn 
o=1 
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Machen wir noch von der Tatsache Gebrauch, daß h,(0) gleich dem o-ten Einheitsvektor 
ist (siehe $ 2.3) und verabreden wir für unseren Spezialfall, daß 


d(0) = 6 
so ergibt sich aus der ersten Zeile von (3.3) 
ai für Pe 5 FORTE 


Aus der zweiten Zeile von (3.3) wird damit 
SG ohld+d)=0 fü vorinm In 
o=1+n, 


Die in $2 angegebenen Lösungsschritte reduzieren ‚sich hier auf: 

a) Man bestimme die Lösung d(s) von (3.2) für die d(0) — 0 ist. Sodann setze man f = 0 
und bestimme der Reihe nach die Lösungen h,(s) mit h,(0) =&, füro =1-+n,...,Aı +. 
(Man benötigt also keine vollständige Lösung mehr!) Aus den so gefundenen Lösungen lese 
man h,,(a) und d,(a) für» =1+n,...,7",+n, ab, setze diese Zahlen in (3.4) ein und be- 
stimme die „(oe=1-+n,...,Aı +n,). 

Darauf löse man (3.2) mit den Anfangsbedingungen h,(0) =0 für »=|1,...,n, und 
h,(0) = ce, fürv=1+n,...,n,+ng In der Rechenschaltung wird dann die Lösung g(s) 
von (1.3) erzeugt. 

Hier werden die Möglichkeiten einer Analogierechenanlage voll ausgenützt. 


(Bemerkung: Sieht man von der Beziehung zwischen (1.3) und (3.2) ab, so wurde hier 
für die elektronische Analogierechenanlage ein Verfahren zur Lösung eines Systems linearer 
Differentialgleichungen mit Randbedingungen (3.3) angegeben.) 

Eine Verminderung des Bedarfs an Rechenelementen kann oft erfolgen, wenn K(s, i) 


GrEENsche Funktion ist für einen linearen Differentialoperator n-ter Ordnung A[u] mit den 
Randbedingungen 


IM» 


(a; u®=D(0) + b;, u®—D(a)) = 0 (el; (3.42). 
k 


Bekanntlich weist diese GREENsche Funktion die Form des in diesem Abschnitt betrachteten 
Spezialfalls von K(s,f) mitn, =n, =n auf. 


Man setzt in diesem Falle 
u= [Kfs, t) g(t) dt 
0 


und erhält g= / + u und nach der Voraussetzung über K 


Mu=fpu „20, m (3.5). 


Es ist also die lineare Differentialgleichung (3.5) unter den Nebenbedingungen (3.4) zu 
lösen. 


Die Einsparung an Rechenelementen macht sich bei der Anzahl der Multiplikatoren 
bemerkbar. 


Für die Lösung von (1.3) nach (3.2), (3.3) würde man in der Regel 4 n Multiplikatoren 
benötigen, während man für die Lösung von (3.4), (3.5) nur n Multiplikatoren braucht, wobei 


die Multiplikatoren nicht eingerechnet wurden, die evtl. zur Erzeugung der k, und /, erfor- 
derlich sind. 


Es ist nämlich nach (3.5) eine lineare Differentialgleichung n-ter Ordnung zu schalten, 
wozu i. a. n Multiplikatoren gebraucht werden; (wenn vorausgesetzt wird, daß die höchste 
Ableitung den Faktor 1 hat). Für die Schaltung von (3.2) dagegen braucht man 2 Multipli- 
katoren zur Produktbildung /,g und 2n Multiplikatoren zur Erzeugung von 

n 2n 
g=/+ Zk,h— P} a RR 
v=1 


v=1-+n 


Die Lösung von (3.4) und (3.5) erfolgt nach dem in $ 2,3 beschriebenen Schema: 


Man bestimme mit dem Analogrechner eine spezielle Lösung u, von (3.5), indem man 
etwa die Anfangswerte u,(0), Uo(0), ..., UR2(0) gleich Null setzt, und n linear unabhängige Lö- 
sungen U, ...,U, des zu (3.5) gehörigen homogenen Problems. 


— 
EZ 
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Man wählt dazu etwa Ä 
(0) =-.-= u 20)=0, 1; Marsa, 9-0: , V-L....D)% 


Dann hat die gesuchte Lösung sicher die Gestalt u(s) = u,(s) + > c;u,;(s) mit noch zu be- 
je 


stimmenden Konstanten c,. Durch Einsetzen in_(3.4a) erhalten wir 


n 


n n 
(au 3 bu fa) g=— bu) für dmheın.. 66). 


a 


Das sind wieder n lineare Gleichungen für die c;, deren Koeffizienten bekannt sind, da 
die u,(s), u,(8),...., u,(s) mit Hilfe des Analogrechners berechnet wurden. 

g9(s) erhält man, indem man (3.5) mit den aus (3.6) berechneten c; als Anfangswert schaltet, 
wobei wegen ; 


u(0) = uk(0) + > cu) = 0 + cur 
j-1 


zu setzen ist u®(0) = c;-ı fürk=0,...,n— 1. g(s) ist dann wegen g(s) = f(s) + u(s) leicht 
erhältlich?). 

Ist der Kern K(s, t) Summe von m derartigen Gr£EENschen Funktionen G,(s, f), so ist klar, 
wie man vorzugehen hat. 


Man setzt 
u = GI) gl) dt, ..., Um = S Gn(s; I) Kt) dt, 
dann ist : x 
ı)=f+ 2 U 
und für die u, u,..., Um gilt : 
/J,luı] = (I... ‚AmlUm] —g 
oder nach Einsetzen von g(s) 


ilel=f+ > re All TE 5 ER SR (3.58). 


Dazu kommen noch jeweils die zu jedem G gehörigen Nebenbedingungen der Form (3.4). 

Nun bestimmt man solche Lösungen von (3.5a), durch die sich alle Lösungen von (3.5a) 
als Linearkombination mit unbestimmten Koeffizienten c,, darstellen lassen, setzt in die Neben- 
bedingungen ein und erhält damit ein lineares Gleichungssystem zur Bestimmung der c,,. Durch 
sie sind dann wie oben die Anfangswerte von (3.5a) festgelegt, für die sich die u, gerade so er- 
geben, daß n 
s=-I+ > u, 

u=1 


wird. 


s 3. VoLTERRASche Integralgleichungen 

e Hier ist B=8,+%, mit , =fs>t 0<st<a} und ,=1s<t; 0<s,t=sa). 
2 Der Kern K(s, f) ist gegeben durch 

„ 

> 

$: 


FR für deL, 


K(s, f) = i,y=1 

24 0 für (,HdEB,. 
u 
& Aus dem Differentialgleichungssystem (2.7) wird hier 
> FOWEHSNEHFISEENHN! en dus et. aut (3.7). 
f (Die Vektoren haben n, Komponenten.) 
“ Die Bedingungen (2.4c) entfallen, weil keine Rechteckbereiche vorhanden sind, und aus 
r (2.8) wird wegen s,ı = 0 und dem identischen Verschwinden der Aıın, Rain; -- 
en TE BR ee (3.8). 
; (3.7) und (3.8) zusammen stellen aber ein Anfangswertproblem für ein System gewöhnlicher 
\ Differentialgleichungen dar, das ohne zusätzliche Rechengänge direkt mit dem Analogrechner 

gelöst werden kann. 
£ 2) Integralgleichungen mit Kernen der hier betrachteten Bauart werden von P. Mössıc in [10] und 


[11] untersucht, 


- 
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9(s) ist damit sofort aus der Rechenschaltung zu entnehmen. Der Spezialfall der VOLTERRA- 
schen Integralgleichung mit obigem Kern ist also der Lösung mit einer elektronischen Analogie- 


rechenanlage am besten zugänglich. 
Auch hier sind — zwar keine prinzipiellen — aber schaltungsmäßige Vereinfachungen 


möglich, wenn K(s, t) in speziellerer Form auftritt. : 
Zunächst nehmen wir an, daß K(s, t) in der Form K(s, f) = 266, ü) vorliegt, wobei 


@ ”- 
die G,(s, I) sogenannte einseitige GREENSche Funktionen von linearen Differentialausdrücken Ay 
der Ordnung m, sind. D. h. G,(s, f) soll als Funktion von s aufgefaßt der Differentialgleichung 


4,16] = 0 
genügen mit den Anfangswerten 
0 Bee Fa Ka 
mu — == 1 0 zZ ————— — 0 x = — = 1 
& os G A? 5 Fa 


fürs=tr+0;t=r; (<Str<a) Wir setzen dabei noch voraus, daß die höchste Ableitung 
in A, den Faktor I hat. Setzen wir hier 


hl) = [EN gOd 2 


dann genügt h,(s) der Differentialgleichung 


A| = HET, = we ne en Re (3.10) 
mit den Anfangswerten 


h0) =, = EN 0) =D. 
Nach Einsetzen von (3.9) heißt unsere Integralgleichung (1.3) 


«= + 2 KÄ3 =. an: a0 7 we 
Tragen wir dies in (3.10) ein, so haben wir 
Alk] =) + 2 he); al... 00 a BE 
mit den Anfangsbedingungen 
h,(0) = h5(0) = «= heT)0) = 0 für e=1,...,r. 


Schaltet man statt (3.7), (3.12), so kann man im allgemeinen mit einer Einsparung an Rechen- 
elementen, speziell an Multiplikatoren rechnen. 


Wenn man nämlich wieder von den Rechenelementen absieht, die zur Erzeugung der 
k,„ 1, bzw. der Koeffizienten von (3.12) nötig sind, braucht man für die 


n=% m, Produktein 9g=/f+ 3 k,h, 
e=1 v=1 


n, Multiplikatoren. Weitere n, braucht man, um , =1,g®=1,...,n,) zu bilden: Insgesamt 
also 2 n, Multiplikatoren. 


In (8.12) dagegen steht in jeder der r Zeilen eine lineare Differentialgleichung m,-ter 
pP #* . * * ” r 

Ordnung, für deren Schaltung man m, Multiplikatoren benötigt, das macht insgesamt ) m, =n, 
Multiplikatoren. te 

(Diese Zahlen stellen den jeweils ungünstigsten Wert dar, weil ja unter den Koeffizienten 
von (3.7) bzw. (3.12) sich auch Konstanten befinden können.) 

Ein weiterer Spezialfall ergibt sich aus Folgendem: (siehe [12]). 

Bezeichnen wir mit 9,(s) das u-fach iterierte Integral von 0 bis s über 9(s), so ist bekanntlich 


8 


A nf. DM | ? 
9,45) = FRIN (OD dt; Bel... Dee 


0 
Hat der Kern nun die Gestalt 
(s — Hr 


Ks; = als) ——— a ET 
2 (5) ... 8.14), 
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so wird unter Verwendung von (3.13) aus der Integralgleichung: 


gs) = Ns) + z Au+1(8) / = — gl) dt, 
0 


ELBE (3.15) 
06) = 16) + 2.06) 946) 
Führen wir noch die Hilfsfunktion &(s) ein mit 
ys) = Im(S) ; 
so wird aus (3.14) die Differentialgleichung 
B — 2" ne Ne (3.16) 
mit den aus (3.13) folgenden Anfangsbedingungen 
| 0) =p0)=-..-=- gm MN) =-0 ... 2.22... (3.16a). 


Die gesuchte Lösung g(s) der Integralgleichung ergibt sich in der Schaltung von (3.16) 
sofort aus 9”%s) = g($). 

Dieser Sonderfall gewinnt Bedeutung, wenn man einen allgemeinen Kern durch einen 
schaltbaren approximieren will. Dann liegt der approximierende Kern etwa in der Gestalt 


R6,) = X ü,6) 1" a Ne (3.17) 
es 


vor. 
Die a,(s) aus (3.14) erhalten wir dann durch: 


RE 
= Zu) a Ye) 


v 
2 gi 
= 26-93 Ya): 
Hieraus folgt durch Vergleich mit (3.14) 
a,(s) = De: yr = Tan ri (VE TE (3.18). 


Liegt also der Kern in der Eos (3.17) vor, so erhält er die Gestalt (3.14), wenn man die 
a,(s) gemäß (3.18) wählt. 

Damit hat man wieder die Möglichkeit, zur Lösung der vorgelegten Integralgleichung die 
Schaltung gemäß (3.16), (3.16a) vorzunehmen. 


$ 4. Abschätzungen des Verfahrensfehlers bei Approximation von allgemeinen Kernen durch 
solche des schaltbaren Typs 


1. Werden Integralgleichungen vom schaltbaren Typ mit Hilfe einer elektronischen Ana- 
logierechenanlage gelöst, so treten nur die durch die prinzipielle Ungenauigkeit des Analogie- 
verfahrens herrührenden Fehler auf. Die Auswirkungen dieser Ungenauigkeiten auf die Lösung 
von mathematischen Problemen zu untersuchen, ist nur bei sehr einfachen Aufgaben gelungen 
und auch hier nur in grober Näherung. Deshalb werden diese „‚apparativen Fehler“ nicht näher 
betrachtet. 

Wie aus den Sätzen 1 und 2 hervorgeht, ist der Typ der schaltbaren FREpHoLMschen 
Integralgleichungen zweiter Art eindeutig festgelegt. Man kann also prinzipiell nicht erwarten, 
mit Hilfe des elektronischen Analogrechners Lösungen von linearen Integralgleichungen zweiter 
Art mit weniger speziellem Kern kontinuierlich als Funktion der Maschinenzeit erzeugen zu 
können. 

2. Liegt nun aber eine nicht schaltbare Integralgleichung zweiter Art vor, so bleibt, um 
zu ihrer Lösung die elektronische Analogierechenanlage einsetzen zu können, nur der Ausweg, 
sie durch eine schaltbare Integralgleichung zu approximieren. 

Wir wollen hier auf jeden Fall voraussetzen, daß (5) stets stückweise glatt ist, so daß 
hier zur Erzielung der Schaltbarkeit kein Ersatz durch eine andere Funktion nötig ist. 
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Es ist also nur der gegebene Kern Ks, {) durch einen vom schaltbaren Typ K(s, t) zu 
ersetzen. 


Damit für die hierbei nötigen Rechnungen auch der Analogrechner eingesetzt werden kann, 


werden wir nur Kerne K(s, t) aus T, (siehe $ 1.1) betrachten. " R 
Nennen wir die Lösung der vorgelegten Integralgleichung (mit dem Kern K(s, )) 96) 
und die der approximierenden Integralgleichung (mit dem Kern K(s, f)) g(s), so gehören (siehe 
$1.1) 9(s) und g(s) zu I. } 
Als Maß für den durch die Approximation bedingten Verfahrensfehler führen wir 


I —IS)|- -- - - et ee ee 


ein. (Dies ist sinnvoll, da 9, ge I, und deshalb das Fehlermaß immer endlich ist. Außerdem 
ist dieser Fehler in der Praxis am einfachsten deutbar; weiter läßt sich daraus der Quadrat- 
mittel-Fehler abschätzen.) 

Damit die Abschätzung von |5— g| günstig erfolgen kann, werden wir als Maß für die 


Abweichung zweier Kerne K und K in einem Punkt (s,, {,) wählen 
| K(s0, 1) — K(s0; bo) | - 


Daraus ergibt sich, daß es dem Problem am besten angepaßt ist, im Sinne von TSCHEBYSCHEFF 
zu approximieren. 


Besonders praktikabel wird der Einsatz einer elektronischen Analogierechenanlage zur 
Ermittlung des approximierenden Kerns, wenn der gegebene Kern K(s, f{) von der Form 


Ks, = F(P(k,6) ; 10) 


ist, wobei P= I k,(s) 1,() ist. 
Man identifiziert hier zunächst P mit der Maschinenzeit (bzw. mit einer dazu propor- 
= N 
tionalen Größe) und hat nun die Aufgabe in 3 c, P® die c, so zu bestimmen, daß F(P) im Sinne 


o-=1 
von TSCHEBYSCHEFF durch 3 c, P® approximiert wird. Diese wird gelöst, indem die Differenz 
r =1 


ES e= 
D(P) = F(P)— 3% c, P® schaltet und durch systematisches Probieren für die c, die passenden 
=1 


e 
Werte ermittelt. 

Als Richtschnur für das Probieren dient das von Hasrısss in [13] verwendete Verfahren 
in sinngemäßer Übertragung. 

Ist F(P) stetig, so hat man ein mit dem elektronischen Analogrechner leicht kontrollier- 
bares einfaches Kriterium für die Bestimmung der c, (siehe z. B. [14]): 

D muß an mindestens r + 2 Stellen seinen Maximalbetrag und zwar mit abwechselndem 
Vorzeichen annehmen. 

Hierauf kann man durch Änderung der C„ das einfach durch Drehen von Potentiometer- 
knöpfen erfolgt, bequem hinarbeiten. Setzt man dann P = N k,(s) I,(t) in NY c, P®? ein, so hat 
man nach Umordnung einen Kern des schaltbaren Typs erzeugt, der den gegebenen möglichst 
gut annähert. 

Im allgemeinen Fall wird man auch darauf achten, den Aufwand an Rechenelementen durch 
einen passenden Ansatz für K(s, {) möglichst klein zu halten. Dazu wird man zweckmäßig die 
in $3 beschriebenen Fälle heranziehen. 


3. Es bleibt nun noch übrig, Schranken für |9 — g|, zu deren Berechnung man möglichst 
auch wieder die elektronische Analogierechenanlage heranziehen können soll, anzugeben. Wir 
beweisen in diesem Abschnitt zunächst folgenden 

Hilfssatz: 

In der Integralgleichung für A}: 


hahrt lea) va, 
seien fi e/\ und L,E I, Weiter sei ‚(s) Ze>0in0<s<aund LS) ZzZ0in0<sst<sa. 
(4.3) habe ein eindeutige nichtnegative beschränkte Lösung h,(s) el. 


Sei ferner /(s) Z |/s(s)| in [0, a] und L,(s, ) > |Lxs, d| in ((<s,t< a), dann gilt für jede 
Lösung h,(s) von 


L ef + el a 
h(s) = |hs(s)| in [0, a]. 


BE EL A ZREB 2 52 


rk 


TERTEFATTLEN 
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Beweis: 
Wir zeigen zunächst, daß der lösende Kern H,(s, !) von (4.3) nicht negativ ist. Dazu führen 
wir die Funktionen 


n—1 
Su na en En ae (4.5) 


ein. (Zur Bezeichnung siehe $ 1.1). 


Ausf, 20 undL, = 0 folgt, daß für festes n jeder Summand der die Funktion g1„(s) defi- 
nierenden Summe nichtnegativ ist. 


Für h‚(s) gilt nach Voraussetzung 
Kein oder Sf = ee (4.6). 


Ersetzen wir /, in (4.5) durch diesen Ausdruck, so erhalten wir 


AR Z n—1 " he: BR n—1l ER ” E 
9 =h—Uh +3 (&h, US h)<h+ ih N &h 


in — h, = uh 1 
Wegen L, > 0 und hy > O0 erhalten wir hieraus 
hı(s) > 9ın6) für Der und re iny 
Lassen wir nun in (4.5) n— oo gehen, so durchläuft (für festes s) 91, eine monoton nicht 


fallende Folge von Werten, da alle Summanden von gı„ positiv sind. Weiter ist nach (4.7), da 


h, als beschränkt vorausgesetzt wurde, die Folge gı„(s) beschränkt, d.h. gı„(s) konvergiert in 
0 <Ss<S a fast überall. 


Wir setzen 


h,6) = lim 91.6) = lim (+ 24) a (4.72). 


NnN>XO 
Aus der Definition von gı„ folgt weiter 


Haar hhn Baal) 2 
Daraus erhalten wir für n — © 


h, = fı + lim 8, 9ın- 


Nn> X 


Da gı„ eine konvergente monoton nicht fallende Folge ist, darf der Limes mit %, vertauscht 
werden (siehe z. B. [9], S. 25{f.) und es gilt 


Behneuh 


h, ist also Lösung von (4.3). Da nach Voraussetzung h, die einzige Lösung von (4.3) sein sollte, 
ist A,(s) = h,(s). Setzen wir dies in (4.7a) ein, so haben wir 


no 00 


h(s) = = Et R gınls) = — lim (ht, 2 >> & h): | 


he) = = fi + lim zZ % fı | 


n>xov=-=1 


Aus den Betrachtungen, die zu (4.8) führten, folgt, daß auch die Lösung Ah} der Integral- 


gleichung 
h=Cch +44 (C = tonst) 
in der Form (4.8) als 


n—1 
Di Ch Blume DC te 9) 
n>@v=1 
dargestellt werden kann. 
Nun genügt der lösende Kern H,(s, t) von (4.3) der Integralgleichung 
ee (4.10). 


Sie ist eindeutig lösbar, weil andernfalls auch (4.3) entgegen der Voraussetzung nicht ein- 
deutig lösbar e Wegen der Beschränktheit von L, und fı Ze existiert eine Konstante C, 
sodaß Ch 2 
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Auf die Folge 
n—1 
Ln=L+2&h 1 
können nun die gleichen Überlegungen angewendet werden, wie auf die Folge gı„, denn Lı„ ist 


wegen L, > 0 monoton nicht fallend und L,„ wird durch hf majorisiert. 


D.h. esist 
=L,+lim > Suhl, 
n—>X0 ı= 
DaL, > 0 sein sollte, gilt 
Hs) =0, Ben 


womit der erste Teil unseres Beweises abgeschlossen ist. 


Aus (4.4) folgt nun 
Sl Hall - 22:4 2.8: u 
wobei |%,| den zu IL,(s, f)| gehörenden Operator kennzeichnen "soll. 
Subtrahieren wir nun (4.11)vonhy=fi + g h, 
h— |) 2 h— Ihl + 9, — [£2| Ihe - 
Zur rechten Seite u wir — Qlha| + X |ha! und erhalten 
4 — |hel 2h— Il + Sl — IRz]) + (& — IE) Ihe] - 


Dies schreiben wir in Form einer Gleichung, indem wir rechts eine Funktion h(s) addieren, von 
der wir nur wissen, daß h(s) = 0 ist: 


h— alu) + h— id + a — (A) + as — Id il. 


Dies ist eine Integralgleichung für (h, — |h,|) mit dem Kern L,(s, {), deren Lösung wir mit Hilfe 
des lösenden Kerns H,(s, f) schreiben können als 


h, Ye [he] =h+(h— rl) le IXe]) Ih, + 9ı h+ 9ı (= I!) =k 9ı er K)) [ha| (4.12). 


Nach unseren Voraussetzungen ist nun i — al 20, L,— |L,| = 0; weiter wissen wir, daß 
h=0. Da außerdem bewiesen wurde, daß H,(s,d) = 0, Pe aus der letzten Gleichung 


a [hl > 0 oder = |R!, 
was zu beweisen war. 
4. Wir geben nun folgende 
1. Fehlerschranke: 
9(s) sei eindeutige Lösung der Integralgleichung g= +89; feT}; XeT,), %s) Lösung 
vn9=f/+89; (ers Ri e I,). Weiter sei 
\K6s, D—Kis, n) sr= const für V(zEs.s> m, 


L(s, t) sei ein schaltbarer Kern mit L(s, ) > IK6s, ) und die Integralgleichung 
a a 
h= Kf Id] dt 2 Us, DD dns Mrs (4.13) 


habe in0O=s=< a die eindeutige positive beschränkte Lösung F(s). 
Dann ist 
Is) — gs)| Z F(s) in (En: 

Beweis: 
Wir setzzeng=g+q,undK=K+H X. 
Dann folgt ug = /+ 89 

ITtn=-/+tR+D)gH+gNM- 
Subtrahieren wirg=f+ 89, so bleibt 
erst tn ha tn... er (4.14). 

Wir schätzen nun ab 


ol SS 1K6,daMar<r | \aoldt. 


x 
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Jetzt wenden wir den obigen Hilfssatz an, indem wir identifizieren 
fi 
Heute { \s(i)| dt ; Fe ee 12: ge nurar. En, 


In eg R,- mit >K: ara. 
Man überzeugt sich leicht, daß die rechts stehenden Funktionen wirklich die im Hilfssatz 
geforderten Voraussetzungen erfüllen. 
Aus der Rückübersetzung der Aussage des Hilfsatzes 


hı > |h,| 


folgt nun 
F2|\nl=9—9]; 


also unsere 1. Fehlerschranke. 
Bemerkungen: 


1. F(s) ist mit der elektrischen Analogierechenanlage kontinuierlich erzeugbar! 
2. Aus dem Beweis in 3. folgt: 


Notwendig für die Existenz einer solchen Fehlerschranke F(s) ist, daß 9 = + g durch 
Iteration lösbar ist. 


3. Damit die Schranke scharf ist, muß nach der Übersetzungstabelle und (4.12) gelten 
>20; Ku=r=0 _ (dies entspricht f, — |fs| > 0) 


ind.) = 14 (was sich aus L, — |L,| = 0 ergibt). 

Die obigen Bedingungen sind notwendig und hinreichend dafür, daß F=g, (in (4.12): 
hh— !h| =. 

Es gibt nun sicher ein K und ein /, so daß die obigen Bedingungen erfüllt sind: 


Sei z.B. K6s, DESUVEr > 0 so klein, daß K = K-r >0undf>0. DaK durch L(=K) 
majorisiert wird und (4.13) nach dem Hilfssatz durch Iteration lösbar ist, ist auch g=/+g 
durch Iteration lösbar, und hat, da /, K >0, als Lösung eine nichtnegative Funktion g(s). 


(Praktisch wird dieser Fall nicht eintreten, da L schaltbar sein sollte und wegen L=K 
auch K schaltbar wäre. Man ersieht jedoch aus Obigem, wann F(s) günstige Schranken liefern 
wird.) 


5. Ein Beispiel: 
g(s) = f(s) + 10:25 log(1 + st) g(l) dt. 


Der Kern K(s, t) = 0,25 log (1 + s £) ist nicht schaltbar. Nach $ 4.1 approximieren wir ihn 
nach TSCHEBYSCHEFF durch ein Polynom in s i, wobei wir das IntervallO < s s 3 zugrundelegen 
wollen. 


3 
Mit einem Polynom K = N c,(s t)? erhält man durch das Probierverfahren für die c, 
oe=0 


0,0155 Er 2 e=1),.0303; & — 0,0161 
wobei [X — K) < 0,015 (=n). 
Nach $ 3.3, speziell (3.18), wird mit 9'%(s) = g(s) aus 
4 
ae ll ar 2 da) 
Ne ’ 

gWX0) = 0 A ER) 
mit 

as)=— tar tas®t+as), 

ESF OS FR T3G3, 

a,() = — (2,2+66%5%), 


243) = CS. 
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Die Schaltung von (4.13) liefert nun für jedes zur Verfügung stehende /(s) sofort g(s). 
Wählt man z. B. /(s) gemäß Bild 2, dann ist das zugehörige g(s) in Bild 3 dargestellt. 


Da K > 0 ist K| = K + r und wir können in (4.13) für L= K + r setzen. Außerdem ist 
K,=0fürt>s, so daß aus r [ Ig(t)| dt wird r E Is] dt. 
0 0 
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Dann kann mit kleinen Abänderungen der Schaltung von (4.15) (c, muß durch Cg + r ersetzt 
werden) sofort F(s) erzeugt werden, wenn statt /(s) die Störfunktion r f |g(d)| dt verwendet wird. 
() 


In unserem Beispiel ist F(s) in Bild 4 gegeben, und nach Bemerkung 3 im vorigen Abschnitt 
wird der wahre Verfahrensfehler nur wenig unterhalb F(s) liegen. 


6. Die 1. Fehlerabschätzung beruht im wesentlichen darauf, daß in gewissen Fällen der 
zum vorgelegten Problem gehörige lösende Kern dem Betrag nach majorisiert werden konnte 
durch den lösenden Kern einer schaltbaren Integralgleichung. 


Wir wollen hier auf anderem Wege eine andere Abschätzung geben, die in jedem Fall 
anwendbar ist. 


Es gibt zwar eine Abschätzung von TrıcomI (siehe [15], S. 76 und die dort angegebene 
Literaturstelle) jedoch ist dazu die Berechnung von FrepHuormschen Nennern D(A) nötig, die 
nicht einfach auszuführen ist. 


Hier soll mit Hilfe der Störungsrechnung aus dem als bekannt anzusehenden lösenden Kern 
H(s, it) vong=f-+ Sg der lösende Kern Es) vong=f+ E27 ermittelt werden. 
Es sei | 
K=K +eK,. 


Für den zuK gehörenden lösenden Kern H machen wir den Ansatz 


und gehen damit ein in 
H=-K+SH=(K+EeK)+RHteR)B: 
Seß-Kt4sK4RDVeC+RaTeHe. 
v=0 v=0 v=0 
Koeffizientenvergleich liefert: 
ER RG,, 
G, = K, + + RG 
= RG + RG W237.) 


Es bezeichne H den zu K gehörigen lösenden Kern. 
Dann erhalten wir 


GEH, 
= [+9 RPTAHNK + KH) al) 22 2a, 


wobei [Operator] ® gleich 1 sein soll. 
Setzen wir noch e = 1, so ergibt sich: 
Der zuK = K + K, gehörende lösende Kern H ist (zunächst formal) gegeben durch 


B=-H MH TUA+HSKRFRH ...:-.- (4.16), 
v=1 
wo H der zu K gehörige lösende Kern ist. 
Zur Anwendung dieser Formel setzen wir voraus, es sei 
IK] Sr = const 
bekannt. Außerdem sei K schaltbar, also H beschränkt. Wenn nur r so klein ist, daß 
(1 Hari: (A = sup f Hs, | di), konvergiert obige Reihe absolut und gleichmäßig und 
0ssSsa0 ” 
die rechte Seite von (4.16) ist wirklich gleich H. 


26* 
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Die Lösungen vng=f+ $g und g= + 8g heißen unter Verwendung von H und H F 
=i+5f. =ehrl 
woraus durch Subtraktion folgt 


9) — | = SH; 
(und mit (4.16)) 


ZA+SRTTA+NAAHD I, 


(dag=/+9N 


zu +9 FH 


PA ES LNZE 


Nach 82 ist H(s, {) als Funktion von  (l Maschinenzeit) erzeugbar, ebenso auch g. Damit 
ist mit dem Analogrechner berechenbar: 


a 


S sd] dt und 7 Hs, t)| dt. 
f) 
Es bezeichne wieder 


H=sup [ |H(s, D)| dt. 
0<s<a 0 


Dann ergibt sich unter Verwendung von |K)'<r 


ZIA+VAr SEA + DA + A tar fowla, 


wo |$| den zu |H| gehörigen Operator bedeutet. 
Weiter umgeformt: 


2 . 1+H)ra 
[9(8) — 965) Sr) ls@| di v +) IHKs, D)| a) k is ee .. . (417). 


Wir fassen zusammen: 


2. Fehlerabschätzung: R 
Es sei g(s) eindeutige Lösung von g= f+ Sg (K schaltbar) und 9 Lösung von 9=f+ 89. 
Ferner sei IK -- K| <r. Dann gilt, falls 1—ar(1 + H)>0, für |d — g| die Ungleichung 
a 
(4.17), wobei H den lösenden Kern vong— + Sg bezeichnet und H y up J |H(s, t)| dt ist. 
sssa 
Bemerkungen: 
1. Falls r genügend klein ist, ist die Fehlerabschätzung immer anwendbar. 
2. Da H nur mit vielen Rechenläufen erhältlich ist, wird man sich damit begnügen, sich 
einen Anhaltspunkt für den Fehler zu verschaffen, indem man Glieder höherer Ordnung in r 
vernachlässigt. (Dann braucht man H nicht mehr.) 
3. Die Schranke ist z. B. scharf für K = c = const. >0 mit 1—ac>0 K= Kr =: 
7 = const > 0 und /(s) > 0, wie man durch Einsetzen in (4.17) und Vergleich mit der bekannten 


Lösung sieht. 
$ 5. Integrodifferentialgleiehung 


Die Betrachtungen von $2 können ohne Schwierigkeit ausgedehnt werden auf die Lösung 
von Integrodifferentialgleichungen, wobei wir uns auf solche der Form 


Al=rFrHr RA]. 2 Se (5.1) 
beschränken wollen. Dabei sei 
),[g] = 9) — 2 a,,6) gl) 50 7 y = 1a SE (5.1a) 
A 


und esseir, =T, 


N we age 


BEE UNE 


we. 2 F- E . } > 
—.. - ö 
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Durch fast. ichs Wiederholung der Schlüsse von $1 ergeben sich als notwendi Be 
dingungen für die Schaltbarkeit von (5.1): E 5 a 
r K muß den gleichen Bedingungen wie in Satz 1 genügen und die Un und f müssen stückweise 
glatt sein. ; 

Um die Lösung von (5.1) festzulegen, sind noch Anfangs- bzw. Randbedingungen für g vor- 
zugeben, die wir abkürzend 

FEr BIN EEE EEE (3.2) 

nennen wollen, wobei R linear in g sei. 


Da das Lösungsverfahren von (5.1) keinen wesentlichen Unterschied vom Verfahren zur 


.. . - 2 
Lösung der Integralgleichung (1.3) aufweist, wollen wir hier nicht auf Einzelheiten eingehen, E 
sondern nur das Prinzip des Vorgehens schildern. r 

Für die Rechteckbereiche (siehe $ 2.1) setzen wir: . ; 
j log) v- i 
e (Es = S LO Alg] dt; a a ee er en (5.38), & 
“ ugple) 


für die Dreiecksbereiche mit Scheitel ins >t: 


Bed eo TR (5.3b), 


bug) 


für die Dreiecksbereiche mit Scheitel ins < Et: 


lo 2) 
= — [ 10%lg] dt; lo) 1 me TR (5.3e). 


Aus den letzten beiden Formeln gewinnen wir wieder ein Differentialgleichungssystem für 
die h,, das hier durch (5.1) ergänzt wird zu einem System für die h, und g: 


h,(s) = 1,(8) Aglg] ; (6) =1,...,m), 
MR m, Mm; Mg a. (0.4)% 
E Al=I+ SKI) — F Kdik)+ 3 ck) = 
b o)=1 YW)=1+m, YW)=1-+m, 


Dies ist ein System gewöhnlicher linearer Differentialgleichungen für die h, und g, das noch 
unbestimmte Parameter c, enthält. Diese und die noch fehlenden Anfangsbedingungen für 
die h, und g verschafft man sich genau wie in $ 2.2. (Man erkennt, daß der einzige Unterschied 
zu dem in $ 2 gelösten Differentialgleichungssystem darin besteht, daß Ableitungen höherer 
Ordnung auftreten und man könnte durch Einführung neuer Funktionen u,(s) = gs) (u = 0, 

., 71 — 1) sofort formale Gleichheit erzielen). 

Wir erzeugen zunächst ein vollständiges Lösungssystem des zu (5.4) gehörenden homogenen 
Systems auf die in $2 geschilderte Weise. Danach stellen wir mit den Anfangsbedingungen 
h, g) = 0 an der Stelle s = O0 nacheinander für die Störfunktionen f(s); k,(s) (aus den Rechteck- 
bereichen) je eine spezielle Lösung her. 

Aus allen diesen (mit Hilfe des Analogrechners wohlbekannten Lösungen) wird durch Linear- 
kombination mit noch unbekannten Koeffizienten, (unter denen sich auch die von (5.3a) befinden) 
die allgemeine Lösung von (5.4) erzeugt. 

Diese setzen wir in (5.2), (5.3a) und (5.3b, c) hier für s = tugu) DZw. Ss = Iogg)) ein. 

Dies Einsetzen ist möglich, weil uns durch den Analogrechner die Funktionen des voll- 
ständigen Lösungssystems bekannt sind. 

Als Ergebnis erhalten wir ein lineares algebraisches Gleichungssystem für die Koeffizienten 
der Linearkombination. Sind diese berechnet, so sind damit solche Anfangswerte von (5.4) be- 
stimmt, für die die Lösung von (5.4) die Nebenbedingungen erfüllt. In der Schaltung wird dann 
die gesuchte Funktion g(s) erzeugt: 

(5.1) ist gelöst. (Ob (5.1) überhaupt lösbar ist, ersieht man aus dem linearen Gleichungs- 
system für die Koeffizienten der Linearkombination der allgemeinen Lösung von (5.4).) 

Erwähnt sei, daß die Lösung von (5.1) ohne zusätzliche Auflösung eines linearen Gleichungs- 
systems mit der elektronischen Analogierechenanlage erzeugt werden kann, wenn K ein VOLTERRA- 
kern ist, und R[g] = 0 Anfangsbedingungen für s = 0 sind. 


Ka Dal a u Zu A 05 oc un nd iu 
u. 1 K e 
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KLEINE MITTEILUNGEN 


Flexure of a Compound Cylindrical Beam From the displacement — strain relations we 
Composed of Non-Homogeneous Material obtain 
Introduction u = 1: 0080 — NE 0050 Bu, ;\ 

That non-homogeneity has gained some impor- 
tance in elasticity and plastiecity is evident from u — EuW + En 0 (3) 
various types of papers read at the International 2K 6K “2 ; 
Symposium on Non-Homogeneity held at Warsaw [1]. W2r 

In this paper formulation of the flexure problem be le cos8 + Wffr)cos0 . . . (4), 
of a heterogeneous ceylindrical beam is given. PoIsson’s ö 
ratio (denoted by 0) is supposed to be constant E(r) W cos® [,, or 1 
throughout the beam and Young’s modulus is taken I Ah ee (l+o) f r) a8 + ®) | 
to be a function of the distance from the axis of the i ä 
beam. Numerical results for a particular case of it gi Eir) W sin [f(r) or: 
are also evaluated. "7 2/l+- 0! | vr art (6), 


where f{r) is a function of r to be determined. 


j Substituting the expressions for stresses from (1), 
We take Z axis along the axis of the eylindrical (5) and (6) in the equilibrium equation 
beam. Let the beam be subject to a force W in the 
plane z = 0 along the radius d=0. Letr = r, and OTer | 1 8705 , Ts | Trs 
r=r(r,>r,>0) denote the boundaries of the MEZ ae Pe Ps IE 
eylinder. Let r=r,(n>r,>r,) separate each x 
cross section into two parts R, and R,. Let R, and We obtain following differential equation 
R, denote the portions of the cylinder for which ! 
nsrsr) and (nSrse,) respectively and Ef+Ef'+ ar _ ar Fur 2Er 

r 


E,(r) and Esfr) be Young’s moduli in regions R, and Rr SEEN e 
R, respectively. 


Solution ofthe problem 


Now in this case obviously no torsion takes place nn N a ER 
due to symmetry. : Let us consider a particular case where E = E, m; 
E, is a constant and n is a real number. 


We assume as usual : 
{ Then equation (8) reduces to 


Wzr 


Tr =W=mw—=), Tp = E(r) cos® (1), 
K Bat dr N 
where K is constant to be found out. 2 
e Hence 
Then e, = EN co80, Hr=eH = — im rcos0 ) = An VRR an 
K K Nena rtirrom, 


We also assume T77,, T9; and hence &r, €, are inde- 


pendent of z. where A, B are constants and p,g are roots of the 


equations m + nm —1=0 


# 


er 
s 


g' 


en 


ENT I ER 
eg Bee 
Let fı(r) and f,(r) be values of f{r) in R, and R, 

respectively. . i 
Bess NT on the free boundaries 77, = 0; hence f;(r) 
Er) 


onr=r, and r,. 


2K 
Displacements and stresses are continuous on the 
_  eommon boundary r = r,. 

Hencee _ 


Bin) In ie = 3,(n) It | 


andfı(r) = hfrı)- 
The unknown constants involved in f,(r) and 
fs(r) are given by the above boundary conditions. 
_ Let \ 
De T,=T%zsinO + 79,c0os0 
and 


> ‘3 


T, = Tr, c0s0 — 79, sin. 
Along the horizontal diameter (Od = + 2) Te=0; 


‚hence T, = 0. Let us consider a particular case for 
Br which r, =0. 
E, = Youne’s modulus in R, is constant and 
2,—=Hrnin R(n sr =n,). 
Let 7 be the value 7, at the centre of the circle 
r=0 where it is maximum and let 7’ denote T, 
at the ends of the diameter of the cirele (r = r,). 


Buchbesprechungen 


Pr 


Following table gives numerical results for T 
and 7’. Young’s modulus is supposed to be con- 
tinuous on the common boundary i.e. BE) = Eyr" 
and o is taken to be equal to 0.3. Ei 


Table 
en 2 3 
n —1| 0 1|2|-—1 o Va? 
Gr r3 |1.71|1.38|1.08|0.81 2.04 1.38 0.86 | 0.49 
si | 
wer 0.82 |1.23 |1.79 |2.52 |0.75 11.23 11.91 | 2.80 


I am indebted to Prof. B. Sen D. Sc. of Jadavpur 
University for kindly suggesting the problem to me. 
My thanks are also due to Prof. C. WEBER, Hannover, 
for showing me the simpler method of solving it. 
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A. Linder, Statistische Methoden für Natur- 
wissenschaftler, Mediziner und Ingenieure. Dritte 
Auflage. 4848. Basel und Stuttgart 1960. Birk- 
häuser Verlag. Preis geb. 58,60 DM. 


Gegenstand des vorliegenden Buches ist die Aus- 
wertung von Beobachtungsergebnissen nach statisti- 
schen Gesichtspunkten. Die Vorzüge des Buches sind 
bereits nach der 2. Auflage 1951 von FR. A. WILLERS 
gebührend hervorgehoben worden (ZAMM 32 (1952)) 
Inzwischen haben die statistischen Verfahren in den 
verschiedensten Gebieten der Praxis Eingang gefun- 
den. Mit der vorliegenden stark erweiterten und in 
wesentlichen Teilen umgearbeiteten dritten Auflage 
(mit der doppelten Seitenanzahl gegenüber der 2. Auf- 

’ lage) kommt Verf. dieser Entwicklung entgegegen. 
{ Allerdings weist das Verzeichnis der 73 Beispiele im 
überwiegenden Teil Anwendungen aus Medizin und 
Biologie auf. Die Neuaufteilung und Erweiterung des 
Inhalts erkennt man am besten durch die Kapitel- 
überschriften: Häufigkeitsverteilung, Durchschnitt 
und Streuung; Schätzungs- und Prüfverfahren; Be- 
urteilen von Häufigkeiten; Beurteilen von Durch- 
schnitten und Streuungen; die Streuungszerlegung; 
Abhängigkeiten zwischen meßbaren Merkmalen; 
Schätzen von Parametern; Numerisches Rechnen; 
Theoretische Grundlagen. 

Durch die Aufnahme moderner Schätzverfahren 
hat das Buch für den Praktiker noch an Wert gewon- 
nen. (Es sei erwähnt, daß das Wort ‚schätzen‘ in 
der 2. Auflage nicht einmal vorkam). In diesem Zu- 
sammenhang wurde der Anhang um vier Tafeln für 
die Transformation von Prozentzahlen erweitert. 


Allerdings sind die Ausführungen im Abschnitt 
„‚Vertrauensgrenzen‘‘ irreführend. Entgegen der üb- 
lichen Terminologie werden die Vertrauensgrenzen 
eines unbekannten Parameters u wie bei R.A. 
FISHER im Sinne der ‚‚fiducial limits‘‘ eingeführt und 
zwar um ‚„‚die Unsicherheit einer Punktschätzung des 
Parameters u auch bei kleinen Stichproben zu be- 
f messen‘. So wird auf $8.55 schließlich behauptet, 


a a Eh in a in En nn una 


az ah 


ii a4 
. 


W 


man könne ‚eine Wahrscheinlichkeitsverteilung des 
unbekannten Durchschnitts u der Grundgesamtheit 
herleiten‘‘, obwohl « in jedem Fall zwar unbekannt, 
jedoch konstant ist. Daß es sich bei der Festlegung 
eines Systems von Vertrauensintervallen um eine 
andere Form der Schätzung von u handelt, wird 
hier ganz verschwiegen. 

Wenn die hier benutzte Punktschätzung von u 
eine diskrete Verteilung besitzt, werden ‚‚Mutungs- 
grenzen‘‘ angegeben (in der 2. Aufl. als ‚‚Vertrauens- 
grenzen‘‘ bezeichnet). Daß es sich hierbei um die 
zweiseitigen Konfidenzintervalle nach CLoPPER und 
PEARSoN zum Konfidenzkoeffizienten ß = 0,9 handelt, 
wird sowohl in der Definition als auch in der Deutung 
nicht klargestellt. 

Dem Referenten fällt es schwer einzuschätzen, 
welche Folgen diese Art der Darstellung in der Praxis 
haben kann. Mit diesen Einschränkungen kann das 
Buch in der vorliegenden vorzüglichen Ausstattung 
sowohl dem praktischen Statistiker als auch dem 
Mathematiker empfohlen werden. 


Dresden W. RICHTER 


Archive for History of Exact Sciences. 
Vol. 1, Number 1. Edited by C. Truesdell. 106 S. 
Berlin/Göttingen/Heidelberg 1960. Springer-Verlag. 
Preis brosch. DM 19,60. 


Das der Schriftleitung vorliegende erste Heft dieser 
neuen Zeitschrift enthält einen Aufsatz des Heraus- 
gebers C. TRUESDELL ‚‚A Program toward Redisco- 
vering the Rational Mechanics of the Age of Reason“ 
und einen zweiten Aufsatz von ÄprAD Szas6 „‚Anfänge 
des euklidischen Axiomensystems“‘. Beide Aufsätze 
sind von hoher Warte geschrieben und lassen bereits 
die gediegene wissenschaftliche Arbeit erkennen, die 
hinter den Veröffentlichungen dieser Zeitschrift steht 
und für die allein schon der Name des Herausgebers 
bürgt. Das neue Unternehmen, durch das eine wesent- 
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liche Lücke in unserem Zeitschriftenbestand ausge- 
füllt zu werden verspricht, wird in weiten Kreisen sehr 
begrüßt werden. 


Dresden H. HEINRICH 


L.D. Landau and E.M. Lifshitz Mechanics. 
Vol. 1 of: Course of Theoretical Physics. 165 S. m. 
55 Abb. Oxford/London/Paris 1960. Pergamon Press 
Ltd. Preis geb. 40 s.net. 


Das Buch enthält eine Deduktion der theoretischen 
Mechanik. Ausgehend vom Hamıtronschen Prinzip 
der kleinsten Wirkung und vom Gariteischen Relati- 
vitätsprinzip wird die Theorie bis in ihre Details ent- 
wickelt. In sieben Kapiteln werden dargestellt: Be- 
wegungsgleichungen, Erhaltungssätze, Integration 
der Bewegungsgleichungen, Stoßprobleme, kleine 
Schwingungen, Bewegung starrer Körper und kano- 
nische Gleichungen. Dieser Stoff, der besonders die 
Fragen der ‚‚reinen‘‘ Mechanik betont, ist vorwiegend 
auf die Bedürfnisse des Physikers zugeschnitten. 

Dies ist kein Lehrbuch der Mechanik. Wer diese 
Disziplin erlernt, hat normalerweise zwar schon etwas 
über Vektorfunktionen, Gradienten, Differential- 
gleichungen, Variationsrechnung, Koordinatentrans- 
formationen und Invarianzeigenschaften gehört, ist 
mit diesen Dingen jedoch noch nicht sehr vertraut. All 
dies aber wird ihm hier bereits auf den ersten vier 
Seiten ohne große Vorbereitung abgefordert. Auch 
wird ihm unklar bleiben, welche speziellen Erfahrun- 
gen jeweils zu den einzelnen Begriffsbildungen geführt 
haben. 

Wer aber dieses Gebiet bereits kennt oder zumindest 
hierfür beträchtliche Voraussetzungen mitbringt, wird 
das nicht so umfängliche Werk mit allergrößtem 
Interesse lesen. Sein deduktiver Aufbau bietet in 
großartiger Weise einen Gesamtüberblick über den 
Stoff und zeigt die großen logischen Zusammenhänge 
der Theorie auf. Fast jeder einzelne Satz ist von be- 
sonderer Tiefe und Aussagekraft. Wer die Autoren 
persönlich kennt, spürt auch durch den englischen 
Text hindurch ihre besondere Art, physikalische Ge- 
danken aufzugreifen und zu entwickeln — ein Hin- 
weis darauf, daß die vorliegende Übertragung ins 
Englische als sehr gelungen bezeichnet werden muß. 


Dresden W. MAcKE 


W. Franz (o. Prof. a. d. Univ. Frankfurt), Topolo- 
gie. Band I: Allgemeine Topologie. (Sammlung 
Göschen, Band 1181.) 144 S. m. 9 Abb. Berlin 1960. 
Walter de Gruyter & Co. Preis brosch. DM 3,60. 


Das in der Sammlung Göschen schon lange erwar- 
tete Bändchen über Analytische Topologie ist nun 
erschienen. Ein weiterer Band über Algebraische To- 
pologie ist in Vorbereitung. Obwohl sich naturgemäß 
die Grenzen zwischen den beiden Gebieten der Topo- 
logie nicht völlig präzise angeben lassen, erscheint dem 
Referenten das Vorhandensein des IV. Teiles (Simpli- 
zialkomplexe, Polyeder, Pflasterdimension, nulldimen- 
sionale Kompakten, ALEXANDER-SPERNERSches Lem- 
ma, Pflastersatz, MENGER-NöBkrLınGscher Einbet- 
tungssatz) im vorliegenden Band I über analytische 
Topologie nicht ganz gerechtfertigt. An Stelle dieses 
Abschnittes wären vielleicht einige Ausführungen über 
Filter und uniforme Strukturen, von denen der Ler- 
nende in der gegenwärtigen mathematischen Situation 
möglichst bald die wichtigsten Eigenschaften erfahren 
sollte, zweckmäßig gewesen. Der einführende Charak- 
ter des Buches wäre damit durchaus nicht in Frase 
gestellt. Abgesehen von dieser Ansicht des Referenten 
muß man aber sagen, daß die ersten drei Teile eine 
außerordentlich gut lesbare und didaktisch hervor- 
ragende Einführung in die allgemeine Topologie dar- 
stellen, die insbesondere den Studierenden der Mathe- 
matik sehr ansprechen und ihn schnell mit der heute 
üblichen Terminologie des Gebiets vertraut machen 


Buchbesprechungen 


wird. Der I. Teil des Büchleins befaßt sich mit 
axiomatischen Ausführungen über den Begriff des 
topologischen Raumes, mit Funktionen, zusammen- 
hängenden und dichten Teilmengen, Basen, Produkt- 
und Quotiententopologien. Der H. Teil geht auf regu- 
läre, normale und kompakte Räume ein. Sein letzter 
Paragraph beschäftigt sich kurz mit lokaler Kompakt- 
heit und Kompaktifizierung. Der III. Teil ist den 
metrischen Räumen und insbesondere den Kompakten 
gewidmet. Besonders hervorzuheben ist hier das Ka- 
pitel über die Metrisierung topologischer Räume, das 
neben den bekannten Sätzen von URYSOHN den vor 
etwa 10 Jahren von Smmrwow, NaGaTA und Bine ent- 
deckten allgemeinen Metrisationssatz enthält. 


Dresden M. LANDSBERG 


Ya. L. Geronimus, Polynomials Orthogonal 
on a Circle and Interval. IX + 210 S. m. 9 Ta- 
feln. Oxford/London/New York/Paris 1960. Pergamon 
Press. Preis geb. 60 s. net 


Das vorliegende Buch beschäftigt sich mit einem 
für Mathematik und Physik bedeutsamen Problem: 
der Untersuchung orthogonaler Funktionssysteme 
und der damit zusammenhängenden Approximations- 
fragen für vorgelegte Funktionen durch solche. Im 
besonderen sind es 3 Fragestellungen, die als die 
wichtigsten herausgestellt und behandelt werden: 

Das asymptotische Verhalten der orthogonalen 
Polynome eines Systems; 

Bedingungen für die Beschränktheit der Funk- 
tionen des Orthogonalsystems; 

Die Angabe der Ordnung für das Wachstum der 
Polynome, falls eine Beschränktheit nicht vorliegt. 

Die Überlegungen werden in der Hauptsache für 
Funktionen entwickelt, die auf dem Einheitskreis 
definiert sind; mit Hilfe einer einfachen Formel 
gelingt die Herstellung des Zusammenhangs zu 
Systemen von Funktionen über einem endlichen 
Intervall der reellen Achse. 

Ein umfangreiches Tafelwerk gibt eine übersicht- 
liche Zusammenstellung der wichtigsten Resultate, 
im besonderen von Abschätzungen und Konvergenz- 
bedingungen. 

Sehr zu begrüßen ist ein Anhang, der die innerhalb 
eines jeden Kapitels benutzten Theoreme aus der 
Analysis beinhaltet; die Lektüre dieses Teiles ist für 
sich schon — losgelöst von den eigentlichen Zielen 
des Buches — sehr wertvoll. 

Das Literaturverzeichnis beschränkt sich auf die 
Angabe der für die Auswahl des Stoffes wichtigsten 
Autoren — man erkennt hierbei den besonderen 
Anteil sowjetischer Mathematiker. 

Das von einem zuständigen Fachmann geschriebene 
Buch, der am Ausbau der Theorie selbst maßgeblichen 
Anteil hat, bringt die modernsten Ergebnisse auf 
diesem Gebiete; an mathematischer Strenge bleibt 
nichts zu wünschen übrig. Die Darstellung des Stoffes 
ist nicht die eines Lehrbuches, sondern viel mehr 
eines Ergebnisberichtes. Alle Überlegungen und Be- 
weise sind wohl vollständig, aber sehr knapp dar- 
gestellt. Daher erfordert die Lektüre ein gründliches 
Eindenken in die Materie und dürfte für einen weniger 
mit den Dingen Vertrauten nicht ganz einfach sein. 

Ein Sachverzeichnis am Ende des Buches würde 
zur Erleichterung der Lektüre beitragen und den 
Nutzen als Nachschlagewerk erhöhen. 


Dresden P. H. MÜLLER 


E. Schörner, Raumbild-Lehrbuch der dar- 
stellenden Geometrie für Ingenieurschulen. 
153 S. m. 60 Anaglyphen und 228 Zeichnungen. 
Tyan 1960. Oldenbourg Verlag. Preis geb. DM 

Im Anfangsunterricht der Darstellenden Geometrie 
treten sehr häufig Schwierigkeiten auf, die dadurch 
bedingt sind, daß dem Studierenden das Sehen und 
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Buchbesprechungen 


Denken im dreidimensionalen Raum noch nicht 


geläufig ist. Das „Lesen“ einer Zeichnung von drei- 
dimensionalen Objekten und umgekehrt, das Ab- 


bilden von Raumobjekten auf eine zweidimensionale 


' Zeichenebene in einer oder in mehreren Ansichten, 


bereitet oft beträchtliche Schwierigkeiten. Anderer- 
seits ist aber für viele technische Zweige ein geschultes 
und geübtes Anschauungsvermögen unumgänglich 
nötig. Um nun die Fähigkeit des räumlichen Denkens 
zu erwecken und zu schulen — ein Hauptziel der 
Darstellenden Geometrie als Unterrichtsfach — wer- 
den häufig räumliche Modelle verwendet. Da beim 
Selbststudium aber selten auf eine Modellsammlung 
zurückgegriffen werden kann, versucht der Autor, die 
Verwendung von Modellen durch Beifügung von 
Anaglyphenbildern zu seinem Raumbild-Lehrbuch 
der Darstellenden Geometrie überflüssig zu machen. 
Tatsächlich erscheinen bei Betrachtung der nach bei- 
gefügter Betrachtungsanleitung waagrecht auf den 
Tisch gebreiteten Plastoskop-Modellen durch die bei- 
liegende Rot-Grün-Brille (Anaglyphenbrille) die dar- 
gestellten Objekte klar und deutlich — fast greifbar 
im Raum — und geben eine bessere Übersicht über 
die räumlichen Zusammenhänge, als sie durch Modelle 
aus Holz, Metall, Gips, Plexiglas usw. gegeben werden 
könnten. Zweifellos ist damit ein Unterrichtsweg 
beschritten, der sehr viel Nutzen bringen wird, doch 
muß noch abzuwarten sein, ob es nicht didaktisch 
richtiger wäre, die Anzahl der gebotenen Raumbilder 
etwas einzuschränken, da diese Fülle an Material und 
die Leichtigkeit der räumlichen Erfassung des Ge- 
botenen den Lernenden zu einer gewissen Trägheit 
in der Bemühung um das räumliche Sehen führen 
könnte. Vielleicht bietet man dem Studierenden mehr, 
wenn man weniger Raumbilder bringt und diese nur 
von wirklich schwer zu erfassenden räumlichen Ob- 
jekten. Diese vielleicht etwas. negative Kritik ist 
aber zugleich ein besonderes Lob über die großartige 
Wirkung derin diesem Werke gebrachten Raumbilder. 


Der dargebotene Lehrstoff entspricht in seinem 
Umfang den Studienplänen für Ingenieurschulen. So 
wird zunächst das Zweitafelsystem (zugeordnete 
Normalrisse) in üblicher Art behandelt, wobei aber 
der Kenner in jedem Kapitel die persönliche Note des 
Autors fühlt. Ein Kapitel über Axonometrie und ein 
kleiner Abschnitt über Perspektive beschließt dieses 
sehr bemerkenswerte Buch, das als Einführung in die 
Darstellende Geometrie sehr zu empfehlen ist. 


Dresden R. BEREIS 


K.-R. Biermann, Vorschläge zur Wahl von 
Mathematikern in die Berliner Akademie. 
Ein Beitrag zur Gelehrten- und Mathematikgeschichte 
des 19. Jahrhunderts. (Abhandlungen der Deutschen 
Akademie der Wissenschaften zu Berlin, Klasse für 
Mathematik, Physik und Technik, Jg. 1960, Nr. 3). 
758. Berlin 1960. Akademie-Verlag. Preis brosch. 
DM 7,50. 


Eine erschöpfende Darstellung der Entwicklung 
der Mathematik im 19. Jahrhundert gibt es bisher 
nicht. Dies hängt zum Teil damit zusammen, daß 
die mühevolle Sichtung und Bearbeitung des um- 
fangreichen Quellenmaterials aus diesem Zeitraum 
noch keineswegs abgeschlossen ist. Es ist nun das 
Verdienst BIERMANNs, durch die Abhandlung ‚‚Vor- 
schläge zur Wahl von Mathematikern in die Berliner 
Akademie“ einen sehr wesentlichen Beitrag für die 
Geschichte der Mathematik der zweiten Hälfte des 
19. Jahrhunderts gegeben zu haben. In der Arbeit 
werden alle fachlich eingehend begründeten Wahl- 
vorschläge für Mathematiker von 1842—1900 an- 
geführt und die für die bedeutendsten Vertreter dieses 
Faches eingebrachten Anträge im vollen Wortlaut 
mitgeteilt. Während man sich früher bei Vorschlägen 
zur Aufnahme in die Akademie meistens mit einigen 
rühmenden allgemeinen Bemerkungen begnügte — 


auch hierfür enthält die Arbeit zwei Proben, die sich 
auf die Wahl DIRICHLETsS und Jacogıs beziehen — 
wurde die Kunst, Wahlvorschläge zu formulieren, 
von WEIERSTRASS und KRONECKER zur Meisterschaft 
ausgebildet. Dabei hatten diese prägnanten Vor- 


schläge nicht nur den Zweck, die Fachgenossen auf 


die wissenschaftliche Bedeutung der Vorgeschlagenen 
aufmerksam zu machen, sondern es lag den Antrag- 
stellern auch daran, die nicht sachverständigen Mit- 
glieder der Akademie über die wissenschaftlichen 
Leistungen der Vorgeschlagenen zu orientieren. BIER- 
MANNs Publikation ist nicht nur für die Geschichte 
dieser Glanzzeit der Mathematik von Bedeutung, 
sondern sie gewährt auch einen Einblick in die Ge- 
lehrtengeschichte des 19. Jahrhunderts. 


Meißen M. MILLER 


R. v. Mises, Mathematical Theory of Com- 
pressible Fluid Flow. Completed by H. Geirin- 
ger and G.S8.S. Ludford. (Volume 3 of Applied 
Mathematics and Mechanics.) XIII + 5148. m. 
174 Abb. New York 1958. Academic Press Inc. 
Preis geb. 15.— $. 


Das Buch, dessen erste drei Kapitel R. vow Mıses 
geschrieben hat, wurde nach seinem Tode an Hand der 
von ihm hinterlassenen Aufzeichnungen von H. GEI- 
RINGER und G.S.S. LUDFORD durch zwei weitere 
Kapitel ergänzt. Es nimmt eine Mittelstellung zwi- 
schen Lehrbuch und Handbuch ein, behandelt die 
theoretischen Grundlagen der Gasdynamik ausführ- 
lich genug, um interessierten Studenten eine Ein- 
führung in dieses reizvolle Gebiet zu geben, und ver- 
steht es ausgezeichnet, den Kenner über den gegen- 
wärtigen Stand. der wichtigsten Probleme der Gas- 
dynamik zu informieren, ohne ihn mit schwierigen 
Konvergenzbeweisen und allzu vielen Einzelheiten zu 
belasten. In über 270 Anmerkungen wird der Leser 
auf wichtige Literaturstellen, abweichende Bezeich- 
nungen und offene Fragen aufmerksam gemacht. 


Besonders hervorgehoben sei im Kapitel I die Auf- 
stellung der “specifying equation”, die im allgemeinen 
Falle eine Beziehung zwischen Druck, Dichte, Ge- 
schwindigkeit und den EurLerschen Variablen dar- 
stellt und aus den thermodynamischen Grundgesetzen 
hergeleitet wird ; sie ergänzt die Bewegungsgleichungen 
und die Kontinuitätsgleichung des Kontinuums zu 
einem System von fünf Differentialgleichungen zur 
Ermittlung von Druck, Dichte und Geschwindigkeit. 
Im Kapitel II wird die Charakteristikentheorie von 
Grund auf entwickelt. Im Kapitel III wird nach der 
Herleitung der Stoßbedingungen für die eindimensio- 
nale instationäre Strömung einer zähen Flüssigkeit 
ausdrücklich darauf hingewiesen, daß diese stoß- 
behafteten Strömungen nicht als ‚„unstetige‘“ Lö- 
sungen der Bewegungsgleichungen einer idealen 
Flüssigkeit aufzufassen sind, sondern allenfalls als 
asymptotische Lösungen im Falle verschwindender 
Zähigkeit. Im Kapitel IV werden die Singularitäten 
der Hodographentransformation, die Grenzlinien und 
die Verzweigungslinien, ausführlich diskutiert und 
eine geometrische Charakterisierung der Schallgrenz- 
linien nach LUDFORD und ScHor bewiesen. Im Ka- 
pitel V werden im Anschluß an das Korrespondenz- 
prinzip von TscHArLIGINn zwischen inkompressibler 
und kompressibler Strömung in sehr übersichtlicher 
Weise die Methoden von LIGHTHILL und BERGMAN zur 
Berechnung ebener stationärer Unterschallströmung 
dargestellt. Nach Herleitung der Stoßbedingungen für 
die ebene stationäre Strömung einer zähen Flüssigkeit 
werden geschickt ausgewählte Beispiele durchgerech- 
net und mit Hilfe sorgfältig gezeichneter Figuren er- 
läutert. Ferner wird die Differentialgleichung für die 
Stromfunktion der stationären nichtisentropen Strö- 
mung einer zähen Flüssigkeit hergeleitet, die nicht 
vollkommen elastisch zu sein braucht. Das Kapitel 
endet mit einem kritischen Bericht über die Proble- 
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matik der transonie flows, mit der sich voX Mıs£s in 
seinen letzten Lebensjahren eingehend beschäftigt hat. 


neuer Versuchsergebnisse. 


Halle/Saale H. SCHUBERT 


- R. Bellman, Introduction to Matrix Analysis. 
XX + 3288.” London 1960. MeGraw-Hill. Preis 
geb. 77s. 6d. 

Die vorliegende Einführung in die Matrizentheorie 
weicht in mehrfacher Hinsicht von anderen Dar- 
stellungen ab. In seinem — übrigens sehr lesens- 
werten — Vorwort stellt sich der Verfasser auf den 
Standpunkt, daß mathematische Lehrbücher den 
Leser leider viel zu oft mit einer solchen Menge 
Stoffes überfluten, daß ihm dadurch der schnelle Zu- 
gang zu interessanten und gehaltvolleren Fragen, wie 
sie vor allem auch in den Anwendungsgebieten auf- 
treten, verbaut wird. Er zieht daraus die Konsequenz 
und beschränkt sich auf einige wenige Fragestellungen 
und verweist viele, auch wichtige Tatsachen der 
Matrizentheorie, die sonst verhältnismäßig breit dar- 
gestellt zu werden pflegen, in die zur Ergänzung bei- 
gefügten Übungsaufgaben. Den Zugang zu der Theo- 
rie der symmetrischen Matrizen, die einen wesent- 
lichen Teil des Buches ausmacht, findet der Leser von 
dem Extremwertproblem quadratischer Formen her. 
Es folgen Kapitel über die grundlegenden Begriffe 
und Regeln des Matrizenkalküls und über die Haupt- 
achsentransformation symmetrischer Matrizen (so- 
wohl für den Fall getrennter Eigenwerte als auch für 
den allgemeinen Fall). Das fünfte Kapitel ist dem 
Extremwertproblem bei quadratischen Formen für 
den Fall, daß noch lineare Nebenbedingungen vor- 
liegen, gewidmet. Es folgen Ausführungen über die 
einfachsten Matrizenfunktionen, über die Extremal- 
eigenschaften der Eigenwerte, über wichtige Unglei- 
chungen und schließlich über die als dynamisches 
Programmieren bezeichnete Methode zur Lösung 
gewisser Minimierungsaufgaben. Die Kapitel 10 
bis 12, bei denen die Symmetrie der Matrizen nicht 
mehr vorausgesetzt wird, knüpfen an die Probleme 
an, die bei linearen Differentialgleichungssystemen 
mit konstanten und variablen Koeffizienten auf- 
tauchen, und dringen bis zur Stabilitätstheorie vor. 
Es schließen sich Ausführungen über MARKOFF- 
Matrizen und ihre Bedeutung in der Wahrscheinlich- 
keitstheorie und über stochastische Matrizen an. Das 
letzte Kapitel befaßt sich mit positiven Matrizen 
(d.h. Matrizen mit nur positiven Elementen) und 
mit ihrer Anwendung bei der Lösung ökonomischer 
Probleme. 

Von dem eigentlichen Charakter des Buches kann 
diese Übersicht keinen vollständigen Eindruck ver- 
mitteln. Als Lehrbuch kann man es nur in einge- 
schränktem Sinne ansprechen. Durch die überaus 
zahlreichen, nicht nur mit ihren Titeln angeführten, 
sondern auch kommentierten Hinweise auf das 
weiterführende Schrifttum erhält es fast einen enzy- 
klopädischen Charakter. Daß jedem Kapitel viele 
Übungsaufgaben beigefügt sind, wurde schon hervor- 
gehoben. Sie sind von recht verschiedenem Schwierig- 
keitsgrad, sind oft als Ergänzung des im Text be- 
handelten Stoffes anzusehen und geben teilweise sogar 
in nur wenigen Zeilen Probleme wieder, denen in der 
Literatur bedeutsame Originalarbeiten gewidmet sind. 
In seinem Vorwort begründet der Verfasser dieses 
Vorgehen. Er kündigt im übrigen an, daß dem vor- 
liegenden Werk als erstem Band einer umfassenden 
Serie weitere Bände anderer Verfasser folgen werden, 
die sich eingehend mit den Spezialproblemen zu be- 
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Dresden INEICH 


L. I. Priwalow, Einführung in die Fonic 
theorie, Teil III. V + 1888. m. 57 Abb. Leipzig 
lagsgesellschaft. 


1959. B.G. Teubner Ver Preis geb. 
12,90 DM. Sr 
hält einige der Gebiete, die in 

Fortsetz gebracht zu werden pflegen. 

— Im ersten Teil (elliptische ) werden in 


iell wird daneben die Abbildung 
von Gebieten, die von Kurven zweiter Ordnung be- 


sind, und von Polygonen behandelt. “ 
den schlichten Funktionen gewidmeten Teil, finden 
sich u.a. Verzerrungssätze und das Wichti über 
das Koeffizientenproblem bei der Entwi ein- 


blättriger Funktionen. — Auch bei diesem \ 
bieten Darstellung und die Üb die 
Gewähr für erfolgreiche Einarbeitung in die genannten 
Gebiete. 

Gießen 


K. MaRrtTHN 


6. Julia, El&ments d’Algebre. VII + 2307S. 
Paris 1959. Gauthier-Villars. Preis brosch. $ 7.93. 


Das vorliegende Buch gibt in sehr klarer, knapper 
Darstellung den Stoff von Algebra-Vorlesungen wie- 
der, die der bekannte Verfasser seit 1956 vor Stu- 
denten des ersten Studienjahres der Ecole Polytech- 
nique in Paris gehalten hat. Inhalt und Darstellung 
des Buches tragen der modernen Entwicklung der 
Mathematik Rechnung, die gezeigt hat, daß die An- 
wendung der fundamentalen Begriffe der Algebra — 
Gruppe, Ring, Körper — auch für die Geometrie und 
Analysis Bedeutung gewonnen hat. 

Die beiden ersten Kapitel vermitteln die für das 
weitere erforderlichen algebraischen und analytischen 
Kenntnisse. Die oben angegebenen Strukturen 
Gruppe, Ring und Körper werden definiert und 
einige erste Aussagen über sie hergeleitet; weiter 
sind hier die Elemente aus der Algebra der linearen 
Operatoren zusammengestellt (Matrizenrechnung, 
Koordinatentransformation, Eigenwerte, Spektrum). 
Den Hauptteil des Buches bilden das dritte und 
vierte Kapitel. Das dritte behandelt die Metrik im 
n-dimensionalen Vektorraum, im vierten Kapitel ö 
werden die hermiteschen und unitären Operatoren 
und Matrizen untersucht. Das abschließende fünfte 
Kapitel ist der multilinearen Algebra gewidmet, wobei 
ihrer Wichtigkeit in den Anwendungen wegen vor 
allem die Tensoren 2. Stufe in orthonormierten 
Koordinaten betrachtet werden. 


| 
| 


Dresden M. Hasse 


6. Sansone — J.Gerretsen, Lectures on the 
Theory of Functions of a Complex Variable, 
Vol.1. 488S. Groningen 1960. P. Noordhoff Ltd. 
Preis brosch. $ 12.—. 

Das vorliegende Werk hebt sich in Aufbau und 
Darstellungsweise eindrucksvoll von der Mehrzahl der 
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z. Zt. vorhandenen Lehrbücher der Funktionentheorie 


' ab; es ist insbesondere keine übersetzte Neuauflage 


des 1947 erschienenen Saywsox&eschen Lehrbuches. 


Nach einem einleitenden Kapitel, das neben den 
Grundbegriffen sich gründlich mit Potenzreihen, den 
elementaren Funktionen (Exponentialfunktion, Kreis- 
funktionen) und deren Umkehrungen beschäftigt, 
behandelt das zweite Kapitel hauptsächlich den 
Gedankenkreis des CaucHyschen Integralsatzes. Die 
Verff. schließen sich dabei wesentlich an die Artıssche 
Notre-Dame-Vorlesung an und schaffen die topolo- 
gischen Grundlagen mit Hilfe des Begriffes der Win- 
dungszahl. Das dritte Kapitel behandelt die iso- 
lierten Singularitäten, die (ebenfalls nach Arrıs) ohne 
Verwendung der LAuRENT-Entwicklung eingeführt 
werden, das vierte bringt den WEIERSTRASSsschen 
Produktsatz und die Partialbruchentwicklung nach 
CAucHY und MITTAG-LEFFLER, das fünfte ist dann 
im üblichen Umfang den elliptischen Funktionen 
gewidmet. Den Abschluß bilden drei Kapitel über 
ganze Funktionen (BoREL, PoINcARE), über DIRIcH- 
LETsche Reihen (einschließlich RısEmansscher Zeta- 
funktion und LArLAck-Integral) und über Summierung 
von Potenzreihen außerhalb des Konvergenzkreises. 


Die geometrische Betrachtungsweise der Funk- 
tionentheorie, insbesondere also die konforme Ab- 
bildung, ist dem zweiten Band des Werkes vorbe- 
halten. 


Das Buch setzt nicht mehr als die Kenntnis der 
elementaren Analysis voraus. Man darf sagen, daß 
die Absicht der Verff., ‚einige Meisterstücke mathe- 
matischen Denkens darzubieten und einem weiten 
Leserkreise zugänglich zu machen“, dank der vor- 
züglichen klaren Darstellungsweise in vollem Umfang 
erreicht wurde. 


Gießen K. MARUHN 


6. Hellwig, Partielle Differentialgleichun- 
gen. 2468. m. 35 Abb. Stuttgart 1960. B. G. Teub- 
ner Verlagsgesellschaft. Preis geb. DM 29,80. 


Das vorliegende, als eine Einführung bezeichnete 
Werk soll der erfolgreichen Einarbeitung in die Theo- 
rie der partiellen Differentialgleichungen für Leser mit 
etwa Vordiplomkenntnissen dienen. Daß es dem Ver- 
fasser trotzdem darüber hinaus gelingt, einige der 
wichtigsten modernen Methoden und auch einige neue 
Ergebnisse zu bieten, muß als wesentliches Merkmal 
dieses ausgezeichneten Werkes herausgestellt werden. 


Verf. verzichtet bewußt auf die Behandlung von 
Gleichungen erster Ordnung in einer gesuchten 
Funktion (sie läßt sich auf die der gewöhnlichen 
Differentialgleichungen zurückführen) und beginnt 
mit den Beispielen der Wellen-, Potential- und Wär- 
meleitungsgleichung, an denen die unterschiedliche 
Problematik bei Anfangs- und Randwertaufgaben 
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herausgearbeitet wird. Es folgen im zweiten Kapitel 
die Typeneinteilung (auf ihre Wichtigkeit wird aus- 
führlich eingegangen), Normalformen und charak- 
teristische Mannigfaltigkeiten bei Gleichungen zweiter 
Ordnung und Systemen erster Ordnung. Im dritten 
Kapitel finden sich für die verschiedenen Typen und 
Aufgaben Eindeutigkeitssätze (auch beim gemischten 
Typ) auf der Grundlage des Maximum-Minimumprin- 
zips sowie der Energieintegralmethode. Das vierte 
Kapitel bildet dann mit der Behandlung von Existenz- 
fragen den Höhepunkt des Buches. Die hier ver- 
wendeten Methoden sind einmal, je nach der Aufgabe, 
im hyperbolischen und parabolischen Fall die suk- 
zessiven Näherungen und die LarLAck-Transfor- 
mation, bei den Randwertaufgaben für die elliptischen 
Gleichungen die Theorie der schwachen Lösungen 
unter Heranziehung des Weyrschen Lemmas in einer 
erweiterten Fassung. Das letzte Kapitel behandelt 
Existenzfragen beielliptischen Differentialgleichungen 
unter Verwendung funktionalanalytischer Hilfsmittel, 
die kurz zusammengestellt sind. Es finden sich hier 
die SCHAuDERsche Beweistechnik (skizziert), Eigen- 
wertprobleme und Randwertprobleme bei elliptischen 
Systemen erster Ordnung in zwei gesuchten Funk- 
tionen. — Die eingestreuten, nach ihrer Schwere 
unterschiedenen Aufgaben unterstützen in glücklicher 
Weise die Einarbeit in dieses gerade auch für den 
Leser der ZAMM wichtige Gebiet. 


Gießen K. MARUTHN 


S. Dans (Ass. Prof. of Economics, University. of 
Copenhagen), Linear Programmingin Industry. 
Theory and Applications. VIII + 1208. m. 6 Abb. 
Wien 1960. Springer-Verlag. Preis brosch. DM 20,—. 


Das vorliegende Buch ist — seiner Anlage ge- 
mäß — eine hauptsächlich für Nichtmathematiker, 
für Ökonomen und Technologen bestimmte Schrift 
über lineare Optimierung. 

Aus dem Inhalt: Nach Definition einer linearen 
Optimierungsaufgabe und einer Erläuterung wichtiger 
Sätze und Begriffe folgen zahlreiche ökonomische Bei- 
spiele, nämlich Produktionsaufgaben bei linearer 
Technologie, Mischungsaufgaben, Kapazitäts-, Lager- 
haltungs- und Transportprobleme. Anschließend 
werden die Simplex- und duale Simplex-Methode in 
ihrer rechnerischen Form besprochen. Ergänzt wer- 
den diese Lösungsmethoden noch durch eine Be- 
schreibung der M-Methode (zur Aufsuchung einer 
ersten zulässigen Basislösung), des Entartungsfalles, 
der Besonderheiten bei der Lösung des Transport- 
problems und des parametrischen Programmierens. 
Im Anhang werden schließlich die Beweise für das 
Simplex-Theorem, das Simplex-Kriterium und das 
Dualitätstheorem nachgeholt. 


Dresden L. BITTXNER 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 
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